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ABSTRACT 


The study defined under Contract NAS8- 28358 consisted of four parallel 
efforts: (1) Modal analyses of elastic continua for Liapunov stability analysis 

of flexible spacecraft; (2) Development of general-purpose simulation equations 
for arbitrary spacecraft; (3) Evaluation of alternative mathematical models 
for elastic components of spacecraft; and (4) Examination of the influence of 
vehicle flexibility on spacecraft attitude control system performance. 

This report includes a complete record of achievements under tasks (1) 
and (3) above, in the form of technical appendices, and a summary description 
of progress to date under tasks two and four. 

Task (1) has provided the basis for the Ph. D. dissertation of Andre Colin 
(see Appendix 3, in Volume 2 of this report). This task in itself required two 
phases of investigation: modal analysis and stability analysis. The modal 
analysis is accomplished for a range of continuum models (strings, beams 
and thin plates with various boundary conditions on spinning spacecraft) by 
means of singular perturbation methods, and the stability analysis is accom- 
plished by using Liapunov theorems with the momentum- constrained Hamiltonian 
as the testing function. 

Task (2) is the basis for the Ph. D. dissertation of Arthur S. Hopkins, 
which is still in progress. 

Task (3) is the subject of two technical papers by the Principal Investiga- 
tor, included here as Appendices 1 and 2. In these papers the range of appli- 
cability of various discrete and continuous models of nonrigid spacecraft is 
examined. It is concluded that there is a domain of engineering applicability 
for each of the models considered, but that finite elements models are 
generally the most valuable for flexible spacecraft simulations. 


iii 



Task (4) is currently receiving primary attention by the Principal 
Investigator and a postdoctoral scholar. Dr. Yoshiaki Ohkami. Results 
will be described in forthcoming documents. 



UNIVERSITY OF CALIFORNIA 
Los Angeles 


MODAL ANALYSIS FOR LIAPUNOV STABILITY 
OF ROTATING ELASTIC BODIES 


A dissertation submitted in partial satisfaction of the 

requirements for the degree Doctor of Philosophy 

/ 

in Engineering 
by 


Andre Daniel Colin 



ACKNOWLEDGMENT 


First, I would like to thank Dr. Willems of the University of 
Louvain who introduced me to space vehicle problems. 

I also wish to express my deepest gratitude to my advisor, 

Dr. Likins, for his guidance throughout the development of this 

v 

dissertation, for his collaboration and his availability. I am also 
grateful to Drs. Lange, Leondes, Cole and Eisenberg for serving on 
my doctoral committee. 

Finally, I would like to thank the European Space Research 
Organization and the National Aeronautics and Space Administration for 
their financial support and the opportunity they gave me to stay in 
this country. 



TABLE OF CONTENTS 


Page 

NOMENCLATURE ..... iv 

LIST OF FIGURES ix 

ACKNOWLEDGMENT x 

VITA . . . . xi 

ABSTRACT xii 

CHAPTER 1 - STATEMENT OF THE PROBLEM 1 

CHAPTER 2 - DYNAMICS OF ROTATING ELASTIC BEAMS AND CABLES. . . 11 

2.1 Introduction 11 

2.2 Dynamics of a Rotating Elastic Beam 11 

2.3 Dynamics of a Taut, Rotating, Elastic Cable. ... 43 

CHAPTER 3 - DYNAMICS OF ROTATING THIN PLATES AND MEMBRANES . . 71 

3.1 Introduction 71 

3.2 Dynamics of a Taut Rotating Membrane in a 

Circular Ring. 71 

3.3 Dynamics of a Rotating, Free, Elastic, Circular 

Membrane ....... 101 

CHAPTER 4 - LIAPUNOV STABILITY ANALYSIS 131 

4.1 Introduction 131 

4.2 Concept of Stability 131 

4.3 Stability Analysis 135 

CHAPTER 5 - CONCLUSIONS AND RECOMMENDATIONS. . 175 

REFERENCES .......... ..... 179 

APPENDIX I - ORTHOGONALITY OF THE MODES ... . . 181 

APPENDIX II - SOLUTION TO- THE REDUCED 'EQUATION OF THE 

TAUT MEMBRANE 185 

APPENDIX III - BEHAVIOR NEAR THE ORIGIN. 191 

APPENDIX IV - DETERMINATION OF THE CONSTANT ^ 195 

APPENDIX V - EQUATIONS FOR THE ROTATIONAL MOTION 197 


2-5 



LIST OF FIGURES 

Page 

Figure 1. Rotating Uniform Elastic Beam . . 12 

Figure 2. Correspondence with Classical Results. ...... 42 

Figure 3. Rotating Uniform Taut Cable. ........... 44 

Figure 4. Eigenfrequencies Corresponding to First Mode . . . 52 

Figure 5. Mode Shapes Corresponding to First Mode. 52 

Figure 6. Eigenfrequencies Corresponding to the Second 

Mode ... . 53 

Figure 7. Mode Shapes Corresponding to the Second Mode ... 54 

Figure 8. Rotating Taut Membrane 72 

Figure 9. Distribution of Stresses 74 

Figure 10. Rotating Free Membrane 101 

Figure 11. Kirchoff Condition 103 

Figure 12. Freely Spinning Spacecraft 136 

Figure 13. Restricted Model ............ 150 

Figure 14. Particular Case. . . 164 

Figure 15. Simple Example . 171 



NOMENCLATURE 


P( ) 
P 

a 

(•) 

(') 


y »s ,x,r 
h i ( ) 

g ± ( ) 

£ i< > 


y.y 

_ * 
x,x 

_ * 
r,r 


y ,x .r 

J V n t ) 


transverse displacement 
external axial load 
harmonic function of the vibration 
derivative with respect to time 

spatial derivative with respect to dimensionless variables 
dimensionless variables 

terms of the expansion of the modes away from the boundaries 

terms of the boundary layer expansions 

terms of the boundary layer expansions 

boundary layer coordinates 

boundary layer coordinates 

boundary layer coordinates 

fixed quantity introduced in the matching process 


u 

h 

m 



<£ 

2 


Ai, Bi 


radial displacement 
thickness of the disk 

coefficient of the argument in the circular dependence of 
the modes 

Legendre functions of the first and second kind, respectively 
Airy integrals 


T 

T 

0 

k 

TST 

El 


built-in tension 
minimum built-in tension 
defined by T = 
transcendentally small terms 
flexural stiffness of the beam 


2-9 



NOMENCLATURE (Cont'd) 


€,U,C 

P.0 

V 

<t> 


a 


E modulus of elasticity in the modal analysis and identity 

matrix in the stability analysis 

£ perturbation parameter 

v^(e),K^(e) asymptotic sequences 

. y (e) *6.(e) asymptotic sequences 

y linear density of the beam or mass per unit volume for the 

membrane 

reference axis for the description of the transverse 
displacement 

polar coordinates of an element dm of the membrane 
Poisson modulus 
mode shapes for the beam 

radial dependence of the membrane mode shapes 

eigenfrequency of the beam vibrations 

eigenfrequency of the membrane vibrations 
nominal spin rate 

eigenfrequency of the beam in the transformed equation 
terms of the expansion of the eigenf requence 
defined by A^ = k^y^ 

eigenfrequency of the membrane in the transformed equation 
tangential and radial strains, respectively 
tangential and radial stresses, respectively 
Laplacian operator 


<{> 


a,m 


0) 


a 


0) 


a,m 


ft 


A (e) 

a 

A. 

■ i 


A (e) 

a,m 

e Q £ 

0 p 

a Q a 

6 p 


V 


£_ 1 _9 _1 

3p 2 p p 2 ,30 2 


m 


Linear differential operator 


dr 


+ 1 X- 

2 r dr 


m 


2-10 



NOMENCLATURE (Cont’d) 


Y 


1 

2 

C 2 

2 


defined by y = 1 + m 
defined by k 2 = 


l+2k 2 +v 


, , , ,2 l+3v 

defined by k^ = 


l+2k 2 +V 


defined by k 2 = k 2 k 2 


V 


F (a, 3, y;x) hyper geometric function 


CM 

N 


Body B 
Body A 
T 
H 
H 

c 

V 

U d 

K 

K 


2 . 

{b} 

{n} 


center of mass 

✓ , - 

location of system CM when steadily spinning 
portion of system identified as rigid (core) 
portion of system identified as flexible (appendages) 
kinetic energy of complete system 
Hamiltonian 

constrained Hamiltonian 
potential energy 
potential energy of deformation 
general spinning stiffness matrix 
spinning stiffness matrix for the planar model 
mass of complete system 
inertial generic position vector 

A 

vector basis fixed in Body B identified as (b^, 
inertial vector basis 
vector from CM to N 

deformation vector and its representation in {b} 


2-11 


|o'> 



NOMENCLATURE (Cont'd) 


angular momentum vector, its representation in (b) and its 
magnitude 

body generic position vector 


orthogonal transformation relating (b) to {n} 

Euler (attitude) angles 

inertial angular velocity vector and its representation in 
vector basis '{b} 

position vector of element dm subsequent to spin and its 
representation in {b} 

inertia matrix about point N of the complete system at 
steady state 

inertia dyadic of complete system about point N and its 
inertia matrix representation is (b) 

inertia dyadic of appendage about point N and its inertia 
matrix representation in {b} 

inertia dyadic of core about point N and its inertia matrix 
representation in {b} 

N 

diagonal elements of Iq 

inertia matrix about point N consisting of first order 
appendage terms 

inertia matrix about point N consisting of second order 
appendage terms 

spin rate at steady state 

diagonal matrix made of the squares of the modal frequencies 
generalized coordinates of deformation 

mode shapes representation in the stability analysis 

column matrix of the modal coordinates 
column matrix of the mode shapes 



M 


m 


A. 


A, 


n. 


NOMENCLATURE (Cont’d) 
symmetric matrix defined by 


symmetrix matrix defined by 
stability analysis 


T 

<t><P dm 


I lp/ 

a column matrix defined by / dm 

A 


T 

<p dm 


in the 


a column matrix defined by J dm 

A 


symmetrix matrix defined by A^A^ 


IL 


symmetrix matrix defined by A_A 


2 2 



2-13 



Page Intentionally Left Blank 



VITA 


Born , 

1966 — Candida! Ingenieur Civil, University of Louvain 
1969 — Ingenieur Civil Mecanicien, University of Louvain 

1969- 1970 — Member of the Scientific Personnel of the University of 

Louvain — Institute of Mechanics and Applied Mathematics, 
Heverlee, Belgium 

1970- 1972 — European Space Research Organization and National 

Aeronautics and Space Administration International Fellow 
at the University of California, Los Angeles. 


1972 — Research Assistant, School of Engineering and Applied Science, 
University of California, Los Angeles, California. 


Page Intentionally Left Blank 



ABSTRACT OF THE DISSERTATION 

Modal Analysis for Liapunov Stability 
of Rotating Elastic Bodies 
by 

Andre Daniel Colin 
Doctor of Philosophy in Engineering 
University of California, Los Angeles, 1973 
Professor Peter W. Likins, Chairman 

The determination of the attitude stability of spinning 
elastic bodies is considered here. A complete study of this problem 
requires in general two parallel investigations. In the first part, a 
modal analysis is accomplished for a class of special problems and 
finally, a derivation for literal attitude stability criteria for 
idealized spinning flexible spacecraft is performed. 

The modal analysis is accomplished for a class of struc- 
tures idealized by beams or disks of small flexural rigidity. In the 
description of the flexible appendage, a continuous elastic model is 
used, and the vibrations in the direction of the spin axis are analyzed 
for a fully constrained base rotation. The study is specialized to 
the cases where the high spin rate is combined to the small flexural 
rigidity to produce a small parameter £. The method of matched asymp- 
totic expansions represents the general framework used in obtaining 
the solution for the several structures studied. The determination of 
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the eigenfrequencies and free vibration mode shapes is obtained by the 
concurrent use of the orthogonality of the mode shapes and the con- 
ditions derived from the matching process. 

Also, the stability of spinning flexible satellites in a force- 
free environment is analyzed. The satellite is modeled as a rigid core 
having attached to it a flexible appendage described by a continuous 
elastic model, such as those mentioned earlier. A Liapunov stability 
procedure is used where the Hamiltonian of the system, constrained 
through the angular momentum integral so as to admit complete damping, 
is used as a testing function. Closed form stability criteria are 
generated for the first mode of a restricted appendage modal lying in 
a plane containing the system center of mass and orthogonal to the 
spin axis. 

/ i 
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CHAPTER 1 


STATEMENT OF THE PROBLEM 

The influence of flexibility on spinning structure has been the 
subject of numerous technical investigations since the unexpected 
behavior of Explorer I. It was, indeed, the deformation of flexible 
antennas that caused the spacecraft to deviate from its prescribed 
motion. The variety of tasks fulfilled by modern satellites implies 
the presence of large appendages whose weights, for obvious reasons, 
are kept as low as possible. As the complexity and the number of 
spacecraft orbiting in space increase, the previous concern becomes 
more and more important. 

In the development of dynamic models, describing the motion of 
flexible spacecraft, many authors employ a system of hybrid coordinates 
wherein discrete coordinates (for the translation and rotation of 
rigid bodies or reference frames) are used together with distributed 
or modal coordinates (for the deformations of elastic bodies) . In the 
modeling of the vehicle appendages , various types of idealizations 
have been adopted: the elastic continuum model, the distributed-mass 

finite element model, and the elastically interconnected nodal body 
model. Common to all types of modeling, the modal analysis requires 
the derivation of the linearized equations of small oscillations from 
the constant state of deformation, induced by spin, and the transforma- 
tion of these equations into a system of uncoupled equations of motion 
in terms of "normal mode coordinates." They represent an infinite 
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system for a continuous model, and a finite set for a discretized model, 
but in either case consideration of only a reduced number of them is 
a practical necessity. 

It is mainly towards the first type of modeling — the elastic 
continuum model — that this dissertation is oriented, and towards the 
derivation of a stability analysis when such a model is used. The use 
of the elastic continuum model presents with respect to the other types 
of modeling a noticeable advantage when applied to some specific 
structures. It is, on the other hand, more difficult to implement for 
a general structure. Our concern in this work is to deal with some of 
those applications where the continuum model prevails. 

The use of elastic continuum model for a rotating structure has 
generally been developed in the context of radial or axial beams . The 
investigation of continuous models for rotating planar structures such 
as membranes is more recent, and has generally been affected by mathe- 
matical difficulties arising in the development of the modes. 

In the development of modal analysis of a rotating radial beam, 

we should mention the work done by R.T. Yntema.^ The author used in 

his analysis of the modes of radial beams a Galerkin method, where he 

expands the modes in terms of the nonrotating beam modes. Interested 

2 

in the same problem, J.E. Rakowski and M.L. Renard, and P.C. Hughes 
3 

and J.C. Fung, used similar procedures to obtain the eigenf requencies 
of a rotating radial beam. They used the linear property of the 
equations of the deformations in order to generate a family of solu- 
tions. The satisfaction of the boundary conditions of the problem gave 
them a procedure to converge to the exact value of the eigenf requencies . 
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The generation of modes for circular membranes using a continuous 

model has mainly been developed for membranes of zero stiffness. Most 

probably, the earliest work was done by H. Lamb and R.V. Southwell,^ 

in 1921, where the authors derived the modes for a spinning disk of no 

flexural rigidity. They, also, used an approximate method to bound 

the lowest frequency in the case where both bending and membrane 

effects are important. Among the other authors, whose contributions 

5 6 

should be mentioned are J.G. Simmonds, and W. Eversman, Both of 
them spent a great deal of effort to solve the problem of a spinning 
membrane when it is clamped at its center. Finally, W. Eversman and 
R.O. Dodson^ studied in 1969 the free vibrations of a centrally 
clamped spinning circular disk where they introduce flexural rigidity 
into their analysis. Here, again, the linear property of the problem 
is used in deriving the general solution. All of the above references 
dealing with beams or membranes of nonzero stiffness employ procedures 
requiring extensive numerical computations. 

In contrast to those past achievements, the problem of the elastic 
continuum model has been approached here along directions guided by 
future applications. It is clear that flexible spinning spacecraft 
will continue to be designed and flown in the future. It is expected 
that many of these satellites will exhibit large flexible appendages 
such as antenna arrays or solar panels. Many of these applications 
are conceived on the verge of instability, and it becomes more and 
more important to develop proper stability criteria to cope with these 
future applications. An important class includes flexible spacecraft 
with high spin and low flexural rigidity. The partial differential 
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equations describing the motion of flexible appendages for the previous 
class of problems are characterized by a small parameter appearing in 
the coefficient of the highest derivative. This common feature will 
constitute the basis for our analytical approach to modal analysis. 

The resulting modal coordinates will be used in the determination of 
attitude stability criteria by means of Liapunov's direct method, 
employing as testing function the Hamiltonian constrained by the angu- 
lar momentum integral. 

8 

The method of matched asymptotic expansions represents the 
general framework and the main mathematical tool used ih obtaining 
modal coordinates. Briefly, the method consists of separating the 
region of interest into "boundary-layer" regions near the boundaries 
where flexural rigidity effects are important and a central region in 
which the terms expressing the centripetal acceleration dominate. The 
solution is then developed in asymptotic expansions in terms of 
the small parameter observed in the highest derivative term. The 
approximation to the solution is then obtained by truncating the 
expansion to a finite number of terms , the error being small for 
sufficiently small values of the parameter. The use of asymptotic 
expansions for the description of the modes of rotating structures 
is not new, and we should mention two references using the previous 
techniques to study the transverse vibrations of rotating structures. 

The first reference to survey is the work done by W.E. Boyce and 

9 

G.H. Handelman in 1961, where they approached the problem of the 
rotating beam with tip mass. This paper is mainly conceived as an 
application of a previous paper published by J. Moser, ^ wherein 
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uniformly valid asymptotic expansions are derived for an equation of 
the type treated there. In this last paper Moser derives a set of 
linearly independent solutions for the problem and, then, uses the 
linear property of the problem to come up with a general solution. 

The linearly independent solutions are taken in the form of 
B(x,r)) exp [ri ^h(x)] where ri is of the order of magnitude of the small 
perturbation parameter of the problem and B(x) and h(x) represent 
functions to be determined. Applying this procedure to the rotating 
beam with tip mass, Boyce and Handelman were mainly interested in what 
is called the zero*”* 1 order solution — the first term of the expansion. 
The solution to this last problem is then obtained through an energy 
method after the observation that one of the boundary conditions in the 
problem represents a natural boundary condition. What is meant by 
natural boundary condition is that the minimization problem arising in 
the Rayleigh-Ritz method gives us as transversality condition for the 
tip mass, the boundary condition of the zero t ^ 1 order problem. In 
Reference 9, the authors do not proceed to the next term of their 
asymptotic expansion. 

Another relevant reference using the method of asymptotic expan- 
sion for rotating structures is given by J.H. Abel and W.C. Kerr,^ 

In their paper, they applied the technique of asymptotic expansion to 

a rotating cable-counterweighted space station in orbit. The problem 

12 

without any flexural rigidity was already solved by V. Chobotov 
in 1963, and represents in fact the zero t ^ 1 order solution for their 
problem. In contrast to the method used In Reference 9, it is by 
matching the central and boundary layer solutions that they were able 
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to come up with the sufficient number of conditions for the determin- 
ation of all unknown constants introduced by the integration of the 
several equations. 

It is by using the conditions coming from the matching of the 
solution valid in the boundary-layer and the solution valid in the 
central region, and also by considering the orthogonality relationship 
between the eigenfunctions, that we plan to obtain explicit expressions 
for the eigenvalues and free vibration mode shapes for several 
elementary structures. For all the structures studied, we concentrate 
on the "transverse vibrations," involving oscillatory motions which 
are parallel to the nominal spin axis and transverse to a plane 
established by the structure, since these vibrations are most critical 
for the stability of a spinning structure. We also make throughout 
the study the assumption that the motion of the rigid rotating base to 
which the elastic appendage is attached is not affected by the trans- 
verse vibrations of the flexible appendage. For a free spinning 
spacecraft, the center of mass of the entire system remains at rest in 
inertial space, but it does move within the core for some modes. 
Similarly, the core body rotates for some of the system normal modes, 
but our modal analysis is based on a fully constrained base rotation. 
This represents in fact a convenient assumption — which is necessary 
if our answers are to have value for hybrid coordinate analysis of 
appendages on arbitrary spinning bodies — and permits the obtention 
of uncoupled modes — conditions under which we can transform the 
partial differential equations into a set of ordinary differential 
equations . 
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In Chapter 2, the derivation of the eigenmodes and eigenf requencies 
for two unidimensional types of flexible appendages is presented. At 
first, we looked at the problem of the radial uniform classical beam. 

The classical assumptions consisting of neglecting the shear deforma- 
tions and the rotatory inertia were made. We also limited our study 
to a radial beam clamped at its root on the spin axis. After the 
completion of this first case — where the spinning motion acts as an 
element stiffening the structure — we looked at the case of a uniform 
cable clamped at both extremities, and spinning about a central axis. 

The spinning motion introduces in this structure a softening effect 
and requires the introduction of a built-in tension. This configura- 
tion seems to be of little direct application, but it has been examined 
in order to familiarize ourselves with the next problem of wider 
application, discussed in the following chapter. 

In Chapter 3, the modal analysis of a rotating membrane is con- 
sidered. Two cases are investigated. In the first one, the flexible 
appendage consists of a circular membrane clamped along its edge, and 
spinning about a central axis normal to its plane. Here, again, the 
rotation introduces a decrease in the eigenf requency of the nonrotating 
structure and justifies the introduction of a built-in tension. This 
last effect is generally referred as the effect of preload and is 
known to be highly configuration-dependent. Here, again, the assump- 
tion is made that the motion of the spinning rigid rim is uninfluenced 
by membrane vibrations. Finally, ,the freely spinning membrane is 
analyzed and it is observed that for this last structure, the eigen- 
frequencies of the rotating membrane are less affected by the flexural 
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rigidity than in the case with outer rim constrained by a spinning 
rigid ring. 

Finally, in Chapter 4, a general derivation is given for the 
stability analysis of rotating structures when the deformations of the 
flexible parts are expressed in terms of the modes of the rotating 
structure when determined from a continuous elastic model. Also, some 
of the restrictions introduced in our derivation are justified. It is 
shown, for instance, that for the class of problems considered in 
Chapters 2 and 3, the wobbling (nutational) motion separates from the 
spinning motion in the linearized equations. For any flexible 
appendage lying in a plane perpendicular to the spin axis, and passing 
through the system center of mass, the linearized equations of motion 
separate into two groups. The wobbling motion consists of the motion 
described by the nutation angles and the transverse vibrations, and 
the spinning motion consists of the rotation along the spin axis and 
what is often referred as the in-plane deformations. It is, in fact, 
this last point that justifies the consideration of only the transverse 
vibration in the previous study. For most of the aerospace applica- 
tions, the attitude stability of the spacecraft is of main interest 
and is only affected by the out-of-plane deformations. It is also 
observed that the attitude stability is affected by the anti-symmetric 
modes of symmetric spacecraft only, and through this, the assumption 
made earlier to consider the translational motion of the central rigid 
core of the system as being not affected by the deformations, is made 
acceptable. The anti-symmetric modes are, indeed, those which keep 
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the center of mass of the total system at rest with respect to the 
rigid core. 

The development of stability criteria for free spinning bodies 
has been the basis of numerous technical papers. Rigid body analysis 
prior to the flight of Explorer I predicts a stable free rotation in 
inertial space if the angular velocity vector is directed parallel to 
a principal axis of either maximum or minimum moment of inertia. The 
analysis following Explorer I led to the general conclusion that for 
a flexible spinning satellite to exhibit stable free motion its axis 
of spin must be restricted to that of the principal axis of maximum 
inertia; this proposition is sometimes referred to as "the greatest 
moment of inertia" rule. However, one would expect that the last 
criterion is not sufficient to assure stability and it is not sur- 
prising at all to find that spacecraft with very large flexible 
appendages are less stable than quasi-rigid ones. It is in the 
examination and development of stability criteria involving the modes 
of vibration and the natural frequencies of the structures that most 
of the recent papers are oriented. 

It would be too lengthy to cover the numerous publications dealing 

with this last task; Reference 13 provides a current bibliography on 

this subject. We should, however, mention here that the problem is 

generally approached through two different procedures. Some authors 

examine the stability problem by using a Routh-Hurwitz analysis, and 

some prefer to use a Liapunov analysis. In our study, we decided to 

employ the latter approach, proceeding in parallel with the work done 

13 

by F.J. Barbera and P.W. Likins. In their work, the authors develop 
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general stability criteria for a flexible appendage described by a 
collection of particles. They used Liapunov’s second method as the 
basic analytical tool and by specializing in the more restrictive case 
wherein the appendage lies in a plane containing the center of mass 
and orthogonal to the spin axis, they are able to come up with some 
analytical criteria. 

The stability analysis developed in Chapter 4 is based on a 

similar approach. Liapunov stability theorems are employed with the 

Hamiltonian of the system, constrained through the angular momentum 

integral, used as a testing function — a method proposed by 
14 

R. Pringle, in order to circumvent problems related to the negative 
definiteness of the Hamiltonian time derivative. In our formulation, 
the flexible appendage is described as a continuous elastic body, as 
opposed to a collection of elastically interconnected particles. We 
kept the derivation as general as possible. However, in order to 
come up with specific criteria, our study has been restricted to more 
particular cases. Finally, the analysis was brought down to a level 
amenable to literal stability criteria by truncating the number of 
normal coordinates to a single mode. The stability criteria emerging 
from this study represents, due to the special form adopted for the 
Hamiltonian, conditions for stability that are sufficient and (except 
for a few recognizable singular cases) necessary as well. 
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CHAPTER 2 

DYNAMICS OF ROTATING ELASTIC BEAMS AND CABLES 

2.1 Introduction 

In the following chapter, we face the problem of a dynamical 
study of two elementary structures characterized by a dominant effect 
of the forces induced by spin over the flexural rigidity. In the first 
part, a modal analysis is done for a rotating beam clamped normal to 
the spin axis. The study is oriented towards the application in 
rotation-stabilized space vehicles having a nondeformable frame to 
which flexible rods are attached. In the second part, the modal 
analysis of a taut cable clamped at its extremities is done. 

In both studies, we will consider vibrations only in the 
meridional-direction. We will also ignore deformations which are 
present only because of the Poisson effect. We will assume, at the 
outset, that the motion of the base to which the elastic structure is 
attached is uninfluenced by the elastic vibrations of the appendage. 

2.2 Dynamics of a Rotating Elastic Beam 

We now specialize in the study of a deformable element which 
consists of a long flexible beam undergoing transverse vibrations. 

When an elastic beam is normal to the spin axis of its inertially 
rotating base, it does sustain deformations in the steady state con- 
figuration in which it remains straight and aligned with a radial line 
emanating from its base. Therefore, as emphasized by P.W. Likins, 

F.J. Barbera and V. Baddeley one must consider nonlinear strain- 
displacement equations if deformation variables are to be measured 
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from the undeformed state. The requirement for the retention of 
second degree terms comes from the fact that steady state deformations 
induced by constant spin are not arbitrarily small, and cannot be 
included with the arbitrarily small deviations from the steady state 
deformations in the linearization process. This great difficulty of 
nonlinear elasticity explains generally why applications in the liter 
ature are restricted to beams. For those particular cases, the equa- 
tions of motion are typically derived by means of procedures which 
rely, from the outset, upon the availability of solutions for the 
steady-state load distribution and deformation of the elastic 
continuum. 

In our analysis, we neglect rotation of the transverse cross 
section of the thin circular rod or beam and consider the linear 
density and the flexural stiffness El to be constant. The rods are 
also considered as being much longer than the frame dimension. We 
therefore assume the point at which the rods are attached coincides 
with the axis of rotation of the frame, as in Figure 1. 



Figure 1. Rotating Uniform Elastic Beam. 



Under those assumptions, the transverse vibrations of the clas- 
sical (Euler-Bernoulli) beam, subject to an external axial load P(n) is 
given in general form by (see Reference 15) 

3 4 w 3 i 3w , , 3^w 

El — T - -r- [P (n) v- ] + y — r = 0 

3 n 4 3 V 3n 3t 2 

where El is the flexural stiffness of the rod 
y is the linear density 
P (r|) is the external axial load . 

The boundary conditions are given for a rotating rod under the assump- 
tions of one cantilevered end and one free end by 


w(0) = 0 = |^ (0) = (L) = (L) . 

' 3n 3'n 

For the rotating uniform radial beam, the influence of steady-state 

centripetal accelerations can be represented by an "effective force" 

or "centrifugal force" given by * 

L 

/ 21112 2 2 

y n n d n = 7 yfi (1/ - n ) 

<1 


so that the previous equation takes the form 


_ T 3 4 w 1 n 2 

EI 74 - 2 ^ 0 

3n 


n 2 2 \ 9 3 w 

(l - n ) — - 2 n — 
3U 3 | . 


a 2 

a. 9 W _ „ 
+ y — r °* 

3t 


The artifice of the "effective axial load" permits this derivation of 
the vibration equation to be accomplished without the reliance on non- 
linear strain-displacement equations noted previously to be required 
for general derivations of small-vibration equations for rotating 
elastic continua. 
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Two remarks should be made about this last equation. First of 
all, this relation expresses the transverse motion of the beam free 
of any external forces or deviations of base motion from simple spin, 
coming, for instance, from the coupling between the transverse vibra- 
tions and any nutational motion of the rigid core. Also, it should be 
emphasized that in the previous relation, the deformations along the 
r| axis represent the result of two separate deformations: the steady- 

state extensional deformation from the undeformed state and the devia- 
tion from that steady-state. This last point could be overlooked in 
the previous relation due to the artifice of the "effective axial 
force," but is an unavoidable fact of the general continuous model. In 
the derivations done in Chapter 3, this last remark will become more 
apparent in the sense that the steady state deformation has to be 
computed first and the equations of the motion are then represented by 
the deviation from this steady-state. 

We use, now, the method of separation of variables: 


w(n,t) = <j> a (n) P a (t) 


Substituting into the previous equation, we have: 


P + oo P =0 

a a a 

- - 1 0 2 

dn 4 2 


2 _ 2 * d \ 


(i* - n ) 


- 2 n 


d<J> 


a 


dn 


dr) 


- oo <p =0 
a Y a 


where (*) stands for derivative with respect to time 


(0 represents the eigenfrequency of the vibration 


and the boundary conditions become 
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d<f> 


*a<°> ‘ 0 ’ dT <°> ‘ 


d 2 <j) 


a 


dn 


(L) = 


d 3 (J) 


a 


dn' 


(L) . 


With the introduction of the dimensionless variable, y = ^ , 
2 

and division by ft , the previous relation becomes: 


IV 1 9 » 2 

e <T V - i [(l-y )<j>'] - X <p = 0 

Y g 2 3 Y a a r a 


( 2 . 1 ) 


where 


(') stands for derivative with respect to the dimensionless 
variable y 
El 


£ = 


yft 2 L 4 

a) 2 

a 

ft 2 


with the boundary conditions 


V 0 ) - 0 - V 0) = 


V'(D 


In the parameter e, the high spin rate and the low flexural 

rigidity combine their effects to produce a small quantity. A typical 

-3 -4 

value for the parameter £ is given by the range of values 10 ...10 

(see Reference 2). As a result, the previous equation is suitable for 
the use of singular perturbation theory by the presence of a small 
parameter in the highest derivative term. 

Let us consider for the central portion of the beam, a solution 
for which analytical dependence on the small parameter £ is given by 


4> a (y,e) = h Q (y) + v^e) h^y) + v £ (e) h 2 (y) + . 


where V^(£) is a small o(e), and the V^(£) constitutes an asymptotic 
sequence of functions. 
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By definition, <}) (p) = o (’i'(lJ)) as y -► y if, given any 6 > 0, 

there exists a neighborhood of y Q such that |<j>| _< 6 | 'F | for y in the 

neighborhood; thus <j>(y) = oC^Cy)) if $/¥ -> 0 as y -*■ y . Also a 

o 

sequence 4> (y) , n=l,2,... is called an asymptotic sequence if 
( nH W * °'(*n (ll) ) 11 * V o • 

By definition, two functions <j),4' of y belong to the same equivalence 
class in a neighborhood of y Q if 

0 < lim wSf < °° * 

y +y 

If this double inequality is satisfied, we adopt the notation 
ord <f>(y) = ord H'(y) 

A partial ordering of equivalence classes is given by 
ord <|>(y) < ord 'i'(y) if lim ^ -»■ 0 

y - y D 

so that if 4>(y) = o ^(y)j then ord $(y) < ord ¥(y) . 

To shorten the writing, the subscript a, characteristic number 
of the mode will be omitted momentarily and we will use this subscript 
explicitly only when needed, in order to differentiate between differ- 
ent modes . 

2 

Similarly the eigenvalue A of the problem will be written 
2 

simply as A and it will be expanded in terms of the small parameter 
e in the following way: 

A 2 (e) = A 2 + k 1 (e) A 2 + <„(e)A 2 +... 
o l l z z 

where. constitutes an asymptotic sequence. 
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The corresponding differential equations defining h^Cy) and h^(y) 


are given by 


f [Ury 2 ) h( ^ +xJb£-0 


and 


\ [(l-y 2 )h']' + \l h L « 


2 K l (£) 

-A? h. if r = 1 
1 0 vAe) 


K,(e) 

if lim — ; = 0 

e - 0 V e) 


where the different possibilities were limited to the case where 
ord v^(e) > ord e, which will be verified later. We also neglect the 
possibility that v^(e) = o^K^(e)^ , because this would imply that 
A^ = 0. 

At this stage, one remark can already be made on the solution for 

h^. There are no apparent boundary conditions to the equation defining 

11q. But this equation is known as the Legendre equation, the solution 

to which is well known in terms of the Legendre polynomials (see 

Reference 17) . Even though it might seem that there would be a com- 

2 

plete freedom in the values of A^, we can see from here, already, that 
2 

the values of Aq are limited to a discrete spectrum based on the obser- 
vation that the point y=l constitutes a singularity for the Legendre 

" 2 

polynomials, and only a discrete set of values of Aq will give the 
value of the Legendre polynomials to be bounded at y=l. 

In order to make those observations more specific, we have to 
concentrate for a while on the boundary layers existing at both extrem- 
ities, whose solutions will give us the. constants of integration we 
need for the previous integration. 
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Boundary Layer Near y = 0 

In order to study the boundary layer near the origin, we have to 

introduce a stretched coordinate 

y = ~ assuming a(e) -*• 0 when e -»■ 0 

The corresponding asymptotic expansion valid near y = 0 is given by 


<J> a (y»e) = y Q (£)g 0 (y) + y 1 (e)g 1 (y) + y 2 (e)g 2 (y) +... 


where the y^(e) constitute an asymptotic sequence. Replacing the 
previous asymptotic expansion into the differential equation of the 
mode shape 4>(y) and considering that y = G(e) y and dy = a(e) dy, we 
obtain: 


y Q (e) d 4 g Q (y) y 1 (e) d 4 g 1 (y) 

a 4 (e) dy 4 a 4 (e) dy 4 


2 ~2 
i-q (e)y 
2 


a 2 (e) 


d 2 g Q (y) y-^e) d 2 g 1 (y) 

— - •+• — - - — 

dy a (e) dy 



+ a(e) y 


a(e) 


dg Q (y) y 1 (e) dg 1 (y) 
dy + a(e) dy + * * • 


” A 0 ^ y o^ e)g o^ + y i (e ^ 8 i^ + * * 

- k 1 (e)X 2 [y 0 (e)g 0 (y) +y 1 (e)g 1 (y) +. . . ] = 0 


A suitable boundary layer coordinate is chosen by the requirements that 
the higher derivatives are of the same order of magnitude near the 
origin or by considering the previous equation, in writing 
ey Q (e) y Q (e) 
a 4 (e) c 2 (e) 

2 

or a (e) = e and a(e) = 
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The boundary layer coordinate is thus given by: 



The dominant boundary- layer equation is given by: 


d 4 g Q (y) _ 1 d 2 g Q (y) 


dy 4 


2 j~2 

dy 


= 0 


Both boundary conditions at y = 0 have to be satisfied by g^: 

d g Q (y) 


d y ■■ g 0 = 0 at y = 0 

Defining d 8 0^ 

2 — = g o ^ y ' * we have 

dy 


d V (y) 1 * ~ 

- T Sn <y) = 0 • 


,~2 

dy 


2 6 0 


y //2 


Using the fact that the exponential growth (e^ ) cannot match as 

y 00 and taking into account the boundary condition, we obtain 
successively 

. ~~/ T y 

—■ c » e 
dy 

and 


g 0 (y) = 2 C o 


72 


-1 + e 


/2 


which represents a solution with one arbitrary constant but doesn't 
depend on the eigenvalue. 
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. We are now in a position where we can try to match this boundary 
layer expansion and the solution we found previously for the central 
portion of the beam. 

An intermediate limit suitable for matching near y = 0 is given 
by y fixed: 

n 


y 8 " — - fa— 

n(e) 


ri(£) -*■ 0 and 


n(e) 

JT 


+ °° as £ -> 0 


so that 


y = n y n -► o 


and y = 


1 — y ->-oo 

/r n 


Matching near y = 0 takes the form: 


lim 
£ -> 0 
y fixed 

n 


h Q (ny n ) + ^(e) h 1 (ny n ) +... 


- p o (e)8 o(“|" y n) - y l (£)g l(^" y r) + -**| = 0< 


Expanding h^ and h^ as Taylor series near the origin and using the 
solution previously found for gg(y), we have: 


lim | h Q (o) + ny n h^(o)+. . .+v 1 (£)h 1 (o) + v 1 (£)h|(o)ny n +. 


y fixed 

n 


- 2y 0 (e)C o 


JL 


T, , . ^ 

1 y -1 + e 

Jl Ji n 


+ . . . j 


= 0 


From here, we can see that the matching is possible only if the follow- 
ing equalities are satisfied: 



Another point can be noticed from the previous equation. The 
term 2 V1 q(£) C 0 appearing in the expression of can not be matched 
except by a suitable value for h^(0) and this requires that 

v, (e) = y.(e) = /e and h. (0) = -2 C 
xu l o 

which satisfies our requirement that v^(e) -* 0 when £ -*■ 0, and implies 

ord V^(e) > ord e an inequality that we still had to prove. 

We are now in a position where we can go back to the solution 

valid in the central region of the beam. 

The equation defining hg 

\ t(l-y 2 )h^]’ + A 2 hQ = 0 

is a Legendre equation. 

Defining the Legendre functions of the first and second kind by: 
r ' 

^ 1 (2A^, y) = 1 -(^ (2A 2 )y 2 + (^) (2A 2 ) (2A 2 - 6)y 4 -... 

■ * 

^2 (2X 0’ y) = y <2A 0 " 2)y3 + -” 

2 

we see directly that only the second expression y) satisfies 

the requirement 11 q( 0) = 0. We thus have: 

V y > “ k o ^2 (2X 0’ y) • 

2 tli 

The set of discrete values Aq which represents the zero order 

of the asymptotic expansion of the eigenvalue are in fact those values 

2 

which truncate the Legendre polynomials. For all other value of A^, 
the serie becomes infinite and unbounded for y = 1, which could not make 
physical sense for our problem. This last point will be shown later. 
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As a result the set of discrete values feasible for A^ is given 


by 


or in general 


Aq = 1,6,15... 


, 2 n(n+l) 

Aq = 1 ^ for n odd. 


Now that an expression has been found for h^Cy), let us look at 
the expression of h^(y). 

Before any computation for h^ , we will need a property of the 
mode of the previous problem, i.e. their orthogonality. In order to 
prove the orthogonality of the modes, we have to consider the following 
equations and boundary conditions: 

e 2 = o 




4> a (o) = - <f>;;(i) = ^"(D = o 

<i> e (o) = <J>£(o) = 4>£(i) - ^"(l) - o 

Multiplying the first equation by and integrating by parts from 
0 to 1, we obtain: 


e 4£"* 6 


- e A" 

J rs Ct 


'^3 d y - j ta-y 2 )*;]^ 


+ i f dy - x l / Vb dy ‘ 0 • 


The boundary conditions and a new integration by parts gives us 

,11 1 1 


Y cr 3 


+ E L dy + 2 / b-y 2 )*^ - x lf Vb dy ' °- 
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A similar expression can be obtained from the second equation and for 
2 2 

^ ^3 We ^ ave ^ tom the difference of the two expressions 
1 

f 0 Vs dy = 0 

The previous result has to be completed by the relation defining the __ 
norm of the modes and for orthonormal modes we also have 
1 

f 0 *■! dy = 1 

Expanding the modes through the asymptotic expansion we defined before, 
and considering that the previous expression has to be true for every 
small e, the last expression becomes, for instance: 

1 

hg(y) dy - 1 

1 

h Q (y)h 1 (y) dy = 0 
and similar expression for h^. . . 

We have to be careful in the use of the previous relation because 
we have to include under the integral sign the contribution coming 
from the boundary- layer solution. For our purpose, we already showed 
that the boundary layer near the origin is of thickness /e and, as a 
result does not contribute in the second integral. We will prove later 
that the boundary layer near the free end does not contribute to the 
previous integral too, and the two integrals we just wrote will, from 
there, be useful. From 
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and 


r 2 

/ h_ dy = 1, we conclude 

•'o 0 

Jo k 0 ^2 ^0* y ^ dy 1 



Considering now the differential equation defining h (y) we have 

*• X 

\ [ (l-y 2 )h^] + X 2 h x - 


o k, (e) 

- X 1»0 1£ -^--1 


V? 

< 1 (e) 

if »- 0 as E -*■ 0 

The second possibility is now to be cancelled if we want to have any 
chance to satisfy the integral 
1 

f h Q (y) ^(y) dy = 0 . 


As a result, the function < 1 (e) is. well defined and we have 


and 


k^Ce) = /e 

2 Ki-y 2 )^]' + -x 2 h 0 


(2.2) 


We also found before in the matching process near the cantilevered 
end: 

r- n.ty; 

M _ /o /i _ i 

or 


hAO) 

h’jCO)- = /2 Co = - -X- 


/2 


k 0 ^2 ( 2X 0 > °) = ^ Co = " h x ( 0 )/ »^2 


h i(0) = -k Q £' z <2XJ, 0)^ 
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We will now determine the value of A^ through the use of the 


differential equation (2) and its boundary conditions and also by 
using the orthonormality condition on h^. 

Multiplying Equation (2) by h^ and integrating by parts from 
0 to 1, we have 


\ [(l-y 2 )h|]h 0 


x x 1 

" 2 X [(lry 2 )h»]hj dy + A 2 f h^dy = ~ X l f Q h l dy 


Assuming h^(l) to be bounded, the use of the boundary condition of h^ 


and its normality gives us 


2 ) Q (1 -y 2) h l h 0 dy - X 0 f 0 h 0 h l dy ■ A 1 • 


Integrating by parts once again, we have 

1 1 t 1 

~2 ^ 1_y )h 0 h l Q “ 2 X t(1 ~ y h]L ^ ~ *° X^ 1 ^ = 


The differential equation defining h^ 
- ± Kl-y 2 ) h']' - + X 2 0 h Q 
brings us the final expression: 


or 


or 


- j h><0) h l( 0) - xl 


A 2 = — — h 2 (0) 
1 n rz 1 


a: = 


2/2 


1 

S2 


[^(2Aq, Oj] 


/ 


(2A 2 ,y)dy 
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Considering the definition of <5? ^(2Xq, y ) , we can see that 
(2X 2 ,0) =1 and 


X? = 


& /" ^2 ( 2 ^0’ y ^ dy 


For the first mode, corresponding to Xq = 1, the previous relation 
becomes 


X 2 = 


I y^dy Jl y 

•'0 


1 Jl 

0 


2.12 


2 

For the second mode or for X^ = 6, we have: 




& ][( y-|y 3 ) 2 dy 


63 

4v^2 


11.18 


0 


We should emphasize at this stage that the previous expression 
2 

found for X^ is simply based on the properties of the differential 

equation and did not use any other relation. It is easy to prove that 

2 2 

the previous values found for Xq and X^ are also those which satisfy 

the orthogonality relationship between different modes. For the zero- 

order terms of the expansion, the relation is simply given by the 

orthogonality property of the Legendre polynomials in the interval from 

o to 1. For the first order terms, the derivation is lengthy, but the 
2 

values of X^ are coming out from those relations. Now that an expres- 

2 

s ion has been found for X^, we can go back to Equation (2.2) 

\ [(l-y 2 )^]' + X 2 h 1 = -X 2 h Q ' . (2.2) 
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The general solution of the previous equation is given by a particular 
solution h^ which satisfies the boundary condition 

h L (0) = -k Q £?• (2 Aq ,0) /2 

and the general solution of the homogeneous part of (2) multiplied by 
a constant k^. This last solution is nothing else than h^, the solu- 
tion of the zero-order equation. So, in general, 
h ][ (y) «.h 1 (y).'+ k^Cy) . 

Out of the family of solutions, we have to choose the solution which 
satisfies the orthogonality relationship between h^ and h^ for one 
mode or 

1 1 1 
0 ■ f 0 h 0 h i dy ■ f 0 h 0 K i dy + k i f 0 h l dy 

or, finally 1 ■ 

k i ' - f 0 \ K i dy 

The knowledge of k^ specifies entirely the function h^ and as we will 
see later, this last step will be needed when we will match the 
central part expansion to the boundary-layer expansion near the free 
end . 

Boundary Layer Near the Free End ' 

We, now, have to look at the boundary layer needed at the free end 
in order to satisfy the boundary condition and confirm what has been 
said before about the value of h^ at this extremity. In order to pro- 
ceed, we have to introduce another suitable boundary layer coordinate 

"k 1“V 

near y=l. Let, then y = where <f>(e) "*■ 0 when e .-*■ 0. We 
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£ £ 

thus have <j>y = 1-y and dy = - (j)dy . The corresponding asymptotic 
expansion valid near y = 1 is taken as 

<f> a ( y » £ ) = <S Q (e) f Q (y*) + 6 1 (e)f 1 (y*)+. . . 

where the 6^(e) constitute an asymptotic sequence. 

Replacing the previous asymptotic expansion into the differential 
equation of the mode, we obtain for the basic equation 


S Q (e) d 4 f Q (y*) 6 1 (e) d 4 f 1 (y*) 


<}> 4 (e) dy* 4 



, *4 
dy 

* 2 

i— (1— 4»y ) 

6 0 (e) 

2 * 
d f o (y } 

6,(6) 

2 

4> 2 (e) 

dy* 2 

<p 2 (e) 


+ • 


i 2 f L (y*) 


*2 


- (1— 4>y ) 


dy 


v e > d f o (y *> ^(e) d 

+ ^ + - 


-■X^[6 0 (e)f 0 (y*) + 6 1 (e)f 1 (y*)+...] 


-K^e)^ [6 Q (e)f 0 (y ) + 6 1 (e)f 1 (y ) + . . . ] = 0 


(2.3) 


* 2 * 2*2 

Rewriting l-(l-<|)y ) as 2<py - <p y in the previous equation, we 
choose a suitable boundary layer coordinate by the requirements that 
the higher derivative terms are of the same order of magnitude near the 
free end. Thus by considering the previous equation, we have 


eS 0 (e) 6 Q (e) 


<!> 4 (£) 


<Ke) 


or 


4> 3 (e) 


<P (e) = e 


e and 
1/3 


so the boundary layer coordinate is given by 


* - ! / 1/3 

y = 1-y/e 
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The dominant boundary layer equation is given by 


d 4 f 0 (y*) * d 2 f Q ( y *) df o (y * ) 

y 


dy 


*4 


d y 


*2 


= 0 


dy 


Both boundary conditions at y = 1 have to be satisfied, by f^(y ) at 
* 

y = o 


2 , * 3 * 

* %<y ) n 0 (i ) , 

= 0 at y = 0 


d y 


*2 


dy 


*3 


The previous equation can be written 


d 4 f o ( y*) 


dy 


*4 


4 . df n 
* 0 
y — 

dy L dy J 


A first integration gives us as a solution 


d 3 f 0 (y ) 


dy 


*3 


* df o 
y — 
dy 


where we already used one of the boundary conditions in removing the 
constant of integration. Defining 

df o * * 

" f 0 (y >’ 


we have 


dy 


d 2 f 


0 * * 
= y f n 


*2 
dy 

^ ^ jc 

with the boundary condition df^/dy = 0 for y =0. 

This last equation is known as the Airy equation, the solution 
of which is given in terms of the Airy integrals. We, thus, have 

f* = A Q Ai(y*) + B Q Bi(y*) 
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The asymptotic expansions of the Airy integrals can be found in the 
literature (see Reference 18) and are given by 

Ai(z) ~ — n " 1/2 z " 1/4 F ( — — • - — ) 

2 2 0\ 6 ’ 6 ’ 25 ' 


and 


w x tt-1/2 -1/4 5 , 1 5 1 \ 

Bi(z) ~ n z e 2 F 0 ( 


where 


5 = § z 3 ' 2 


and 2 ^q stands for the generalized Hypergeometric series. 

A direct observation of the expression for Bi(z) shows that it 
is inappropriate for matching due to its exponential growth, and only 

& * 

f Q = A Q Ai (y ) 

has the exponential decay property. 

We also have to satisfy the boundary condition 
* 


df 


0 * 

— y = 0 for y = 0 

dy 

and the power series expansion of Ai*(z) near z = 0 is given by the 
following expression: 

Ai’(z) - -b () F 1 ( ; | ; 5 2 /4) + \ a z 2 ^(j j ; ? 2 /4) 


where 5 has the same meaning as previously and the constants a and b 
take respectively the values 0.355 and 0.258. 

As a result, the satisfaction of the last boundary condition 
requires that A^ = 0 and we are left with 
* „ * 

f 0 (y ) = o 

or f o = V 
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If we spent some time in the development of the solution for f^, 
it is due to the fact that the previous development will be found 
helpful in the search for the solution of the other functions in the 
boundary-layer. 

To investigate a little further the boundary layer solution, we 
have to try to match the solution valid in the central region of the 
beam to the expansion valid near the free end. In order to accomplish 
this last objective, we have to expand the solution found in the 
central region near the free end. We have: 

4> a ( y »e) = h Q ( y ) + ^ h 1 ^y)+... 

Defining the variable y' = 1-y, we have: 

<f> a (y',e) = h 0 (l-y') + »/e h^l-y'H. . . 
or 

y . 2 

<J) a (y\e) = h Q (l) - h^(l)y' + h”(l) “+... 

± /z h^(l) - y/e h^(l)y'+... 

For matching, an intermediate limit suitable for matching near y = 1 
is given by y^ fixed: 


y n n(e) 

n(e) -*■ 0 and 


where y' 


1/3 


e 

* n y 


so that y' = r) y^ -*• 0 and y ’= 


is always positive. 


->- OO 


-> OO 
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Matching near y'=0 takes the form: 

r 2 2 

r) y 

lim h Q (l) " h o^ 1 ^ n y n + h (j^ “ Jj +•••+ ^ h 1 (l) 

y n flxed '- .6 h*(i)n y n +... 

- 6 0 (e)D 0 - 6 x (e) f 1 (^ 1 /3^ “ 6 2 (e)f 2 ( 


From here, we can see that the matching suggests the use of the follow- 
ing equalities: 


V e) 

o 

= e 

and h Q (l) = D Q 

6 1 (e) 

- £ 1/3 

. , v 2/3 

63 (e) = e 

6 2 (e) 

- £ 1/2 

* / \ 5/6 

64 (e) = £ 


As a result, we can see that the determination of seems to be free 
and the conjecture made earlier about the value of hg(l) is perfectly 
valid. 

The boundedness of the value of h^(l) is the real boundary con- 
dition that should be used at the free end, and the value of is 
thus known! 

Now that the matching process has suggested an order of magnitude 
1/3 

of £ for <S^(e) in order to match the linear term in r) y^ appearing 

* 

in the previous relation, let us compute the expression for f^(y ) . 

* 

The differential equation that f^(y ^ roust satisfy is obtained 
by picking the term of order of magnitude 

E 6 1 (£) 

4 £ 

<f> (e) 

in the general expression (3) 
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d 4 f 1 (y*) 


2 * 2 * 
* d 2 f(y ) „*2' d Z f,(y ) 


dy 


*4 


- y 


dy 


*2 + 2 


dy 


*2 


+ y 


* dVy*) 


dy 


dy 

with the boundary conditions 

2 

d f 1 <y ) 


d f i(y ) 2 * 

* X 0 f 0 (y } = 0 


*2 

dy 

3 

d f i<y ) 


dy 


*3 


= 0 


= 0 


for y =0 


Using the constant for fg(y )» we are left with 


- y 


df 


*4 


*2 


I = A 2 D 

* A o °0 


dy dy dy 

In this last linear differential equation, the homogeneous part 
is exactly the same as the one defining f^. 

A particular solution is given by taking a linear relation for 


f^Cy ) such as 


f T^y > “ -X n D ny 


The boundary conditions on f^(y ) cancel once more the two Airy 
Integrals, solutions of the homogeneous part and as a result 

f i.(y* } = - x o K o (1)y * • 

Going back to the expression computed for the matching near y’=0 
we have to show that 


-h^(l) + AqIIqU) = 0. 
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The differential equation defining is given by 


\ (l-y 2 )h|j - y + X 2 h Q = 0, 


evaluated at y=l, we have: 

-hj(i) + x 2 h 0 (l) =0. 

In the matching, valid near the free end, we observe that the 
limit process suggests to introduce a term of order of magnitude 
6 2 (e) = in order to take care of the presence of h^(l). We will 
thus introduce here into the asymptotic expansion valid near the free 
end the expression 

f 2 (y*) = D 1 

where is equal to h^(l). There is, as expected, no new information 
added to the problem, the value of h^ being already completely 
specified. 

Proceeding to the next term, the matching process of the quad- 

2 2 2/3 

ratic term in n y^ suggests for ^(e) an or der of magnitude of e 

* 

The equation defining ) is given by considering the terms of order 

of magnitude 


e 6 3 (e) 
4> 4 <e) 


“ 2/3 = e 1/3 


.4/3 


in the expansion near y = 1 or 


d f 3<y ) * d ‘ £3 


*2 d 2 f 


dy 


*4 


y *2 + 2 

dy dy. 


1 

*2 


d ^3 * 2 * 

— + y — - A 0 f x (y > 

dy dy 


= 0 
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with the boundary. conditions 

* 

for y =0 

We already found that f^(y ) ~ ~^q h^ClJy . So the previous expression 
becomes : 

d f 3 (y ) * d f^ df^ 4 * 2 * 

*4 ~ y *2 * " X 0 + X 0 h 0^ y 

dy dy dy 

The solution of this equation is given by a particular solution and the 
general solution of the homogeneous part. The particular solution is 
taken as 

f 3 p(y*) = °3 y * 2 

where 9 

D 3 ’ -f - < X 0 - « h 0 (1) ' 

The general solution of the homogeneous part has to satisfy the 
requirement on the third derivative due to the choice of the particular 

solution and as a result, the general solution is given by 

* 

y 

£ 3h (y * )=A l jfAi(5)dC + D 4 

as can be seen from previous development , where the other Airy Integral 
has been cancelled for its exponential growth. 
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The final expression for f^ is then 


* y r a 2 

f 3 (y*) - A 1 J 1 Ai(£)d£ + -jr- (A 2 -l) h Q (l)y* 2 + D 4 


The determination of is given by the use of the second boundary 
condition 


or 


or 


d 2 f. 


dy 


*2 


= 0 


for y = 0 


A 2 

* A 0 


1 Ai'(y ) + ~ 2 ~ ( X 0 - D ^q(I) = 0 for y = 0 


A i = - 


Xg(> 2 - 1) h 0 (l) 

2 Ai ' (o) 


where the value of Ai'(o) can be found in tables. (Reference 18). 
Going back to the matching procedure, we note that the term 

v“’ 


/ 


Ai(OdC + D. 


presents the exponential decay which allows us to neglect its contri- 
bution in the process, and we are left with 
2 2 2 

^ y n A n 2 2 2 

h o a> -it - 4 <x o - « h o (1)n *n ■ °- 


This last relation represents an identity as we now show. 

Considering the expression defining h^ and expressing this rela 
tion in terms of the variable y', we obtain: 


2x!.-.x 


,2 


K - (l-y')h' + A 2 h Q = 0. 
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Expanding fig, tig and hg in a Taylor series expansion valid near 


y'=0, we have 


h 0 (1 > 2 

h o( 1-y ') = bgd) " h^(Dy' + ~Yj — y'+.. 

h o' (1) 2 

h*(l-y*) = h'(l) - h” ( 1) y * + y* 2 + 


hg (1-y 1 ) = h^(l) - hg' (l)y '+. . . 


Introducing the last expansions into our previous equations, we obtain 
the following identities: 


. -hj(l) + Ajhg(l) •- 0 

2h^(l) + hj(l) - Ag h* (1) = 0 . 

This last expression can also be written: 

2hg(l) = (Ag - 1) A l hg(l) , 

which represents precisely the equality coming from the matching 
process . 

Pursuing the matching even further, the next term to consider is 
f^(y ) corresponding to <5^(e) = £ . The equation that f^(y ) must 

satisfy is represented by the terms of order of magnitude: 


e6 4 (e) 


dll - e V2 

£ «/3 - 


in the development (3) , or 


d f 4 (y ) * d f 4 (y ) * 2 d f 2 (y ) df 4 (y ) 


*\ 

* ^2^ '2*2* 

+ y -Sr - - f 2 (y > - A 1 f 0 (y > = 0< 
dy 
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Taking into account the fact that 

f 2 (y*) = h 1 (l) and 

f Q (y*) = h o (1) » 


we have 


d\(y > 


dy 


*4 


- y 


* d V df 4 (y*) 


dy 


*2 


dy 


Aq h 1 (l) + A^ h 0 (l) 


The use of the boundary condition for 4 limits the solution for 
f 4 (y*) to 

f 4 (y*) = h i (1) y* " x i h o (1) y * • 

In the matching process, we are left with 
-h’(l) + Aq hid) + A* h Q (l) 

which is identically equal to zero if we consider the equation 
defining h^ evaluated at the point y=l or 

-4(1) + Aj 4(1) = -A* h Q (l) 

To complete the study of the boundary-layer solution near the 

free end, let us discuss briefly the term f,. of the development of the 

boundary- layer solution. For this purpose, we have that < 5 ^(e) = e, and 

* 

the equation defining f -(y ) is given by 


d 4 f 5 (y*) * d^f 5 (y") ' y *2 d ^y") df 5 (y") 


2 * 
1^. - / N 


. 2 . , * 


dy 


*4 


- y 


dy 


*2 


2 , *2 
dy 


dy 


* df 3 <y > 2 * 

+ y — ^r- - x o f 3 (y > = 0 
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with- the Boundary Conditions: 


d 2 f 5 (y*) 


*2 


= 0 


dy 


3 * 

d f 5 (y ) 


dy 


*3 


= 0 


for y = 0 


Only the asymptotic behavior as y ■> 00 is essential for matching and 
this is easily found. 


Writing for the expression already found for f^(y ) the following 


expression . 


a: 


f 3 ( y } (A 2 - 1) h 0 (l) y* 2 + TST 


where TST stands for transcendentally small terms, always neglected 
for matching, the equation defining f is given by 

d f 5 < y ) * d f 5 (y J df 5 (y ) A 0 '2 ' m *2 

- y = — (A 0 -l)h 0 (l)y 


dy 


*4 


dy 


dy 


3A 


0 2 *2 

(A q -l)h 0 (l)y + TST. 


The equation for f,_(y ) tan be integrated once and gives us: 


3 * 

d f 5 (y ) * df 5 

*3- - y — 

dy dy 


A 0 2 ^ A 0 2 

(A 2 -!) h 0 (l) ~ ^ (A 2 -l)h 0 (l) 


*3 

I 

3 


+-p 5 + TST. 

The calculation for f^ is carried out simply from the differential 
euqation, noting that 


d 3f 5 (y*) . df 5 

« y — * as y -*■ °° 


dy 


*3 


dy 
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The rationale beyond this last argument is exactly the reason why the 


matching works , so that we have 

df 5 _ < X 0 ‘ 1) 4 2 

* 4 <*0 ~ 3 *0^ h 0^ 3 

dy 


*2 D 


-5 + I_ 
* * 
y y 


d 3 f 5 (y*) 


dy 


*3 


+ TST . 


It follows that 


d 2 f, (X 2 -1) . , * D_ . ' d 3 f c 

5 ? 0 n 4 ,,2*^ n s y_ , 5 _ _ _1_ 5 

*2 2 k 0 0' 0 V ^ 3 *2 *2 *3 

dy /. y y dy 

1 d f 5 

+ — — + tst, 
y dy 


and 


d 3 f t 


dy 


*3 


<\> ~X) 42 1 

(Xj -3X 2 )h 0 (l) + CX-ijH... 

y 


Finally, we have: 


df, 


dy 


: < x o -*> 

12 


4 2 *2 

(XJ -3X‘) h 0 (l)y - 


4 2 

V C^-3xJ)h 0 (i) 


+ °( -Jj) + 1ST , 


£ 5 * 


(^0 "X) 4 o *3 

36 (X 0 -3X 0 )h 0 (l)y 


<^0 X ) 4 2 

D 5 + 6 ~ 3 \P h 0^ 


Jin v* 


+ d 6 + °( -Jj) + TST . 

y 


The constants of integration and are then found from the 
matching. The other constants of integration for each solution are in 
the transcendentally small terms and are found from the boundary con- 
ditions at the origin. 
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In the matching conditions, expanding the various terms and 

neglecting transcendentally small terms, we have: 

(»0-« 4 2 

D 5 = - -S-g- <X 0 - 3X 0 ) h 0 (l) 


D 6 = h 2 (l) . 


We also have to show the identity: 


h"'(D x 2 4 2 

3! + 36 (X 0 3X 0 ^ 


= 0 


All the other terms omitted vanish more rapidly than those matched. 

By going back to the development we already used for the identity found 
for f^» we also have: 


3 hjj'(i) = (X; -3) hgd) 


or 


3 hJJ'd) - 


(^0 -1) 4 2 


2 (A q -3A q ) h Q (l) 


The last identity is thus proven. 

Having now completed the modal analysis of the uniform classical 

beam, we summarize the results by recording the final expression for 

the eigenf requencies , or 

a ) 2 = A 2 ft 2 = fi 2 (A 2 + v'e A 2 + . . .) 
a a 0 1 

,2 n(n+l) 

A q = — for n odd, 


where 

and 


X l 2 = 


± • 

/2 (2A 2 ,y)dy 


This section concludes with a comparison of previous work and 

19 

that of the author using the results presented by F.R. Vigneron. 


In 


Reference 19, the author presents the relation for the lowest loaded 
natural frequency of a radial boom as 

W 1 = V + 1:193 fi2 ’ 

2 

where u)^ R stands for the unloaded lowest frequency for a uniform beam. 

2 2 2 4 

The quantity is recognized as = (3.515) EI/yL , so the expres 

sion of co 2 becomes 

w 2 = (3.515) 2 —p + 1.193 ft 2 . 

yL 


The definition of e gives us 


El n 2 
— t- = eft , 


so this expression can be written 


u) 2 - ft 2 (1.193 + (3. 515) 2 e) . 

2 

Comparing the above expression with our development of or 


w 2 = ft 2 (l + 2.12 Jz ) , 
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In the interpretation of these results, we have to keep in mind 

that the result presented by Vigneron represents a better approximation 

for increasing values of £ while our approximation becomes more precise 

as £ decreases. The limiting behavior of a beam when £ goes to zero is 

2 

given by the cable and the true value of the ratio (oj^/Q) is given, 
for that case* by 1 (as opposed to 1.193). 

After having considered the case of the classical rotating beam, 
we will now consider the case of a rotating cable of small flexural 
rigidity when clamped at both ends. 

2 . 3 Dynamics of a Taut, Rotating. Elastic Cable 

In this derivation, the same assumptions as those made earlier, 
i.e. to consider the motion of the spinning rigid ring to which the 
cable is attached as not being affected by the transverse vibration of 
the cable, and also to ignore the rotation of the transverse cross- 
section, will be made. The total system is also considered as being 
in rotation with an angular velocity around the axis of symmetry of 
the ring , to which the cable is clamped (see Figure 3). 

For the special type of configuration under investigation, it 
must be emphasized that the spin has essentially a destabilizing 
effect. It is known that the influence of spin on flexible spacecraft 
is characterized by the presence of three effects: preload, Coriolis 

coupling and centripetal acceleration. The effect of preload is of 
main concern for this structure, because of its dependence upon the 
configuration (orientation of the flexible appendage with respect to 
the rigid core). For this particular structure, the preload modifies 
seriously the stiffness properties of the structure and justifies the 
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introduction of a built-in tension that we will assume to be big 
enough to ensure a tension everywhere within the cable. 



Figure 3. Rotating Uniform Taut Cable. 

In what follows we do intend to invest the cable with some small 
flexural stiffness, so that the governing equations of motion are the 

• V ' ‘ ' ' ' 

fourth order partial differential equations of the beam, rather than 
the classical second order partial differential equations of the taut 
string. Under the previous assumptions, the transverse vibrations of 
the classical (Euler-Bernoulli) beam subjected to an external axial 
load P(£) is given in general form by 



where El is the flexural stiffness of the beam or cable 
U is the linear density 
P(£) is the external axial load. 

We should notice here that the equilibrium position of the taut 
cable is given when the total system is rotating at an angular velocity 
$2, with the cable straight along the £ axis. No displacement can occur 
either in the r| direction or the C direction. The position vector £ 
appearing in the last relation is in fact the position vector at 
equilibrium. If we stress this last point, it is due to the fact that 
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the equilibrium state is not a zero-stress state but there exists a 
steady state stretching into the cable, steady state stretching which 
is already included into our position vector £. We assume that this 
last fact doesn't alterate the constancy of the linear density or the 
flexural rigidity. 

The boundary conditions are given for this problem by: 
w(±L) =0 = J| (±L) 

where L stands for half of the total length of the cable. 

For the rotating uniform taut cable, the steady state axial force 
is given by the constant built-in tension T applied to the cable and 
by a compression force R resulting from the integral of the "centri- 
fugal forces" applied to the elementary masses. If we consider a 
cross-section of the cable, located at a position vector £, the system 
of stress applied to the whole section is given by 

P(5) = T - J yfi 2 5'dC 

0 

or 

2 £ 2 

no = T - 2 • 

One remark should be mentioned here. One of our assumptions, 
when we introduce the built-in tension T, was to make sure that every 
elementary element along the cable was in a state of tension. With 
the definition of the axial load, it is the same constraint as assuming 

P(Q > o 


that 


for every £ or 



Introducing the value just found for P(£) into the previous 
relation, we have: 


etIji.J 

9S 4 " 




_ L 


The introduction of the dimensionless variable s = ~ transforms the 

La 


previous equation into 


El 3 4 w ufi 2 9 L 2 2. 9w 1 

* L (k - S) ^J 


+ y 


9 2 w 

9t 2 


= 0, 


2 2 

where we introduce the definition T = k T . The constant k is a 

o 

known quantity which has to be greater than one. 

The last relation can also be written: 


9 w 


9s 


4 9s 


L _ §1 ) 

V k 2 / 


+ _2_ i* = o 
2 2 2 * 
ftV 9t* 


where e stands for ^El/yflW . 

Once again, we observe that the low flexural stiffness is com- 
bined with the high spin rate to produce the small parameter £. 

Written in a slightly different way, the parameter £ is also 
El 


£ = 


l 2 t 


and is a measure of the ratio of the flexural stiffness to the built-in 
tension. 

, The boundary conditions for the problem are also 
w(±l) = (±1) = 0, 

conditions which express the physical clamping of the cable at its 
extremities . 
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Now that the dynamic equation of the rotating cable has been 
obtained through the use of the well-known relation of a classical 
beam subjected to an external load, we could in a completely similar 
way have derived the previous relation by solving first the steady 
stretching of the cable, and then determine the transverse vibration as 
being a perturbation with respect to that steady-state of deformations. 
The last derivation would give in general a partial differential 
equation with nonlinear coefficients, and the first approximation of 
the coefficients would give exactly the same relation as the one we 
derived before. This last feature was mentioned before as being a 
recognizable fact of continuous modeling, and already included into the 
general derivation of the transverse vibrations where the substitution 

mi _ e 3 u 

A 

has been done into the derivation, where 

u expresses the deformation along the cable axis 
A expresses the cross-section area of the cable 
E expresses the modulus of elasticity. 

The procedure we just outlined was the approach used by Abel and 
Kerr,^ when they derived the dynamic equation for the transverse 
vibration. 

We now solve the transverse vibration equation, by using the 
separation of variables 

w(s,t) = <J> a ( s) P a (t) . 

Substituting into the previous equation, we have 
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a 


a a 



where (•) stands for time derivative, w represents the eigenfrequency 
of the vibration. The boundary conditions become at the same time: 

d<j> 


As announced in Chapter 1, the method of matched asymptotic 
expansion will be used in this problem, but a suitable change of vari- 
ables is recommended first, so we define 
s = kx. 


We then have 



7T tU-* 2 >VJ 



2 2 
Ok 



where (') stands for spatial derivative with respect to x, of 


c tv 2 * 2 

— <j> - [(l-x^)4> •] - <p = 0 

, 2 1 r a a a 

k 


(2.4) 


2 2 2 
where A replaces 2 U) /f2 
a v a' 


The boundary conditions become: <j> a (±k =.tj) ’ (±k "*") = 0. The para- 

meter k is always bigger than one and consequently k ^ is < 1. 

Let us consider for the central part of the beam, a solution for 
which analytic dependence on the small parameter E is given by: 

4> a (x,e) = h Q (x) + v 1 (e)h 1 (x) + v 2 (e)h 2 (x)+. . . 
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Once again, the subscript a, characteristic number of the mode, will 

2 

be omitted when not needed. Similarly the eigenvalue A of the problem 
can also be expanded in terms of the small parameter e in the following 
way : 

A 2 (e) - X 2 + K 1 (e)X^ + k 2 (e)A2 +. . . 

where k (e) constitute an asymptotic sequence. 

The differential equation defining h^Cx) is then given by 

[<l-x 2 )hj]' + X 2 h 0 = 0 (2.5) 

and is a Legendre equation. 

Defining the Legendre function of the first and second kind by: 

^(X 2 , x) = 1 - (^y) \ 2 q x 2 + (^i) A 2 (A 2 - 6 )x 4 +... 

^ 2 (A 2 , x) = x - (yy) (A 2 -2 )x 3 +..., 

we see directly that corresponds to the even modes and corres- 
ponds to the odd modes. 

We will now. discuss briefly some of the characteristics that we 

til 

are able to deduce from the zero'" order expansion of the modes based 

on the fact that we are expecting boundary conditions for the Legendre 

polynomials of 'the form hg(±k ^) = 0. This last point can be seen 
physically in . the sense that we are not expecting a jump for the modes 
near the extremities and will be shown later by the process of matching 
the asymptotic expansions. 

Based on the boundary conditions hg(±k = 0, we see that the 
eigenf requencies are the solution of the transcendental equation: 

^(A 2 , k' 1 ) = 0 and SB 2 (A 2 , k" 1 ) = 0 . 
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One remark should be made here. Because k can take on only positive 
values less than one, we observe that the singularity of the Legendre 
polynomials at x=l is avoided. 

2 2 2 

We first study the even modes. Defining the ratio A^/k = y^, 

2 2 

we will try to plot the variation of y^ with respect to k . The 
2 -1 

relation k ) = 0 becomes 


1 2,1 2,2 6 s 12,2 6 w 2 20 . 

1 " 2! y 0 + 4! y 0 (y 0 “ .2 6! y 0 (y 0 2 ^ (y 0 " 2 * 

k k k 


Some discrete values for the plotting are easily found like: 

2 

2 6 : y f) ' o o 

y Q = — gives us 1 j = 0 or y Q = 2 and k = 3 

k 


Pq = 14.8 and k 2 = 1.35 

y 2 = 2.32 and k 2 = 8.6 

2 

We remark easily here that the value k =8.6 corresponds to the third 

2 

mode and the value k =1.35 corresponds to the first mode. By first 
mode, we mean that the corresponding solution has no point of zero 
displacement except the extremities and by third mode, we mean that 
along the span of the cable, there is two points with a zero transverse 
displacement. 

Some interest lies also in the limiting case given when the 
2 

parameter k goes to infinity or when it goes down to one. In the 

2 2-1 

limiting case, when k -*■ the expression k ) = 0 becomes 

2 

the expansion of cosPq, and as a result the corresponding value of y^, 
valid for the first mode, is given by 
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2 

This last result is reasonable, because the case k -*■ 00 , corresponds 

to the ordinary vibrating string. The effect of the distributed forces 

induced by rotation becomes indeed negligible. 

2 

The other limiting case k + 1 is more difficult to handle due 

to the singularity at x = 1. In order to study this last limiting 

case, it is indicated to look for the intersection of the curve we are 

2 2 

looking for with the hyperbola defined by y^ = £g/k where is a 
small parameter. 

2 —1 2 e 0 

The intersection of ,5 ?^(Aq, k ) = 0 with ^ is given by 



Neglecting the higher order terms in e^, we have 



which can be written as 


1 + -f log(l - ~ ) = 0 or 1 - ~ 
k k 


-2/e, 


so 


k 2 -■ 


-2/e, 


-2/e, 


= 1 + e 


1-e 


As a result, we thus see that when y^ -*■ 0, we have that k differs 

from one by a transcendentally small quantity. 

We finally have for the first mode the following graph (Figure 4) 

2 

The previous graph shows also the asymptote when k -* 00 , which also 
represents the value of the eigenfrequency of the nonrotating taut 
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Figure 4. Eigen Frequencies Corresponding to First Mode. 

string. The hashed part of the graph represents then the decrease in 

the eigenfrequency due to the rotation, decrease in the eigenfrequency 

which is often referred to as the destabilizing effect of rotation. 

As expected by the previous remark, the destabilizing effect decreases 
2 

when k increases, which corresponds to the decrease of the effect of 

the forces induced by spin. We could also show how the first even 

2 

mode varies as a function of the parameter k (Figure 5). 



Figure 5. Mode Shapes Corresponding to First Mode. 


A similar development could be done for the odd modes. 

2-1 222 
The relation ^ ) = 0 becomes, after defining A^/k = Pq 
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We could once again derive some discrete value of y^ and k . The 

2 

limiting case, when k -*■ °°, reproduces here again, the result we might 

2 

expect from the taut vibrating string. In the case, where k -*■ 1, the 

2 

corresponding value of Pq is once again obtained by looking for the 

2 2 

intersection of the graph with the hyperbola y^ = 2+e^/k and as we 

2 

might expect, the resulting value of y^ is 2. All the results are 

summarized in the following graph (Figure 6) , which represents the 
2 2 

variation of y^ with respect to k for the second mode. 



The destabilizing effect introduced by the rotation is again 

apparent on the previous graph. By second mode, we mean, that there 

exists along the span of the cable one point which doesn't experience 

any transverse vibration. We could also show how the second mode 

2 

varies, as a function of the parameter k (Figure 7). After this 
discussion on the zero-order solution of the eigenvalue problem, we 
will now turn our attention towards the next order of magnitude of the 
solution. However, before any attempt could be done in this sense, 
we will have to look at the solution valid near the extremities of the 
cable where we still have to prove that hQ(±k = 0. 
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0 1 

Figure 7. Mode Shapes Corresponding to the Second Mode. 

Boundary Layer Near x = k ‘ 

We now plan to look at the boundary layer solution near x=k 
It is obvious that the boundary layer near x=-k ^ could be obtained in 
a completely similar way and will thus be omitted. 

In order to construct the boundary layer expansion, a suitable 
bound ary- layer coordinate has to be chosen, such that the higher order 
derivative terms dominate the equation. 

where 4>(e) 0 when e -* 0 

and dx = -<}>dx 

The corresponding asymptotic expansion near x = k ^ is chosen such that 
<}> a (x,e) - 6 Q (e) f 0^ x ^ +<s 1 C £ ) f • • 

where the form an asymptotic sequence. The basic equation 

becomes with the previous substitutions: 


Let 


x = 


k 1 -x 
(j)(e) 


so 


<j)x = k ^ - x 
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r« 0 w i 4 f 0 , 

S 1 (£) 




! <j> 4 d ; 4 

T 

4> 4 

d; 4 ' ' 



- [i - (k _1 - 

■<}>x) 2 ] 

V e > 

d2f o + 

6 1 (e) d 2 i 

<t> 2 

dx 2 

<$> 2 dx" 

- 2(k 1 -<j>x) 

rv e: 

_ fo 

6. (e) 

d£i + 1 

<p 

dx 

5 

dx + ‘“ 


- A 0 [6 0 (e)f 0 (x) + 6 1 (e)f 1 (x)+. . . ] 


- K 1 (e)X 1 [6 0 (e)f Q (x) + 6 1 (e)f 1 (x)+. . . ] = 0 


( 2 . 6 ) 


Expanding the coefficient of the second derivative, we have 

1- (k ^-<J>x) 2 = 1 - k 2 + 2<j>k ^x - <j) 2 x 2 . 

We now have to specify some information about the order of magnitude 
-2 

of 1-k . We will assume throughout the rest of the study that 


i - -4 - o(i) . 

k 

By assuming this, we consider the built-in tension to be greater than 
T Q , by a "sensible" factor. The thickness of the boundary layer is 
then obtained by considering the highest derivative terms or 


e6 Q (e) 6 Q (e) 


so 


<P 

<j> = Jz 


and x = 


k~ 1 -x 


The dominant boundary layer equation is given by: 

1 d4f o 1 d2f o 

2 ~4 - ^ 2 ^ ~2 0 ’ 

k dx k dx 
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The boundary conditions are: 


f Q (x) = 0 


df 


8 (x)= 0 


for x = 0 


dx 


Defining d^fg/dx^ = f*(x) and (k 2 -l) = k*^ , we have 
2 * 

d f 0 2 * - 
—3- - k Z f Q (x) = 0 

dx • •• . 


Cancelling, among the two independent solutions of the last equation, 
the exponential growth which is unfitted for the matching process, 
we are left with 


d 2 f 


dx 


0 * r -k x 

~2 r 0 o 6 


A first integration and one boundary condition gives us: 
df n C ir 

ss -T- (1 - e > • 

k 

A second integration and the second boundary conditions gives us: 

C, 


f 0 (x) *2 


o „ * ~ , , -k x s 


(k x - 1 + e ) 


The matching condition with the solution previously found for the 
central portion of the cable enables us to define the constant C . 

We now expand the outer solution, near x = k \ We define the variable 
x' = k d - x. We have for the even modes 

^(A 2 , k _1 - X 1 ) = ^(A 2 , k' 1 ) - se ± ( a 2 , k _1 )x’+... 


and for the odd modes: 
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5? 2 Uq, k -1 - x') =^ 2 (A l, k X ) - ^ & 2 (\ 2 q , k 1 ) X , +. . 


An intermediate limit suitable for the matching near x 
defined by 


k ^ is 


x = 


X 


with x fixed. 

n 


n n(e) 

Taking into account x* = 4 > ( £ ) x = 77. x' , we also have 

Se X , _n_ 


x n n (e) 


-1 


where — 1 ► 00 when e 0 . 


The behavior near x = k takes the form: 


lim 
e -*■ 0 
x fixed 

n 


hp(k ^-x') + V^(e)h^(k 1 -x')+... 

- - V e > £ 0©- 


= 0 


Replacing h^ and f^ by the solution previously found, we have 


lim 
e ->• 0 
x fixed 

n 


D 0^1 ^0’ k ^ " D 0 dx °^i <V k ^ nx n + '"‘ 


+ v 1 (e)h 1 (k 1 ) - V 1 (e) ~ h 1 (k "^nx +. . 

. nx ■ 

* n 

-k t" — 


c 0 f * nx, 


77 


■ 6 ° (e, ?l k \T ' 1 + e 




The matching is then accomplished by taking 

(Aq, k _1 ) = 0 


y = o 
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which is the relation we used before but was still to be proven and 
also 


6 Q (e) = Jz 


* d 
C n = -k D n -p- 
0 0 dx 


^ ttj, k 1 ) 

2 


where 5? 

2 

Now 
V-^(e) has 


stands for both Legendre polynomials. 

that the linear terms have been matched, we can see that 
to be of the same order of magnitude as 6 q(e) so 


v 1 (e) = 71 


and also 


v. k. ) ^2 


With the information obtained through matching, we are now able 
to look at the next order of magnitude for the expansion in the central 
area. Applying a procedure similar to the one used for the rotating 
cantilevered beam, we plan to consider the orthogonality property of 
the modes . 

Let us consider two different modes <J> and <p a defined by the 

u p 

following relations. 

r 


~r <|> IV - [ (1-x 2 ) <(> ' ] ' - A 2 <p =0 
. 2 Y a a a a 

k 


fa - la-x 2 )^]’ - x* 4> B - 0 


and the boundary conditions 

<f> a (±k _1 ) = ^(±k _1 ) = <j>g (±k _i ) = ^’(ik -1 ) = 0. 
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Multiplying the relation defining <J> a by <f>g and integrating over the 


whole span, we have 
, -1 


-1 


-1 


e 

1 


f Y Ob d * - y l( i-* 2> *iJ+e dx - x a J *, 

-k -k -k 


A 


= 0 


The second term can be written: 

,-l 


f [•(i-x 2 )4£] , <t> 0 dx = [(i-x 2 )<i>;]<i> 6 
-k' ;L 


.-1 1 -1 

k k 


-k _1 -k 


J [(l-x 2 )^’ dx. 

i ”1 


The use of the boundary conditions limits the previous expression to 
the last term. Integration by parts and the use of the boundary con- 
ditions transforms the first terms into the following sequences of 


expression 

, -1 


-1 


-1 


f 4> IV| 

0 “ 


♦p dx ' 


- W 


-k 

-1 


I d>'" 

-i J. “ 6 


dx = 


-k 

-1 


! + J x *a*s 

-k -k 


dx 


We thus have finally 


-1 


-1 


_e 

2 


/ <p”<p" dx + / (l-x 2 )^'^' dx = A 2 / <p (p 

I . r crS / , r a Y B . a / , r a T 

-k" 1 -k" 1 -k _1 


,dx 


The symmetry of the left hand side suggest that it is also equal to 


-1 


-k 


V* 


as may be confirmed by a parallel development. Combining those two 
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expressions, we obtain the expected result i.e. 
, -1 


° 


<p a dx = 0 for a 4 3 
p 


The orthogonality relationship between the modes has to be com- 
pleted by the relation normalizing the modes. We are thus free to 
choose the modes such that 



Considering now the asymptotic expansions used for the modes 


<P a (x,e) = h Q (x) + /z h 1 (x) +..., 
the previous relation becomes, for instance: 
k _1 

f hj - i 

-k _i 

k' 1 

f h ohl dx = 0 . 

-k' 1 

Now that the orthogonality has been proved between the different modes 
and also between the different terms of the asymptotic expansion, we 
can come back to the main problem. 

Once more, it should be emphasized that in the expansion used 
for the modes in the orthogonality relationship, a uniformly valid 
expansion should be used, a uniformly valid expansion which takes into 
account the contribution of the boundary layer into the mode. But 
for the problem at hand, the contributions of the boundary layer do 
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not affect the relations found before, because we limit ourselves to 
the expansion up to . 

The differential equation defining h^(x) is obtained by replacing 

2 

the asymptotic expansions valid for X^Ce) and <J> a (x>£) into the differ- 


ential equation of the modes . We thus have : 

t(l-x 2 )h^] + X 2 h 1 = j-X 2 h Q if ^(e) = 

K,(e) 

0 if — = -> 0 

The relation of orthogonality between h^ and h^ found before, limits 
the possibility to the case where K.(e) = -St . The equation we have 
to look at is now given by 


[(l-x 2 )hj/ + X 2 h x = -X 2 h Q . 


( 2 . 7 ) 


From this equation, we notice that an even correction h^ corres- 
ponds to an even function h^ and reciprocally an odd correction to an 
odd function. With the remark noted, we can now use the relation 
obtained with the orthogonality property. The relation 
gives us 


X-1 h o 0*- 1 


-1 


/. 

-k 


B 2 x\ (Xq, x)d X 
-1 2 


= 1 


or 





,x)dx 
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Now, considering the equation defining , Equation (7), and multi- 
plying this relation by h^ and integrating over the whole span we have 


-1 


i - 1 


-1 


-1 


[ (l-x 2 )hj] h Q dx + A 2 


f_ h l h o dx ■ - 4 J h 0 dx “ A l 

-k -k 


where the normalization relation has been used. 

An integration by parts and the use of the relation 


£ 


^h^hg dx = 0 gives us 


[(l-x 2 )h']h 0 


k" 1 k' 1 


tv n. - 

- j a-^hj »• = -xj 

-k -k 1 


The use of the boundary condition and a new integration by parts gives 
us 


-(l-x 2 )h^ h 1 


-1 -1 
k X kj- 

+ I h l 

-k -k x 


t(l-x 2 )h^] dx = - A 2 


Considering now the differential equation defining hg. Equation 

(5) , and using the orthogonality between h_ and h , we obtain 

k' 1 


Aj = (l-x 2 )h^h 1 


-1 


We mentioned before that to even modes corresponds even corrections 
and vice-versa. As a result the expression h^ h^ takes on different 
signs when evaluated at the end points k X and -k X and we are thus 
left with 


A 2 = 2(l-k’ 2 )h^ ) (k" 1 )h 1 (k X ) . 
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Considering now the information we found in the matching near 

-1 


x = k we have that 


hf)(k ^ D 0 dx ^1^0* k ^ 

2 

hi<k ) = - 7* dx ^0* k ^ 

k k 2 


-U 


We thus have 


A l = 2 




The final expression for A^ is then 


2 v 

Xi‘ 2j S 


- 'df »o- k_1 >’ 

1 L 2 


-1 


/4 


(Aq, x) dx 


2 2 

The values found for Aq and A^ are also those values which 

guarantee the orthogonality between different modes. This derivation 

is less straightforward and is presented in Appendix I, so to lighten 

2 

the presentation. Now that the value of A^ has been found, the 
determination of h^(x) is completely determined with the help of the 
boundary conditions we obtained in the matching process. 

To complete the study of the present problem, we will now look 
at the next term of the asymptotic expansion used in the boundary 
layer. The matching near x = k A suggests as the order of magnitude 
for 6^(e) the value e. The differential equation defining f^ is given 
by considering the terms of order of magnitude 
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e6 1 (e) 

4> 4 


£ O 

~2 =£ 

£ 


in the boundary layer expansion, Equation (6) , or 

1 d4f l „-2. d2f l , k -l- d % , -1 df 0 „ 

, k 2 di 4 di 2 dS 2 dx 


The expression previously found for f^ is 

C 


* 

0 -fir — 1c x 

f q (x) = —^2 tk x - 1 + e ] 
k 


where Cq is a known constant. We thus have: 


d4f l 2 d2f l 

— r - -1) — 2 “ 2k C 0 

dx dx 


* 

-k x 


~ -k x , 1 e 

x e + — - * 

k k 


It would be instructive to compute the complete solution valid 
for f p . To reach this objective let us define 

* 

— a f 

,~2 *1 * 
dx 

*2 2 

We have, taking into account the definition of k = k -1, 


d 2 f* 

dx 2 


*2 * 

k f, 


= 2k C r 


* ,, 

,* ~ . -k x 

~ -k x . 1 e 

x e + -* - ~ 

k k 


The particular solution f pp is taken under the form 

£ ^ £ 

* -k x, ~2 -k x. 

f lp = 2k C Q [a 0 + a 1 x e + a 2 x e ] . 


Introducing this development into the equation defining f p , we obtain 
the following three relations: 


a 0 *3 * “2 


a„ = - 


4k 


and a p = 


4k 


*2 
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The total solution is now written as: 


d 2 f *~ 

f * = lil =c e “k x + 2k C 

f l ,~2 C 1 e + 2k C 0 
dx 


1 , _1 ~ -k*x ~2 -k 

JLO ' JL r\ X 0 “ r\ X 6 

L k 3 4 k* 2 4 k 2 


The boundary condition for f^, are given by: 


df l 

-j=o— (x) = f (x) =0 for x = 0. 
dx 1 


Integration of the previous equations gives us: 


df l C 1 -k*J 2k C 0 - , k C 0 -k*J 

. e -—*3-* + _e 

k k 2k 


x 1 

* *2 

-k k 


k C n *- 
0 -k x 
* e 

2k 


x 2 2x _ 2 _ 

L-k* " k * 2 ' k* 3 J 


+ c. 


The use of one boundary condition transforms this expression into 

df i 2 k C n kC n kC. , , *~ 

1 „ -k x 0 ~ . 0 -k x . 0 ~2 -k x . 

— = " c 2 e " ~* 3 ~ x + —£T x e + -~*2 X 6 

dx k J 2k 2k 


Integrating once, more we have: 


C * 

, z -k x 


kC n , k C n * 
— 0 x 2 4 - ^ -" k x 


*3 


2k 


*3 


-k 


1 

*2 


k c n 

. U -k x 
H — e 


2k 


*2 


-k 


2x 
, *2 


*3 

k 


+ C2X + 


The boundary condition gives us 

k C 


S k^ 

* *5 

k 2k 


if + C 3 


0 or = 3 


k G 0 * 
— —■ - kC. 
, *4 3 
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So we finally have 


k C 0 ~2 -k x 3k C 0 ~ -k x -k x k C 0 ~2 

f l = ~ “ *3 x e - *4 x e - C 3 e - -jpr- x 

2k 2k J k 


3k C 0 , \ ~ 
k *4 - k c 3 ) x + c 3 • 


We can see that in the previous relation, only one constant is left 
unknown, C^, and this last constant will now be determined through the 
matching process ! ! 

The matching near x = k ^ is given by the following expression 


lim h (k 1- x '^ + v i (e)h. (k 1 - x')+. .. 

e -*■ 0 


fixed 

n 


/ nx \ /'nx \ I 


Some of the terms appearing in the previous expressions have already 


been matched, so let us consider the remaining ones; we then have 


lim 
e ■> 0 
x fixed 

n 


_ d « ,.2 ,-i v 1 n . . 1/2 

• • *Dq 2 °^i^0 yk ^ 2 ^ |:h^(k 

dx o _ L 

_ k * 

] c ^ 

- Jl -2j -i + a ^ 


k (k_1) - h V k " 1,r " 

*- -1 




— £ _ L -J- 

*3 e 

k 


- ~rr ~ - k C n + C, 

k* 4 °1 si 3 


k c. x n 
o n 

*3 e 


2 -k 



1 


Ji 

3kC o^ " 

c 

Ji 


' 2k* 4 6 

> 


* nx 
-k*-a 



- C„ e ^ 

+... 1 
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The matching is obtained by the satisfaction of the following 
equalities 

1 D £— v a 2 k' 1 ) + — ° - n 

2 u ,2 Z ! U 0’ k ' + *3 " ° 

dx 2 k 

-1 C 0 
h l< k - > 


3 )-» 


The last equality defines the value of and the last constant to be 
determined is thus given by 


c 3 ■ 

J k 


a 1 3k c n 

s h i< k > + ~*r 


The second equality was already found before and used. We are thus 
left with the last equality or replacing Cq by its expression 


2 D 0 2 ^l (X 0’ k “ *3 k D 0 dx ^i (A o ,k * 0 

dx o k 9 


This last relation can also be written 

,2 


c * 2 _sL 

dx 


2 , -1 


2 ,-l. 


2 «l«0' k > - 2k £ ) ■ 0 • 


and this last relation represents an identity as we can see from the 
differential equation defining h^ or 

(l-x 2 )h^ - 2x h' Q + A 2 h Q = 0 . 


2-85 




which is the last identity. 

With this, we have completed the study of the modal analysis for 
the transverse vibration of a rotating cable of small flexural rigidity 
clamped at both ends. This last problem has been introduced mainly to 
become more familiar with those types of problems and gain more insight 
into the physics of the problem. The experience is mainly to be used 

in the next step of our study which will be dealing with a two dimen- 

sional approach for similar problems. 

We conclude this chapter by summarizing some of the results of 
this last section. The final expression for the eigenfrequencies is 
reported as 

2 1 n 2 .2 1 n 2,,2 ,/ Ty 2 , 

oj = — A = t Si. (a. + /e a. +...) 

a 2 a 2 0 1 

2 

where Aq is the solution of the transcendental equations 
^(Ag.k -1 ) = 0 and S? 2 (AQ,k -1 ) = 0 
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The determination of the constant k 2 is obtained by recalling the 
definition T = k 2 T Q , where T Q - V ^ 2 L 2 /2. 
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. CHAPTER 3 


DYNAMICS OF ROTATING THIN PLATES AND MEMBRANES 

3.1 Introduction 

In this chapter, we concentrate our attention on the problem of 
a dynamical study of two elementary planar structures characterized by 
a dominant effect of the forces induced by spin over the flexural 
rigidity. 

We focus our attention mainly on two problems. In the first one, 
a modal analysis is done for a rotating very thin circular plate or 
membrane with some flexural rigidity, clamped along its edges. This 
type of study is oriented towards the possible use of large spinning 
membrane-like disks as optical and radar reflectors for space vehicles. 
This last use of spinning membranes has led to a renewed interest in 
the problem of calculating the transverse vibrations of an elastic disk 
rotating at a constant speed — a problem which traditionally has been 
studied in connection with gas and steam turbines. In the second part 
of this chapter, the modal analysis of a free spinning membrane is 
accomplished. In both of those studies, only vibrations usually 
referred as "out-of-plane" are considered. 

3.2 Dynamics of a Taut Rotating Membrane in a Circular Ring 

The deformable element consists of a large flexible circular 
membrane undergoing transverse vibrations, while rotating with an 
angular velocity ft with respect to the axis of symmetry, defined when 
the structure is at equilibrium. 
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The following equations of the motion are given in a reference 
frame located at the center of mass of the total system and rotating 
with the system at the constant angular velocity ft. The direction of 
the E, axis will be chosen in such a way that there will be no need for 
a phase angle in the equation of the motion. We are free to choose 
such an axis due to the symmetry of the particular structure under 
investigation. The physical characteristics of the membrane will be 
assumed constant throughout the whole membrane. 

It must be emphasized here that the effect of the spin on the 
membrane is a destabilizing effect: the rotation not only introduces 

a decrease in the stiffness but also instability can occur due to the 
creation of a state of compression into the membrane. To take care of 
this last problem, we introduce a built-in tension — a similar pro- 
cedure was used for the rotating cable — such that the membrane, when 
rotating, is everywhere under tension. 



Figure 8. Rotating Taut Membrane. 

In the study, we will disregard the stresses, which occur on 
account of the mutual pressure of horizontal layers of the disk and, 
consequently, the strains in the direction of its thickness; the 
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problem is thus simplified and the generalized Hooke's laws give us the 
following relations: 


°p ' 777 [e P + 

°e ‘ 7^72 te e + V 


where E stands for the modulus of elasticity, 

V stands for the Poisson modulus 

e Q and e stand for the tangential and radial strains 

O p 

respectively 

(Jq and 0^ stand for the tangential and radial stresses 
respectively. 


Also, because the thickness of the disk is small in comparison with its 
radius, the variation of radial and tangential stresses over the thick- 
ness can be neglected. Also, a hypothesis analogous to the hypothesis 
of plane section in a rod will be used. A linear element within the 
disk perpendicular to the middle plane remains straight and normal to 
this plane after its deflection during bending. All those previous 
assumptions are part of the approximate theory of thin plates. 

In this study, the problem of the stretching of the thin disk is 
solved first and the dynamical equations for the transverse vibrations 
are then derived as deviations with respect to the steady state. This 
last remark constitutes a characteristic of continuous modeling and 
will be materialized in the following pages. We also assume that the 
stretching of the disk does not alter the constancy of its physical 
characteristics . 
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The stress distribution in the membrane must satisfy the follow- 
ing expression, as an equilibrium condition for the forces acting along 
the radius . V°0 


°d 

Figure 9. Distribution of Stresses. 

^ (pa p ) -a e + y fl 2 p 2 =0 (3.1) 

where y is the mass per unit volume of the material of the disk. 

Because of the symmetry, there is ho dependence in the angular variable. 
This relationship and the following results from classical thin plate 
theory can be found in Reference 20. 

For a circular disk, the strain components, in the case of sym- 
metry, are well known and given by 




and 



where u represents the radial displacement. The generalized Hooke's 
laws and the expressions for the strain components give us: 



Combining those expressions with the equilibrium condition (3.1), we 
obtain the equation for the radial displacement of the disk u.: 
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2 d u . du 


1-v 


,2 3 


2 + p — - u = - y ^ p . 


dp 


The general solution of this equation is 


1 

U = E 


C 2 

(l-v)Cp - (1+v) -i - ~~ y Q 2 
P ° 


-•] 


where C and are arbitrary constants. The corresponding stress com- 
ponents are now found from: 

r 

1 3+V 2 2 

° P = c + "7 ' IT e 

p 


_ „ 1 1+3V n 2 2 

a 0 = C - ~ - 1 y p ‘ 

P 


Because we are dealing with a complete disk, we must take = 0 to 
have u=0 at the center. 

In the determination of the constant C, we have to be very 
cautious . The reason is coming from the interpretation of the boundary 
condition for the problem 

u = 0 when p = a. 

This last boundary condition is, indeed, only valid after the membrane 
has been stretched statically by the built-in tension T. As a result, 
the constant C contains two parts, one coming from the built-in tension 
and the other from the satisfaction of the previous boundary condition 
for the radial displacement originated by the radial forces induced by 
spin, only. The last reasoning is perfectly valid if we consider that 
the generalized Hooke's laws represent a linear dependence between 
stress and strain. We know that in a circular membrane under a 
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constant tension T, the stress components are the same and equal to the 
applied tension. We, thus, have 

C = T + C’ , 


where T is the constant built-in tension. The determination of the 
constant C' is now determined by using the previous boundary condition, 
but by considering in the constant C appearing in the expression of 
the radial displacement only the contribution of C*. The total dis- 
placement is then obtained by the summation of the displacement 
produced by the built-in tension and the displacement we just derived. 
So 

2 

(l-v)C'a = y ft 2 a 3 


and 


C' - y n 2 a 2 


The final expressions for the radial and tangential stresses are given 
by 


a 

P 

a e 


T + 
T + 


1+v 

8 

1+v 

8 


n 2 2 
y fi a 

n 2 2 
y a 


3+v 

8 

l+3v 

8 


n 2 2 
y n p 

n 2 2 

■ y Q p 


(3.2) 

(3.3) 


One comment should be added here. Namely, we assumed before that the 
built-in tension must be such that a positive tension is present every- 
where in the disk or 


= T + 


1+V 

8 


o 2 2 
y a 


3+v 


n 2 2 
y 0, p 


> 0 


for every p such that 0 _< p a, or 
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„ , 1+v n 2 2 3+v .22 

T + — — y a - — r— y ft a .>0, 


which gives us 


T > 


n 2 2 
y a 


4 


For simplicity, we will write the tension T under the form 

. 2 n 2 2 
T _ k y Q a 

4 


where the positive constant k is greater than one. With this last 
remark, we have 

1 + 2k 2 + V a 2 2 3+v ^2 2 

°p = ' — 8" y a _ ~8~ V ^ P 


1 + 2k + v n 2 2 1+3V ^2 2 

= y n a - — g- y n p 


We thus have illustrated here the fact mentioned several times earlier 
that the equation of the transverse vibration for a continuous model 
represents in fact the perturbation of the system with respect to its 
steady-state deformation. 

Now, that the distribution of the radial and tangential stresses 
is known, the dynamic equation for the transverse vibration of a cir- 
cular plate subject to radial and tangential stresses is given in 
general form by: (see Reference 20) 


2 4 

E h V w 

. £ 

. 12(l-v 2 ) 





= y 



wher e 


h 



represents the thickness of the disk and 

is the Laplacian operator which expressed in polar 

coordinates represents 


2-95 



v 2 ,-ii + I -» + -I 

P 3p 2 P 3 P p 2 30 2 


The boundary conditions for this problem are given by expressing the 

clamping condition along the edges or 

w = 0 | when p = a 

3w 


3p 


= 0 


Let us now introduce the dimensionless quantity r, defined by the 
relation: 

p = ar. 


The previous expression becomes then: 
Eh 2 


12a 4 (l 


„4 1 3 3w 

— r V w + —z~ 0 r 7— 

2. 2 3r r 3r 

-V ) a r L J 


+ 


2 2 
a r 


a 2 

o w 

30 Z 


= y 


a 2 w 

3t 2 


(3.4) 


where 


and 


2 2 

°r = [(1 + 2k2 + v)-(3+v)r 2 ] 


2 2 

o 0 = [(1 + 2k 2 + v)-(l+3v)r 2 ] . 


Dividing through by y fl we obtain: 

Eh 2 a 4 _1 _3 

2 4/, .2, V w 8r 3r 


12yfi a (1-v ) 


[(1 + 2k 2 + v) r - (3+v)r 3 ] |^ 


+ [ ( l+2k 2 +v) - ( l+3v) r 2 ] 

8r 30 ST 3t 


This last equation is appropriate for the use of separation of vari- 
ables, providing the initial conditions are also separable. Let us 
then use for the solution an expression of the form: 
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i(a) t + m0) 
. , N a,m 

w = 4 (r)e 
a»m 


This last expression is completely general and there is no need for a 
phase angle for the reason explained earlier. 

The equation of main interest for our purpose is given by the 


differential equation defining (p (r) : 

a,m 


Eh 4 (l+2k +v) _d 

12(l-v 2 )aV ” “> m Sr dr 


r 1 - 


3+v 


l+2k 2 +v 


2\ ^a.m 


dr 


(l+2k 2 +v) mi 

8 2 
r 


1 - 


1+3V 2 

r 


l+2k 2 +v 


(0 


<f> 


a,m 


a.m 

« 2 




where oo^ m represents the eigenfrequency of the vibration 


m 


is the linear differential operator defined by 

n 2 = d 2 - I 
-m 2 r dr 2 
dr r 


We now define the following quantities 

2 


c = 


2 E h 


3(l-v 2 ) (l+2k 2 +v)a\iQ 2 


8 w 


q,m 


a,m (l+2k 2 +v)Q" 


k 2 = 3+v 

K 1 2 
l+2k +V 


, 2 l+3v 

^2 2 

l+2k +v 


The final expression for the differential equation for (J> a ^(r) is thus 


given by 
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0 

(3.5) 


-eV <j> 

tn ' , 


m a,m 


r 


_d 

dr 


tr(l-kjr 2 ) 


d<{> 


dr 


^ [l-k 2 r 2 ]<p +X 2 <p 

2 2 JT a,m a,m a,m 


with the following boundary conditions: 


<p 

a,m 


(r) 


= 0 


dcf> 


a,m 


dr 


(r) = 0 


y 


when r=l 


Here, again, we notice that the small flexural rigidity of the 

disk combined with the high spin rate produces the small parameter £. 

2 2 

We also notice that the parameters k^ and k -2 represent positive quanti- 
ties less than one. 

The problem has now been formulated in such a way that we are in 
a position where a procedure similar to the procedure used before can 
be used. 


Expansion Valid in the Central Region 

When r ^ 0, the expansion valid in the central area takes the 

form 


m ( r > £ ) = V r > + V i( e )h^(r) + v 2 (e)h 2 (r) + . . . 

A similar asymptotic expansion is also used for the eigenvalue or 


X 2 (e) = X 2 + k. (e)X 2 + kAz)X 2 0 +... 
a,m u l l 2 2 


Once more, the subscript a,m will be momentarily dropped in order to 
shorten the notation. It will be explicitly used later on when needed 
to differentiate the modes. 

2 

Introducing the asymptotic expansion valid for (J) and X 

a y HI 0 1 y IQ 

into the differential equation (3.5), we obtain the following expression: 
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- £ V m (h 0+ v 1 (e)h 1 +...) 


+ i / 
r dr 


n 2 2 . dh 0 

r(1 " k l r > 17 


, , . 1 d 

+ v. (e) - 3— 

1 r dr 


,, .2 2. dh i 

r(1 - k l r > "drJ 


2 2 

m r , . 2 2,, .. . m r. 1 2 2,, /v 2 , .,2, . 

- — [l-k 2 r ]h Q - v 1 (e) — [l-k 2 r Ih^ (X Q +K 1 (e) . .) 

r r 

' (h Q (r) + v 1 (e)h 1 +...) =0 

The differential equation defining hg(r) is then given by: 

1 d I , 2 2, dh 0 1 m 2 M ,2 2, ,2 , _ n 

r dr ^ 1-k l r ^ dr J “ ^2 ^ " k 2 r -* h 0 + *0 h 0 0 * 

This last expression can be written, when we introduce the new variable 


k i r = r i 


as 


r l dr l 


„ 2. dh ol m 2 .,,2 2,,, , X 0 

r i <1 - r i ) 33" ' ~7 tl_k 3 r i ] h o + 72 
ljr, k. 


h 0 = 0 


(3.6) 


where k? stands for ~ = ~~~ V ~ . 

J . / JTV 

k i 

Let us write down the solution of the previous equation under 
the form 

h 0 = r l H 0 ’ 

By substituting the previous expression into the equation defining 


h n (3.6) we can eliminate in the coefficient of the terms in r 


n-2 


0 

by taking for n the solution of 


2 2 - n 

n - m = 0 

or 


n = ± m . 
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The solution n = -m is to be eliminated for the singularity that such 

a solution introduces at the origin; we thus have 
, m .. 

h o ■ r l H 0 • 

The differential equation defining is then obtained 

, 2. d H 0 : 1 . 2. dH 0 , dH 0 2m . 2. dH 0 

<1 - r i ) TT * -T <1 - r i ) £7 - 2r i d77 + (1 - r i> d77 

dr, 1 1 .11 1 


, _ . 2 2.2 0 
-l 2m + m - m k_ - — — 

J k 2 
fc l 


H o = ° 


or 

2 

9 d H n i ? dH n / ') ')') “n 

(l-r 1 ) - — 2 ‘ + — t(l+2m)-(3+2m)r 1 ] ^ (2nri-m -m k 3 ~ “ J H q = 0 • 

dr i 1 1 \ k l; 

(3.7) 

To bring this last equation to a form more familiar, a new change of 

variable is recommended and we define 
2 

r ] _ = x. r 

We should emphasize here that x can only take on positive values lying 

2 2 

in the interval between 0 and r^ max = k-^ which is less than 1. With 
this last change of variables , the previous equation can be written 
under the form - 

dH 0 i / 9 2 2 A 0 

x(l-x) — — + [(nri-l)-(nrt-2)x] — — - ( 2mfm -m k^ 2 J^q~ 0*8) 

dx \ k i/ 

This last expression can be reduced to the hypergeometric equation 
which is written under the general form 


d 2 H, 


x(l-x) 


0 


dH 


dx 


+ [ y-(a+g+l)x] 


0 


d X - aBH o- 0 
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So we have to solve the following substitutions: 

Y = 1 + m 

a + 3 = 1 + m 2 

a3 = j \2m + m 2 -m 2 k 2 - -j 

\ k I 

The range of value for x lies within the range of convergence of the 
hypergeometric function. 

For the previous equation, we know that the roots of the indicial 
equation at the regular singular point x=0, are given by 0 and (1-y) . 

As a result, one of the solutions of the previous equation can be com- 
puted by considering 


H 


10 


= X) a x n 


n=0 


A direct substitution into the equation yields the identity 


w w — 

J'j n(n-l)a n x n 1 - n(n-l)a n x n + 


n=0 


n=0 


n=0 


ny a x 
' n 


n-1 


- VI n(a+3+l)a x - ^ ] a3 a x = 0. 
n=0 n=0 


We then group the corresponding summations and find 
00 <» 

^ n(y+n-l)a x 11 ^ [n 2 + n(a+3)+ a3]a x n = 0 

r\ n a n 


n=0 


n=0 


Now this may be written 


00 00 

^ n(y+h-l)a n x n ^ (os+n-1) (3+n-l) a .x 11 ^ = 0 

n=l n n=0 n ~ 
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From there, we find that the recursive relationship for all determin- 
ation of the a 's are 
n 


(g+n-1) (B+n-1) 
a n (y+n-1) • n a n-l 


for n > 1, with a^ arbitrary. At this stage, we must insist on the 
fact that y has to be different from zero or any negative integer to 
make sure that the previous recurrence is defined, but 

Y = 1 + m which represents always a positive 

quantity. 

Finally, the solution of the previous equation is written with 
the help of the hypergeometric function 


F(a,e, Y ;x) = l + ^ 


x a(arfl) g(g-H) 

1! + y(Y+D 



a(ori-l) . . . (ori-n-1) g(g+l) . . . (3+n-l) x 11 

Y(Y +1 ) • • • (Y+n-1) nj 

under the form H^q = agF(a,j3 ,y;x) . 

The solution we just developed represents one solution of the 
previous equation. This last solution is convergent for |x| < 1, 
which is our case of interest. 

In our particular case, it turns out that the difference between 
the roots of the indicial equation at the point x = 0 is zero or a 
positive integer. As a result, the second solution generally used for 
the hyper geometric equation, more explicitly 

x 1 ^ F(a-Y+1, $-Y + l> 2-y;x) 

is not linearly independent of the solution previously derived. How- 
ever , it can be proved that the second linearly independent solution 
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is to be canceled. The derivation of the second solution is quite 
lengthy and is presented in Appendix II in order to lighten the text . 
The final expression for hg(r) is then 

h Q (r) = a Q k“ r m F(a,3, 1+m; k 2 r 2 ) 


or 


where 


h Q (r) = b Q r m F(a,3, 1+m; 


3+v 2' 

r i 


l+2k 2 +v 


ot+3 = 1+TO 


2 2,2 


a3 = T I 2m + m -m . k 0 - 


0 


3 k 2 
R 1 


and 


if we write 


8 A 


,2 

0 


(3+v)S^ 


i 2 = A' 2 + k.(c)X1 2 + K 9 (e)AI 2 +... 
a,m u 1 1 2 z 


Boundary Laver Near r=l 

We now have to study the boundary-layer near r=l in order to 
pursue the problem. 

The separation of variables being valid for the partial 
differential equation for the problem, the boundary layer problem is 
limited to radius dependence of the modes . 

To study the boundary-layer expansion, another boundary layer 
coordinate must be selected, so that the higher derivative terms 
dominate the equations and are of the same order of magnitude. Let 

•s* ] — y 

r = where (p(e) 0 when e -»■ 0. 
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The corresponding asymptotic expansion valid near r=l is taken as 

Vm (? ’ £) = V £) f 0 (f) + S 1 ( e ) f 1 (?) '+... 

where the 6^(0) constitute an asymptotic sequence. 

2 

The linear differential operator becomes now 


V 2 = d 


m 


m 4> 2 df 2 (1— ) (j)df (l-<J>r) 2 


The basic differential equation ( 3 . 5 ) becomes with the previous substi- 
tutions: 


- £ 


m 


4> 2 d? 2 (l-$r) <})dr (I-*) 5 ?) 2 


2 , 
m < j > 


(l-<J>r ) 2 


(l-<f>r) 4>dr 


d 2 <J) 


2 ~2 
4> dr 


q.m 1 _ 


dcj> 


q.m 


(l-<j>r) (J)dr 

r\ dd> 

2- q.m 


0H>?) [l-k^ (!-<}>?) 1 


(j)dr 


m 


(l- 4 >r) 


- [l-k^a-cf*) 2 ]^ + A 2 4 j = 0 

2 2 T q,m Y q,m 


( 3 . 9 ) 


or 


- e 


i d ^< j > n 

1 q,m _ 2_ 

,4 a-* 

<p dr 


3 2 

d 4» d <J> 7 , 

q,m a ,m / 1 


(l-4>r) 4>^dr 


+ 


2 m 


<}> 2 df 2 \ (l-<J>r) 2 (1-4)?) 2 ; 


d(f) 


q,m 


1 . 2 m 

+ - 


4>df( \(l-4>r) 3 (1-cJ.f) 


n I + < j ) 

3 / Y a,m 


m 


4 m 


(l-<j>r) 4 (l-(J>r) 4 


(l- 4 >r) 4 ) 


(l-4>r) fl-rk 2 (l-<J>r ) 2 \ — — (l-k 2 (l-<f>?) 2 ) 

' ' 4>dr ' ' 


d(f> 


q.m 


dr 


+ 2k 2 (l-<j)r) 2 


— 2 r (l - k 2 (l-<j)r) 2 )4> + A 2 6 = 0 

( l -< J ) r ) 2 ' 2 ’ a,m a,m 


( 3 . 9 ’) 
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The thickness of the boundary layer is determined by balancing the 

2 

higher derivative terms or, if we assume that 1-k^ is not of an order 
of magnitude of £, by considering: 

e ~ 1 

x4 ~ * 2 
<P <P 

or 

2 

e = (p and <J> = /e 

2 

The previous remark, about 1-k^, is considered now 


1 _ k 2 _ i _ ..3+v = 2 (k - 1 ) 

12 2 
l+2k +v l+2k +v 


The previous assumption is of the same nature as the assumption we 

made while studying the taut rotating cable. We want the parameter k* 

2 

to differ sufficiently enough from one, so that the difference 1-k^ 
differs from zero by an order of magnitude of £ . 


Replacing in the last development 4> a m (r,E) by its expression 

defined previously, we obtain for the differential equations defining 

f_, the dominant part in the development (3.9') or 

u 4 2 

d f 0 2 d f 0 
- + (1-kf) = 0 


d? 4 


1 ,~2 
dr 


The boundary conditions for the problem are 


f 0 <?) = 0 


df 


dr 


| (?) = 0 


when r = 0 


Defining 


A> * 

. ~2 t 0 
dr 
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and considering the fact that l-k 1 


as a solution for f 


0 


* • 

f 0 = C i 6 


-vVkj r 


is a positive quantity, we obtain 


where the exponential growth has already been disregarded for its 
problem in the matching process . 

Integrating the previous expression twice, and taking into 
account the boundary conditions for f^, we obtain finally: 



-v'l-k 


1 + ^1-k 2 r 



where the positive sign is always chosen in the previous square roots . 

We are now in a position where we can try a matching between the expan 

sion valid in the central area and the boundary layer expansion. The 

2 

equation defining the eigenvalue A^ should come out of this matching. 

We need to expand the solution valid in the central area near 
r=l. We define the variable 


r ' = 1-r 


and we have 


dh. 




where 


h 0 (X^ 2 ,l) = b Q F(a,e, 1+mj 


3+v 


l+2k 2 +v 


and 
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dh 


TT (A’ 2 ,l) = b 0 m F(a,3,l+m; — ^ — 

\ l+2k +v . 


. , 2(3+v) dF , „ , . 

0 ^ (a.B.Hm; x) 

l+2k + V 


x = 


3+v 


l+2k +V 


For matching, an intermediate limit suitable for matching near r=l is 
given by r^ fixed: 


r n n(e) 


where r' is always positive. 


What we want in r 00 when e -> 0 and this is reached if 
—2 >- oo when . e ->- 0 

/e 


The behavior near r=l takes the form: 
r 


lim 
£ 0 

r fixed 

n 


2 ^^0 2 

h 0 (X 0 » 1} " “dr (X 0 2,i) V-' 


+ v i (e)h 1 (l)+. ..-6 0 (e) 


A-k 2 




- 1 


' 6 i <E)£ i&) + - 


The matching is accomplished if we take 


h 0 (A^ 2 ,l) = 0 


dh 


6 Q (e) = and - (A£ 2 ,l) = C ± 


The first equation is nothing else than the equation defining the 
2 

eigenvalue Aq of the problem or 
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Fla, 1+m - a, 1 + m; 


3+V 

l+2k 2 +v 


)-°- 


The second equation gives us 


c = -h dF(a,e,l+m;x) 

1 °l + 2k 2 « d * 


X = 3+v/l+2k +V 


Coming from the matching, we also have that 


v 1 (e) = /z 


and 


^(1) + 




= o 


or 


h, (1) = 


_0 2(3+y) dF 


1 / 2 

/^2 1+2k +V 


2 . ... dx 


(a, 3, 1+m; — ) 

\ l+2k +v / 


This information about h^ enables us to obtain the differential 
equation defining h^. Considering the development of the basic 
differential equation (3.5), we have: 


Id r , , 2 2, dh ll m 2 

7 37 ra-ki r ) -57 - - 

L Jr 


[l-k 2 r 2 ]h 1 


X 0 h l 


- - x J h o 


(3.10) 


This last expression implies also the choice of Jz for K^(e) . This 
value has been considered for the same reason as in earlier development. 
There exists an orthogonality relationship between the modes — shown 
later — and this last constraint suggests as an order of magnitude of 
K^(e) , the value /z . We see, indeed, that for A^ » 7z , the solution 
for h^ would be the same as the one found for hQ and the orthogonality 
would be violated. 
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In the differential equation (3.10), the only freedom left, to 

2 

define the problem completely, lies in the determination of and this 
constitutes our next task. 


Orthogonality Relationship Between the Modes 
The complete expression for the modes is given by 


<J> (r) cosm9 . 

a,m 

Let us now prove that the orthogonality relationship between two 
different modes is given by 



r)cosm6 cosnB drdB = 0 


This last relationship is straightforward if we consider two modes for 
which angular dependences are different. The direct integration with 
respect to dB, for m ^ n gives us directly 
2 tt 

cosmB cosnB dB = 0 



As a result, the case of main interest is given for two modes having 
the same angular dependence i.e. with m=n. Let us now prove that 
1 

f X <p _(r)<f> R m ( r ) dr = 0 • 

J Q ct , m P » m 

The differential equations defining <J> and <p a are respectively 

u,m p,m 


-.4 . 1 d T , 2 2 V ^a,m| m r , , 2 2, , 

7 4> + — — r(l-k, r ) — » ■ r [l-k„r ]d> 

m Y a,m r dr (_ 1 dr J ^2 2 JY a, 

7 4 . 1 d f y, . 2 2^ d< ! > 3,ml m 2 , 2 2,, 

7 m %,m + r dr [ r(1 k l r ^ dr r 2 [1_k 2 r ^0, 


2 2 2 
-kS ]<j> + X <j> =0 

2 a,m a,m a,m 


+ a r Jr m =0 

.m p y m p y m * 
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We now have to distinguish between the cases where m=0 and m^O. 

V 

Case A: m ^ 0 


The boundary conditions are given by m = ® 


d(j) 


^ = 0 


dr 


V for r=l 


and similar boundary conditions for (j)^ . 


Developing the expression of we obtain: 

, d 4 (J) 0 d 3 cj> „ 2 .. d 2 <{> 2 dcj> 

„4 , q,m _2 T q,m _ 2m +1 q,m 2m +1 q,m 

m ^q,m . 4 + r , 3 2 ,2 3 dr 

dr dr r dr r 

4 . 2 
, m -4m , 

H : ; ® • . 

4 - q,m 


Taking the differential equation defining (j)^ ^ and multiplying this 
relation by r 4^ m > an integration from 0 to 1 gives us 


-e / r V 4 (p <p a dr + 
J m q,m 3>m 


f 0 dr [ r(1 ^ r2) ■3H*6,m dr 

1 2 1 

J ' — [l-k«r 2 ]cj> <j>„ dr + /A 2 (J) <J>„ dr=0 

n r 2 q,m 3,m q,m q,m p,m 


The first term can be developed as 

d 4 cj) r 

, r O. Jr J r = / r 

'0 


£ 


q,m 


0 dr 


Qa Jr + 
p,m 


1 d 3 ^ 


/ 


<Po „ dr 


1 a s 2 * 

2,i d <p 

+1 q,m 


f 2S_ 
"* / 0 r 
1 

£ 


0 dr 
1 


q,m 

3 Y 3,m 
2 j.i d<j> 


dr 


/ O 2. 1 a? 

2m +1 q,m j. , 

— — <p dr 

2 dr p ,m 

u r 


4 , 2 

+ / — — 7” — <j> 4>_ dr 

/q J Y a,m Y g,m 


Integrating by parts and using the boundary condition, the right hand 
side can be written 
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-L 


1 3 

d <j> d<J) Q 

ShSL . ft . m dr + 


1 d 3 (J) 


"0 dr" 


dr 


/ 


a.m , , 

<5> fi dr 


0 dr 


3 ^.in 


1 

/ 


dr 


4 ■ 2 

, m -4m , 

+ I ; q> <p„ dr 

T /"V m T Q Tri 


'0 Ct,m Y 6,m 

One remark should be added here. For the case where m^O, the nature 
of the solution is such that there will always exist for those modes, 
a nodal diameter with a zero transverse displacement. But the exist- 
ence of a nodal diameter of zero displacement implies, that the dis- 
placement of the origin should be zero and we thus have for the case 
where m^O 

<f> _ = 0 = ' for r - 0 . 

a,m p ,m 

An integration by parts performed on the three first terms, combined 
with the use of the previous remark transforms the last expression into 


12 2 

d 4> d z <f>_ , 

'‘-f dr + <2« 2 +l) 

r 0 dr dr 


r d 

f r — % 


( I. , \ 

\r dr ^3»m/ 


+ (2m 2 +l) f - — dr 
r dr dr 


+ (m 4 -4m 2 ) / — y <f> cj> Q dr . 
* / 0 r 4 a > m 


A similar derivation could be done starting with the differential 

equation defining 4> R , and subtracting this similar expression from 

p ,m 

the expression we just found, we obtain 
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-e [2m + 1] 


_1 ^^a.m , 
r dr ^3 ,m 


1 d< ^ft n , 

- — — <j> 
r dr T a,m 


+ ( X 2 - X 2 ) f r 4> (j) =0 

y a,m 3 >™/ Jq a,m 3,m 


Considering that <p and <p a are analytic in the neighborhood of the 

a,m p ,m 

origin, we can expand <p^ m and cj>g m in the following way 


dcj) 

cj> (r) = (J) (o) + 

T a,m v ' Y a,m dr 


d 2 cj> 2 

(o) r + 2 ^ — (o) jj- + . . . 

dr 


and 


dtj> R d 2 (p 2 

<P a „(r) = <P Q „(o) + ^ ,m (o) r + 1* 51 (1) ~+... 


rQ T o ' 

p,m • 3 ,m 


dr 


where the first terms are zero by virtue of the remark done before. 

We thus have 

. dcj) .. dcj) Q d(j> / dcj> d 2 cj) \ 

I a ,™ a _ 1 a = (o) + ( 0 )-4 

r dr %,m r dr 9 a,m dr \ dr K } dj .2 K) JT'j 


dcj) /dcj) d 2 cj) \ 

..li° (_m (o) + — a*a (o yi_+ ) 

dr \ dr w 2 '' o; 2! "'J * 


This last expression evaluated at this origin gives us zero. We 
finally obtain the predicted orthogonality relationship 
1 


/ r <p (r) cj) R (r) dr = 0 
J a,m 3 ,m 


'0 


Case B : m = 0 

This case corresponds to radially symmetric vibrations, and as 
a result the physics of the problem introduces as value for the first 
derivative at the origin: 
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when r = 0 


dr dr 


A development similar to the one done before for m#0, where the last 
remark has to be substituted gives us 


- e 


r <j> 


a,0 dr 


JL d4> B.O , 

r dr J " r 


d_ 
0 dr 


JL 

r dr 
1 


+ <X a,<f X 6,0 ) f 0 r *a,0*B.O dr ' 0 ' 


Let us consider the expression 

d 


r <p 


a,0 dr 


1^0 


r dr 


d 2 <p 


= <t> 


£» 0 


d(j> 


&lQ_ 


a,° dr 2 


r a,0 dr 


in the neighborhood of the origin. Considering that 4>„ _ is analytic 

d<Po o 

in the neighborhood of the origin, we can expand — in the follow- 
ing way 

d<P a n d<P Q n d\ n d\ n 2 


6.0 


dr 


dr 4 dr 


dr 

d4>, 


6.0 


Taking into account the remark done before about — (o) , we have 

d<l> c 


. _d 1 
r ^a,0 dr r 


r £>0 


dr 


d \ 0 

^a,0 dr 2 " < * >a »° 


d 2 (p 


3^o 


(o) 


dr 


d 3 (J) 


i*0 


dr' 


(°) Y\ + * • • 


Evaluated at the origin, the last expression becomes zero and we 

finally proved also that 
1 


/ r *a,0 (r) Vo W dr ‘ ° 
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Now that an orthogonality relationship has been established for the 

modes, the substitution of the asymptotic expansions valid for each 

til 

of them brings us the following relation for the zero order terms. 


1 211 

a 


r h~ (r)h_ 0 (r)cos m0 cos n0 drd0 = 0 . 
0 9 CL y IH (J y fj y XI 


Again, our main interest will be oriented towards the case where 
m=n. The case m^n is, as mentioned before, automatically satisfied 
by performing the integration in 0 first. The differential equation 
defining h^ ^ and h^ ^ are given by: 


Id f . 2 2. dh 0,a 1 m^ n , 2 2 , , , ,2 , n 

r dr [ r ( 1-k i r ) dr J f 1-k 2 r ] h 0,a + X 0,a h 0,a 0 

id 

* h - 


m /, ,2 2, , . , 2 . n 

^2 1-k 2 r ] h 0,B + X 0, 3^0,3 ° 


or for r^O 


d r 

dr 


r(l-k 2 r 2 ) 


dh n 1 2 

0,a _ tn_ 

dr “ 2 

J r 


»- k 2 r X.a 


+ X 


n r h~ 
0,a 0,a 


= 0 


and a similar expression for hg g. But these last expressions repre- 
sent particular cases of the more general eigenvalue problem of the 
Sturm-Liouville type 

4- fp(r) ■ 1 - g(r)h n + A. 2 r h A =0. 

dr ' dr J 6 0,a 0,a 0,a 

If we now consider the boundary conditions of our problem or 

h_ =0 when r = 1 , 

0,a 

combined with 
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or 


=0 when r = 0 for the case m 4 0 
0,a r 


dh 


0 ot 

, — 1 — = 0 when r = 0 for the case m = 0 , 
dr ’ 

our problem satisfies the requirements of the Sturm-Liouville problem. 

We can then deduce that for our problem, there exists a countable 

sequence of eigenvalues and a countable sequence of eigenfunctions 

satisfying our problem. In addition to that, we have that 

1 _ 

r h 0,a (r)h 0,6 dr " 0 


L 


and the orthogonality of the zero order is automatically satisfied. 

The orthogonality relationship we derived before has to be 
completed by the relationship normalizing the modes. We then define 
i 2Ji 


/ ( r <}>^ „(r) cos^m© d0 = 1 . 


j2 ( ' 2 

/ Q J Q a,m 


Integrating with respect to 0, we finally obtain 
1. 


f f <p 2 


m (r) dr ’ If ' 


Introducing in the last expression, the asymptotic expression used for 

cb , we have the following relations 
a,m ° 


/ r h o (r) dr - n 


J r h Q (r) ^(r) dr = 0 . 
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The first relationship will be used in order to define the con- 


stant bg, still to be determined in the expression of h^, and the 
second relationship will be used in order to determine the value of 


X^. We had 


so we have 


and 


h Q (r) = b Q r m F^a,3,l+m;- 

/ 2m+l ,22 / . _ . 

r b Q F ^a,g,l+m; 


3+V r 2 


l+2k 2 +v 


3+v 21 

r 


l+2k +v 


) dr = n 


b 2 = — x 

o n 


/ 2uH- 1„2/ „ ' 3+v 2\ , 

It F (a,B ,1+m; — r ) dr 

■Q \ l+2k 2 +v / 


In order to determine X^, let us consider the differential equation 
defining h^(r) (3.10) or 

d f ,, ,2 2, dh ll m 2 M , 2 2. , , ,2, ' 2 . 

— [rd-^r ) ^7 J - — [l-k 2 r ] h^ r X^ = -X^ h Q . 

Multiplying this last relation by h^, and integrating by parts the 
first term, we have: 


2 2 d ^l 

r(1 ‘ k l r ) “~d7 h 0 


C , 2 2. dh l dh 0 , /*m Z . 2 2., , 

0 " J Q r(1_k i r } 17 17 dr - J n T (1-k 2 r )h i h o dr 


dh, dh, 


■ X 
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2 

Considering the fact that h^A^ >1) = 0, and integrating by parts once 


more , we have 


1 1 


-jtU-kjr 2 ) ^ hj] o + jf 37 1?] dr 

- f' 4 a-^r 2 )hih 0 dr - - h- 


or 


but 


The use of the differential equation (3.6) defining h^, gives us 


2 ^0 

17 h i 


-fr 


1 0 


x l h 0 h ! dr - - ■¥ 


2 2 ^0 
■ x i-- n(i“ k I) 17 h i 


^(l) = - 


A-k? 


dh o 

and — (1) - - C ] _ 


so 


dr h l 


1 4-k? 


4? 


2 ( 3+v) dF / „ , . 3+v \ 

■ - ■ j- (a, 3 ,l+m; 5 — 

4.9V 2 ^, dx \ l+2k 2 +v 1 


l+2k +v 


Finally, we have 




t2 


2 (3+V) dF / „ . . 3+V \ 

— N j- a,3,l+m; j— 

-l+2k 2 +v d \ l+2k 2 +v / 


/ 2mfl 2 / Q 

r F ^a,3,l+m;- 


3+V 2 , , 

7 — r dr 


l+2k +V 


) 


With the determination of A^, the obtention of h^(r) becomes an easy 


task and a well specified problem. 
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In order to complete the derivation done for this membrane, we 

should look at the problem that could arise at the origin. We know 

that the origin r=0 represents a singularity for the linear operator 
2 

and, as a result, the fact of neglecting the contribution of this 

operator, close to the origin, might seem incorrect. When we looked 

4 

at the reduced problem, we overlooked the term -e <j> ^ ^ appearing in 

Equation (3.5). Close to the origin, this term could represent a main 

contribution. In Appendix III, we address ourselves to the problem 

near the origin and we find that the solution previously found for 

r#0, is also valid through the origin. This result might be expected 

in the sense that there are no boundary conditions imposed at the 

origin and the solution found for r^O represents an expression that, 

4 

when introduced into -eV d> represents only a bounded function at 

m a,m r J 

the origin. Also, the presence of the singularity at the origin is 
simply originated by the choice of polar coordinates. For a rectang- 
ular plate described by cartesian coordinates , the center of the plate 
doesn't contain any singularity. 

With this last remark, showing that the solution previously 
found is valid in the complete interior of the membrane, we have solved 
the proposed problem. To conclude this section, we summarize briefly 
some of the main results. The eigenfrequencies of the taut membrane 
with small flexural rigidity has been taken under the form 


2 1+2 k+V ^2 ,2 

a) = s Q A 

a,m 8 a,m 


l+2k +v ~2 ,, 2 


8 


(Aq + JZ A^+. ..) 


where 


£ = 


2 E h 


3(l-v 2 ) (l+2k 2 +v)a 4 yft 2 
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The value of Aq is obtained by solving 


F(a^ 1+m - a, 1+m; k^) = 0 


and 


1 / 2 2 2 0 
a(l+m-a) = — ( 2m + m -m k^ - — 

\ k l 

2 

The relation defining A^ is also repeated here 


A 2 = v l-k^ 
1 1 K 1 


0 , 2 dF , , . ^ ,2,1 2 

2 k^ — (a, 1+m-a, 1+m; k^) 

__ 

/ r 2 - 1 F 2 (a, 1+m-a, 1+m; k 2 r 2 )dr 


In these relations, the following definitions have been used 


, 2 3+v 

k l 2 

l+2k +-V 


and 


2 = l+3v 

3 3+v 


3.3 Dynamics of a Rotating. Free. Elastic. Circular Membrane 
We now end the chapter, 3, by considering another type of 
elementary structure. It represents the last analysis of modes under- 
taken in our research using continuous models. This time, the "object" 
under investigation is a rotating disk of small flexural rigidity, free 
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The basic physical assumptions made in this study are the same 
as those made earlier for the taut rotating membrane. This disk is 
assumed to rotate with constant angular velocity around its axis of 
symmetry and the equations of the motion are given in a reference 
frame located at the center of mass of the disk when at steady state 
deformation and rotating with the system at the constant angular 
velocity fi. In contrast with the previous example, the spin has a 
stabilizing effect on the structure and loosely speaking, the 
structural stiffness of the disk is increased by the sometimes called 
geometric stiffness, induced by spin. The origin of our angular vari- 
able is chosen such that no phase angle is needed in our equations. 

The generality of our problem is not affected by this assumption but 
the text is noticeably lightened, Even though the membrane will 
undergo a steady state radial displacement, we assume the physical 
characteristics of our system to remain constant throughout the whole 
disk. 

Other classical assumptions, such as the planar state of 
stresses, are also adopted. There will be no variation in the radial 
and tangential stresses over the thickness of the membrane due to the 
assumption of small thickness. This last assumption merits one more 
comment related to the boundary conditions valid at the free edge. We 
have, indeed, to apply as boundary conditions at the free edge the 
Kirchoff condition and the zero bending condition. Let us develop 
those considerations for a while. 
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For the free edge, no restrictions are imposed on the deflection; 
however, there should be no bending moments, no twisting moments and 
no shear forces. There is no reason to disregard any of those as 
secondary effects. Kirchoff suggested the following way out of this 
difficulty. 

Calling the bending moment M, the twisting moment H and the 
shear forces N, he considered the edge x = a of a plate as being acted 
upon by distributed torsional couples, as in Figure 11. Their 



Figure 11. Kirchoff Condition. 


distribution was assumed to be nonuniform, in general, and the magni^ 

/ . 

tude of the moment of these couples per unit of length of the face, at 
a given point, is denoted by H^. The magnitude of the moment per 
length dy of the face is obviously equal to H^dy. It is noted that 
from the viewpoint of statics, the distributed torsional couples are 
equivalent to a certain shearing force. Indeed, the couple with the 


moment H^dy may be applied by means of two equal and opposite forces 

H 1 acting at the edges of an area of length dy. The couple 

9 ^ 

(Hi + —7^ dy)dy on the adjacent area of length dy may also be applied 
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3H. 


as two opposite forces ^ dy with the arm dy. Having done this 

for all areas of the face concerned we see that the forces applied at 

points m and n on the boundary of two areas reduce to a single force 

3H, 


3y 


dy per length dy. We observe, however, that there will remain two 


A B 

nonvanishing finite concentrated forces and at the edge A and B 
of the face. 

Hence, we conclude that the distributed torsional couples of 

intensity H are statically equivalent to the distributed shearing 
± 3^ 

forces which has the intensity — and two concentrated shearing forces 
and at the corners of the plate; if the contour is smooth, and 
without corners, these concentrated forces are absent. Proceeding 
from this, Kirchoff proposed that the three boundary conditions at the 
free edge be combined into two by equating to zero the bending moment 
M and the shearing force N and by adding to the latter the term 3H/3y 
which reflects the influence of the twisting moment H. We now arrive 
at the following two conditions at the free edge: 

M = 0 

N + |2- = 0 . 

dy 

Written in polar coordinates, those previous boundary conditions 
become: 


3 

3p 

_3 

3p 


w (l 3w ' 1 3 2 w \ _ n 

2 p 2 38 2 ' 

f* 3 2 w 1 3w 1 3 2 w I 

1- So 2 P 9p n 2 3fl 2 j 


+ (1-v) 


P30 




o 


3p H p 36 

valid at the free edge of the disk, where V is the Poisson modulus. 
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As stated previously we know that the dynamic equation for the 
transverse vibrations of a circular plate, subjected to radial and 
angular stresses is given in general form by 


_E 

12(1 


h 2 n 4 1 3 / 3w\ °Q 

t— V w + — t— la p t— )+ ~ r 

i-v 2 ) p p 3p \ p 3p / p 2 


2 2 
3 w 3 w 

2 “ p 2 
36 3t 


(3.11) 


where and Og represent respectively the radial and angular stresses 
h stands for the thickness of the disk and 


7^ is the Laplacian operator which expressed in polar coordinates 


is 


2 2 
V 2 = -2I+I-3+- -± 11- 

p 3p 2 p 3p p 2 30 2 


y is the disk mass per unit of volume. 

For this case too, we found that the deformation equation repre- 
sents a perturbation with respect to the steady state deformation, and 
once more, the radial displacement has to be determined first — a 
similar procedure was adopted in the previous case. 

The equilibrium of the forces applied to an element of the disk 
along the radius is given by 


df (p 0 p > - o e + p aV - 0 

where there is no dependence in the 0 variable due to the symmetry of 
the problem. The use of the generalized Hooke's laws expressing the 
relationship between stress and strain, give us as differential equa- 
tion for the steady state radial displacement u: 


2 d^u 

5 a 2 
dp 


+ P 


du 

dp 


- u = - 


1-V 


n 2 3 

y ft p 
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The general solution of this equation is 


u = 


(l-v)Cp - (1+v) 


1— V n 2 3 
— ET VJ . n P 


where C and are arbitrary constants. 


For a disk without central hole, has to be zero, to keep the 
displacement bounded at the origin. The corresponding stress compo- 
nents are now found from: 

_ 3+v n 2 2 

a p = e - — y ft p 

„ l+3v _2 2 
a 0 = e - -g — y Q p . 


The constant C is determined from the condition existing at the 
periphery (p=a) of the disk. If there is no force applied there, we 
have 

(a ) =0 

p p=a 


or 


_ 3+v n 2 2 

C — g y a , 


from which we have 

_ _ 3+v 
P 8 


n 2, 2 2 ^ 

y ft (a - p ) 


2 

a Q = [(3+v) a 2 - (l+3v)p 2 ] 


Let us now introduce the dimensionless variable, p = ar, so that 
Equation (3.11) takes the form: 

Eh 2 „4 . 1 9 r . 9w, . a 0 9 2 w 9 2 w 

- 7 x— V w + ~x— -r— [r a x— ] + — r—r — x =y — x 

n 4, 1 2. 2 9r 1 r 9r 2 2 _ Q 2 ~ 2 

12a ( 1 — V ) a r a r 90 9t 
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where 


n 2 2 ,3+V Wl 2 N 
O r = p a (-g-) (1-r ) 

2 2 

°e = H T ~ t<3+v) - Cl+3v)r 2 ] . 


Dividing through by p Q , we obtain 


E h 


12yfi 2 a 4 (l-v‘) 


_ v ^w + ± -1 [ r (l-r 2 )^ ] 

2, 8 r 3r 1 ^ ; 3r J 


. 1 r ,_, s /-Ii-a \ 2, 9 2 w 1 3 2 W 

+ — T [(3+v)-(l+3v)r ] — - = — — - 
8r 90 z ST 9t 


This last equation is suitable for the use of separation of variables, 
providing that the given initial condition of the problem is also 
separable. We will then use an expression of the form 

i (w t + m8) 


w - (r)e 

a,m 


a,m 


where ^ is the eigenfrequency of the vibration. 

In this last expression, no phase angle has been introduced for 
a reason explained earlier and this doesn't affect the generality of 
our problem. 

The equation of main interest is given by the differential 
equation defining ^(r) • This last expression can be written 


E h 


4- [r(l- r 2 > 


12(l-V 2 )aW m a ’ m 8r dr 


dr 


2 0 co 

^2 [<3+v)-(l+3v)r 2 ]<J> =- 

8r 2 ' a ’ m « 


2 y a,m 


where V m is the linear differential operator defined by the following 
relation 
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2 2 
^2 _ d 1. d_ _ m_ 

m , 2 r dr 2 

dr r 


We now define the following quantities: 

2 


e = 


2 E h 


3(1-V 2 ) (3+v)aVft 2 


a,m 


k 2 = 


8 a) 


cum 


(3+v)Q 

l+3v 


3+v 


The final expression for the differential equation for <p^ ^(r) is. then 


given by 


- evS + i f - [r(l-r 2 ) ^=1 - 4 U- 
m a,m r dr [ dr J 2 1 


2 2 

k‘ r ]4> + 

a,m 


+ X <f> = 0. 

a,m a,m 


(3.12) 


Applying the previous separation of variables to the boundary con- 
ditions of our problem, we obtain the boundary conditions valid 


for $ (r) or 

r a,m v 


d 2 4> 


a,m 


d_ 

dr 


/ 1 ^a,m m , \ n 

\r ^ *<»,■) " 0 

- ™ M ■ 0 


dr 

d 2 <p . d<p 2 

a,m JL a.m _ m_ j. 

j 2 r dr 2 ™a,m 

dr r ’ 


when r = 1. 

Here, again, the small flexural rigidity of the disk is combined 
to the high spin rate to produce the small parameter £. The problem 
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is now formulated in a way where a procedure similar to the one used 

before could be applied. 

2 

The operator V , defined earlier, contains a singularity at the 
origin r=0 and in the use of a perturbation method, some attention is 
required when considering the limit process. However, the solution 
found in the central area remains valid in the neighborhood of the 
origin. There are, indeed, no boundary conditions at the origin and 
by a procedure similar to the one adopted in Appendix III, we can show 
that no particular behavior occurs near the origin. 

Expansion Valid in The Central Area 

When r f 0, the expansion used for m (r) take the form: 

<t> a m (r,e) V r > + v 1 ^ e > h 1 (r) + v 2 ^ e ^ h 2^ r ^ + " ’ ’ 

The subscripts a,m will be momentarily dropped, and used explicitly 
only when needed in order to differentiate between different modes. 
Introducing for the eigenvalue the expansion: 


A (e) = 
a,m x 


A" + 


K 1 (e)A 1 + 


k 2 (e)A^ 


+ • 


the basic differential equation for the mode (3.12) becomes: 


-e V (h + V (e)h +. . .) 
m 0 1.1 


’ f h [ r(1 - r2) "df] + V £> r dr [ r(1 - r2) “d^] 

2 2 
m r . , 2 2., , N m ri , 2 2,, 

■ — j [1-k r ]h Q - v^e) — [1-k r [1^ 

r r 

■^Aq + K 1 (e)A^+...^ ^h Q + v 1 (e)h 1 +. . .^ =0 . 
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The differential equation defining hp(r) is then given by: 

2 


1 

r dr 


,, 2 , dh 0 

r (i-r ) 


m , 2 2. ,2 _ 

— (1-k r )h Q + A Q h Q = 0 

r 


(3.13) 


or 


dh 


2 ^ ^0 1 2 
(1-r ) — + - (l-3r ) , 

,2 r dr 

dr 


0 


1 r 2 , 2. 2. , 2. 2,. 

-j [m -(m k + A Q )r ]h Q 


The point r=0 is a regular singular point, for this last equation, and 
adequate for a solution in power series around r=0. Instead of trans- 
forming the previous equation to obtain a solution expressed in terms 
of the hypergeometric series, it is indicated to expand the solution 
directly under the form: 



n+c 

a r 
n 


where c is the index of the previous equation and a^ an arbitrary 
constant. 

Introducing the last expansion of hg(r) into the differential 
equation, (3.13), the indicial equation gives us 



Only the positive root is feasible for the problem, for its boundedness 
of the displacement at the origin, so c=m. Also, the recurrence 
relationship between the a_^'s is easily obtained and we report the 
result 

a (n+m) (n+m-2)-(k^m^ + A!?) 
n _ . ; u 

a n-2 n^ + 2nm 
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By direct observation, we can see that the denominator is always 


different from zero and the solution hg(r) contains only odd or even 
powers of r. In order to pursue our study, it is interesting at this 
stage to look for a while at the boundary layer valid near the free 
edge of the disk. We need then to define a suitable boundary layer 
coordinate or 

~ 1-r 

r <Ke) ‘ 

2 

The operator becomes : 



2 -2 
<f> dr 


1 d 

l-<f>r. <j>dr 


2 


m 

(l-<j>r) 



So the basic differential equation for cj)^ ^(r) (3.12) becomes in the 
new boundary layer coordinate: 


’ d 2 

1 

d 

2 

m 


d 2 cf> 

Y a,m 

1 

dcj> m 2 (f) 

T a,m a ,m 

cj> 2 dr 2 

l-4>r 

4>dr 

(l-cf>r) 2 _ 


2 ~2 
<f> dr 

l-(pr 

(1-0,5) 2 


_a_ _d_ r (1 . 

l-(pr (j)dr |_ V 


2 

4>r)[l - (!-<*>?/ J 


cj> 

a,m 

<j>dr 


— r [l-k 2 (l-^) 2 j <f> + A 2 4> = 0 

,, 2 a,m a,m a,m 

(l-q>r) 


(3.14) 
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or 


-£ 


_i 

dr 4 


d 3 <j> 


<P 


<*+*> <j) 3 dr 3 


_ ( 2m 2 +i) — i 


d 2 <J) 


~ 2 2 ~2 
(l-cf>r) <|) dr 


2 1 
- (2m +1) 


deb d) 

q,m . , 4 , 2^ a,m 
7 ~ Tj~ — + ( m - 4m ) — 7" 

a -*;> 3 ^ o -* f ) 4 


<1-Ur) 2 ' ,1 


o o dd> 

(3<j> 2 ? 2 - 6<j>r + 2) 3 *- - 


dr 


„ „ . d 2 cb 

i /i2~3 «-> ^,^2 , ^ ~ \ ex t m 
+ (q> r -3 <pr +2r) — ■- 1 

dr 


m 


(!-<}>?) 


2 [l-k 2 (l-<|>?) 2 ]<J> m + X 2 4 = 0 

^ 06 ^ IQ 06 ,m 06 y m 


The thickness of the boundary layer is then obtained by balancing 




e ~ 1 

4 <j> 


i . e 1/3 


If we define now the asymptotic expansion of d> valid near the 

q,m 

free edge in the following way 


4> a> m^» £ ) = <S 0 (e) f 0 ^) + <S- L (e)f 1 (f) +..., 


then the differential equation defining fg(?) is given by 


d V E) d f o . A 

-J— + 2 it + 21 tt 


= 0. 


dr 


dr 


On the other hand, the boundary conditions become in the new boundary 
layer coordinate 


d 2 <J) 


q,m 


i 2 ,~2 

9 dr 


-v 


d<4> 


(!-<}>r) <J>dr 


q ,m ^ m 


t 


(l-(pr) 2 a,m 


= 0 


(3.15) 
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and 


d 3 <j> , d(J> 

g.m 1 g.m 


d<J) 


<j> 3 d? 3 (l-<pr) 2 4>dr 


g.m 


2m 


(l-«J>r) (j)dr ' (1 _ ( j )? )3 T a,m 


m 


d<}> 


g,m 


(l-(j)r) 


2 <j)dr 




(1-v) 


m 


(l-cj)r ) 


d(j> 


g.m 


(l-4>r) (f)df 


(l-4>r) 2 a,m 


= 0 


( 3 . 16 ) 


when r = 0. 

For the first term of the expansion they become: 


d 2 f 


dr 2 


° (r) = 0 


d 3 f 


when r = 0 


dr 3 


° (?) = 0 


The integration of the differential equation defining f^ and the use 
of the boundary conditions just derived limit the solution for fg(r) 
to a constant Cq. The determination of this last constant is 
determined through the matching process near the free edge. For 
matching, an intermediate limit suitable for matching near r = 1 is 
given by 

r = , r ; - where r is fixed 

n n(e) n 


and ri(e) 0 


which implies 


when £ -*■ 0. 

1/3 ~ 

£ r 

n(e) 


n/e 1/3 - - 


What we want is r 


00 when £ -*■ 0, or 
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The behavior near r=l takes the form: 



but h^Crir^) can be written 

h 0^ r n^ = h 0^ " dr h 0 nr ri + "* 


and fQ(r) = Cq. So we finally have: 


Ub h 0 (l)- ^ h o a) nr n +..; 
£-*■0 1 
r fixed ' 

n 


-v 1 (e)h 1 (nr )+. . .-6 0 (e)C 0 



The matching is then giving us 

6 0 (e) = e° 

and 

h 0 (1 > - c 0 • 

2 

The determinations of hg(r) and Aq are now possible. The con- 
dition h Q (l) = C Q — given through matching — appears, at first, to be 

unsufficient. But the consideration of the recurrence relationship 

2 

between the a^' s shows that a sequence of value of Aq truncates the 
expansion for hgCr) , namely: 

2 2 2 

(n+m) (n+m-2) -km = Aq . 

2 

Only the sequence of values we just generated for Aq truncates 
the series and gives a bounded displacement near the edge. The value 
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h^Cl) is indeed given by 7 , a . By comparing the ratio of two con 

n=0 n 

secutive terms of the last sequence to the ratio of two consecutive 
terms of the diverging sequence defined by 

i + JL_ + _L_ 

p p+2 p+4 • * * * 

it is easy to prove that the ratio generated for the a^'s is greater 
than the ratio defined for the diverging sequence. The boundary con- 
dition for hg(r) is in fact 

hg(l) = bounded 


It is concluded that for every value of different from those 
truncating the series, the displacement along the edge is unbounded. 

Having now completely defined the main contribution of the dis- 
placement valid in the central area, let us go back to the boundary 
layer expansion valid near the free edge. The consideration of the 

expression developed earlier for the matching suggests as an order of 

1/3 

magnitude for 6^(e) the value e 


Going back to the development of the differential equation 
defining the modes expressed in the boundary layer coordinate (3.14) 
we find that the differential equation defining f^ is given by 
df. 


4 2 

dV df d f. 


d 3 f 


dr 


dr 


0 ~ df 0 , ~~2 d f 0 
= -2 r + 6r — -p? + 3r r 

df 3 dr a ? 2 


, , 2 .22 , 2 . 

+ [m - k m -AqJ f q . 

Taking into account the expression already found for fg or 


f 0 (r)'-h 0 (l) . 
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/ 


We have 


dr 4 


2 

^1 ~ ^ ^1 2 2 2 2 
+ 2 -^4+ 2r — ± = [in - kV - A*]h (1) 

dr 1 


(3.18) 


The expansion of the boundary conditions (3.15) and (3.16) expressed 
in the boundary layer coordinates also give us 


d 2 f 


i (r) = 0 


dr 


A 

dr 3 


when ? = 0 


= 0 


The integration of the previous linear differential equation and the 
use of the boundary conditions defined for f^, gives us as a solution 

for f the following expression: 

,2,22.2. 

(m -km - A_) 

f 1 (r) = h o(l) ? + C ± . 


Introducing the last expression into the matching process , we have 
d 


lim 
e -+ 0 
r fixed 

n 


- I7 h o (1) nr n + ---+v 1 (e)h 1 (nr Ti ) 

r^2_,. 2 . 

.1/3 


- e 


, 2 , 2 2 ,2. 

(m -k m - A ) nr 

Q_ — n~ 

2 n o u ' 1/3 


+ c. 


+ . .) =0 


The matching process should determine the value of C^. 

directly that the expression 

2,22 ,2 

, m - k m - A n 

-d7 h o (1) 2 — h 0 (1) 


We observe 


is identically zero as seen by evaluating the differential equation 


defining hg(r) , Equation (3.13), at the point r=l or 
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2 


= 0. 



r 2 2 , 2 , 2 t . 

[m - m k - A q ] h Q 


The determination of the constant seems to depend on the value of 

Contrarily to our expectation, the value of v^(e) is not of the 
1/3 

order of magnitude of £ . The boundedness of the displacement along 

the rim of the disk and the orthogonality relationship — valid for 

1/3 

this case too, as shown later — cancel the choice e for <^(e) and 
v^(e) in their respective asymptotic expansion. As a result, the only 
feasible value for is 



In Appendix IV, it is briefly presented why the choice k^Ce) = 

1/3 

£ in the asymptotic expansion of the eigenvalue has to be rejected. 
We now prove the orthogonality relationship between the modes. This 
last relationship is used in the determination of the next term of the 
asymptotic expansion valid for the eigenvalue. The orthogonality 
relationship between the modes will once more be given by 



r (b d > 0 cos m0 cos m0 drd0 = 0 
CL ) m p j n 


It is obvious that, when m ^ n, the relationship previously written is 

satisfied, as a direct integration with respect to 0 will give. The 

case of major interest is in fact given for the case where m=n. For 

that case, we thus have to prove that 

.1 

/ r <J) <p„ dr = 0 . 

a,m p,m 
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The differential equation defining 4> a m is given by 


-E V 


m 


, , 1 d ,, 2 S a.m m ,..22,. . ,2 

<p + ~ T~ r(l-r ) — t j — - [1-k r ]d> + A d> =0 

a,m r dr dr 2 Y a,m a,m Y a,m 


with the following boundary conditions: 


d 2 <j> 


a.m 


dr 


+ - 4 + ) - o 

\r dr 2 Y a,m/ 


:a , 2 i 2 , /* \ 

— + i -Jwa - =T t . a-v, a. o 

[_ dr 2 r dr r 2 a » m J r dr \ r / 


d_'. d * 
dr 


when r=l. Similar expressions could be written for <b B . Taking the 

p,m 

differential equation defining ^ and multiplying both sides by 

r<j) R and integrating, we have 
p 9 m 

1 1 r dd> i 

-£ f r <t> 4> ft dr + f r(l-r 2 ) -- a> - L dr 

J Q m Y a,m Y B,m dr L dr J Y B,m 


-I 

0 

*1 


— [1— k 2 r 2 ] <f> <j> dr 

r a,m 


A r <J> dr = 0 

a,m a,m B,m 


Expanding the first term under the form 

/ 

1 1 


-£ 


/ r (f> <j> dr = -e / r <j> Q [V' - <j) ] dr 

m Y a,m Y 3»m B,m Ar /. m Y a,m 


0 

1 


d__ r „2 
dr' 


-e f 4> q "T^ [V 2 (f) ] dr 

Y 3,m dr m Y a,m 


+em 


L r ^B, 


V 4> dr . 
m m T a,m 
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An integration by parts and the use of the boundary conditions trans- 
form the previous expression into 


-e(l-v) m 


2 -A 

dr 


feus),* + e f r ~ [V 2 4> ] dr 

\ r / ^ Jq dr dr m r a,m 


+e m 


f*JL 

*'o r 


m [V 2 4> ] dr . 

m T a,m 


Integrating by parts once more, and considering the boundary conditions, 
we obtain for the previous expression 


-e ( 1 -v) 


2 ^a,m^3> 

d <J> 


dtp 




d(j) d(J > 0 2 

a.m . A iS + 5L_ a 

m dr dr r y a,m dr 

1 


<1> 


a,m 


r 3 ,m dr 


-e 

II J 0 


f r[V 2 (j) ][V>. ]dr 

J n m T a,m J m r B,m 


Introducing this last expression into the expression found previously, 

we observe a perfect symmetry between ^ m and <J>g m * As a result, a 

similar procedure applied to the differential equation defining (})g ^ 

would give an expression that would differ only by the last term. 

Substracting both expressions, we are left with 

.1 


(X 2 - A 2 ) f r <p <p 0 dr = 0 

a,m 3»m J o a,m T 3,m 


which contains our required expression. 

To complete the study of this problem, we will need the next 
term of the asymptotic expansion valid near the free edge. The match- 
ing process valid near the free edge suggests for 62 (e) the quantity 
2/3 
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Going back to the differential equation defining the inodes, 
expressed in the boundary layer coordinates. Equation (3.14), we find 
that the differential equation for f 2 is given by: 

df. 


4 2 3 2 

d% df d% d J f df d Z f n 

f + 2 —4 + 2 r f = -2 7 + 4r -Hr + 3r 7 

-4 dr ,~2 ,~3 dr ,~2 


dr 


dr 


dr' 


dr 


, r 2.22 . 2 . . _~2 df 0 , 0 2 ~ 

+ [m -k m -A^Jf. - 3r — rsr + 2 m r f_ . 

0 1 dr 0 


Taking into account the expressions already found for 


and 


f 0 « - h 0 (l) 


, 2 .12 ,2. 

(m -km - 

fi (~) _ - hpC!) r , 


the last expression becomes: 
df 


2 2 

d f„ df_ ^ f n o 9 9 9 

-£ +2-^+2? — = 2(m z - km - Ag) h Q (l)f 


dr 


dr 


+ (m - k m - Aj) 


2 2 ,2. 2 h 0 (1) ~ . • 2~ 


2 r + 2m r h Q (l) . 


The expansion of the boundary conditions (3.15) and (3.16) expressed 
in the boundary layer coordinate give as boundary conditions for f 2> 
the following expressions: 


dh 2 df 

— I - v - d ? - " v f o ’ 0 

dr 


d 3 f 


dr 3 


2 2 d ^0 
^ + (1-V)m — ^ = 0 


dr 


when r = 0. 
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when r = 0 . 


Those last expressions can also be written: 

f 

A 2 f 2 » 2 2 , 2 

d f„ m - k m - A_ „ 

—A „ v - u - h 0 (l) + m V h (1) 

dr ^ 

dr 3 


= 0 


The equation defining can also be written 


d f 2 _d 

,~4 " z d? 
dr 


^2 

dr 


= -2(m 2 - k 2 m 2 - A 2 )h Q (l)r 


2,2 2 , 2 S 2 h 0 (1) ~ . 2 


- (m-km - Aq) 


2 r + 2r m h Q (l) . 


A particular solution is given by considering 


where 


D o - 


f 2p = D 0 * 


\ (m 2 -k 2 m 2 -A 2 )h 0 (l) + ~ (m 2 -k 2 m 2 -A 2 ) h Q (l) + ^ h Q (l) 


The solution to the homogeneous part is obtained, after a first 
integration, and the use of the second boundary condition by solving 


d 3 f 


2h 


dr 3 


- 2 


[*%■]-• • 


Defining df„, /dr = f„, , we have 

Zn Z n 


2 * 

d f 2h . * 

TT.- 2r f 2h = °* 

dr 


£ I / O 

The change of variables r = 2 r transforms the last equation into 
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d 2 f 


2h 


* * 


dr 


*2 r f 2h °» 


which general solution is expressed in terms of the Airy functions as 


f* h = D 1 Ai(r*) + D 2 Bi(r*) 


The asymptotic expansions of the two Airy Integrals are given in the 

A 

literature and show that Bi(r ) is inadequate for matching because of 
its exponential growth: we are thus left with 

f 2h D 1 Al(r * 


or 


f 2 h = D x Ai(2 1/3 r). 


The general solution for f 2 is: 


- V 2 + »x l 


Ai(2 1/3 £)d£ + D 3 . 


The determination of is done by looking at the first boundary con- 
dition 


d 2 f 


or 


2 , 22,2 

„ m -k m - A 0 ? 

_ v h Q (l) + m V h Q (l) when r = 0 

dr : 

2,22 ,2 

1/3 m • m . 2 

2D Q + D l 2 1/J Ai'(o) = V 2 “ + m v h 0 (1) 


From this last expression, is determined. The determination of 
is then obtained through the matching process by ensuring that the 
quantity 

r 

D x J Ai(2 1/3 Od? + D 3 
0 
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represents a transcendentally small term away from the boundary. The 
introduction of the expression found for f ^ (?) into the limit process, 
considered in the matching, leaves us with the requirement: 

1 ~~2 h 0 (1) ~ D 0 = ° 
dr 


where 


D Q = i(m 2 -k 2 m 2 - A 2 )h 0 (l) + (m 2 -k 2 m 2 - A 2 )^!) + f- h Q (l) 


The differential equation defining h.Q, Equation (3*13), expressed in 
the r ' variable becomes : 


2 d ^0 i 2 ^^0 

(2r'-r ' ) ~ + 37-,(-2+6r'-3r ,Z ) 

dr 


1 r 2 \ 2 \ 2 
r [m -m k - A 

I \ ^ v 


(1-r ' ) 


+ 2(m\ 2 + A 2 )r’ - (m 2 k 2 + A 2 )r’ 2 ]h n = 0. (3.19) 


O' 


0 


2 2 ■> 

Expanding d h Q /dr , dh Q /dr and h Q in Taylor series in the neighborhood 

of r 1 = 0, we obtain the following identities: 

d “n 2 2 2 2 

- 2 (1) - (m -mV - A Z ) h Q (l) = 0 

2 

d h» 2 2 2 2 dh 0 2 

4 (1) + (mVk Z " X q + *> (1) - 2m Z h Q (l) = 0. 

dr 


This last expression is also: 


,2, 2 V x 1 

d h. m -m k - A_ 

(i) _ : y. 

• 2 2 
dr 


, 2 2 , 2 

(m -m k 


h Q^ + 8 


2 2 

A q ) dh 


0 
dr 


( 1 ) 


+ T h 0 (1) 


2 2 

or d h^/dr (1) = 2 Dq which represents the predicted identity. 
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Having previously concluded the orthogonality relationship 
between the modes, we have to complete the specification of the mode 
by writing down their normality or 



A direct integration with respect to 0 gives us 
1 


/ 


.2 , 1 

r (p dr = ~ . 

a,m II 


This relation represents the relation needed for the determination of 
ag. Introducing the asymptotic expansion of the previous equal- 

ity becomes indeed 
1 


r a,m’ 


f 0 r h o dr * ff 


or by writing the development of hg(r) under the form 


h 0 (r) a 0 n=0 “n 


£ a' r 


n+m 


= a Q h'(r), 


we obtain finally 


a„ = 


j. 


(r)dr 


With the determination of hg(r) completed, we are now able to 

2 

proceed and determine the value of The orthogonality relationship 

between the different terms of the asymptotic expansion of the modes 
and the considerations introduced previously suggest for v^(e) the 
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order of magnitude of e. A suitable choice for K^(e) is also given by 
< 1 (e) = e. 

The differential equation defining h^ is then given by: 

■1 if,, 2. dh l | m 2 M .2 2,. . ,2 , _ n 4 ,2 , 

r dr [ r ^ 1_r ^ drj" _1 ^ k r -* h l + X 0 h l V m h 0 X 1 h 0 ' ( 3 - 2 °) 


Multiplying this equation by r hg(r) and integrating from 0 to 1, we 
obtain the following expression: 

/ h 0 (r) df [ r<1 - r2) I?] dr - L ' 7 U-k 2 r 2 )h 1 h 0 dr + X 2 f r h^dr 


■jf 


4 1 

r V h_h_ dr - 
m 0 0 IT 


where the normality relationship of h^ has been used. Integrating by 


parts the first term, we have: 

1 1 


2 dh l 
h Q (r)r(l-r ) 


0 ■'o 


f 2, dh l dh 0 . 
J r(l-r > "57 -*r dr 


f ~ (l-k 2 r 2 )h,h„ dr 


TO 


-L 1 ^ 

+ X 2 f r h h dr = x / r h h dr - — h . 
0 J n 0 1 Jq m 0 0 JL 


Integrating by parts once more, and taking into account the differen- 
tial equation defining h Q , Equation (3.13), we have: 


2 dh 0 

-rU-r ) h 1 


1 1 
■ + 


J 0 \ |/7 tt- k2r2 > h 0 - r A o h o. 


dr 


1 

' f Q *7 < 1 -k 2 r 2 )h 1 h 0 dr h^dr 

y 4 

= / r V h n h n dr - — =• 

J n m 0 0 n 
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or 


1 

A, '=* II f r ^ h L dr 

1 Jq m 0 0 

2 

The previous determination of Sq and this last relationship com- 

2 

pletely define the value of A^. In this problem, the flexural rigidity 

of the disk modifies the value of the eigenfrequencies only through 

an order of magnitude £. In a procedure similar to the one used for 

2 

the rotating cantilevered beam problem, the determination of A^ has 
been accomplished without the need of the orthogonality relationship. 

Here again, the general expression valid for h^ is given by an 
expression of the form 


h l = h l + k i h o 


where h^ is a particular solution which displacement at r=0 is zero 
and tig is the general solution of the homogeneous part. The determin- 
ation of the constant is accomplished by considering the ortho- 
normality ’condition. However, we have to be careful when doing this 
step. We already found that the boundary layer expansion contains a 

transcendentally small term in which order of magnitude is given x 

2/3 • 1/3 

by £ . The boundary layer being of thickness £ , we might conclude 

that the contribution of the TST in the integral appearing in the nor- 
mality relationship gives us a contribution of order of magnitude £. As 
a result, the uniformly valid expansion has to be considered in the 
normality relationship. 
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The uniformly valid expansion for 4> a m is given by 


\,m = h 0 (r) + £ h i (r)+ * • * + z 2 '\ 


where has been written 


A. 

f Ai(2 1/3 Od? + 


+ . . • 


l 


Ai(2 1/3 Od5 + 


represents a transcendentally small quantity away from the boundary. 

The normality condition gives us then: 

1 L. 1 

/ r (j)^ dr = = / r h 2 dr + 2e It hJi-dr +. . . 

Jo a ’ m ■ 11 Jo 0 Jo 01 


+ 2e 


2/3 




Ai(2 1/3 ?)d? + D' 


dr +. . . 


Let us concentrate on the last term: 


V 


/ r h o 

/ Ai(2 1/2 Qd? + D' 

/o 

Jo 


dr . 

The mixed appearance of r and r has been kept on purpose as the rest 
of the development will show. We replace r by 1-r' in the previous 
expression or 

1 


1 r r. 

\ y (l-r')h 0 (l-r') j Ai(2 1/3 Od£ + D' 


dr 1 


t 1/3 

The expression Ai(2 represents a transcendentally small 

quantity away from the boundary, when considered in the r 1 variable and 
this remark suggests an expansion of hgCl-r’) in Taylor series is the 
neighborhood of r' = 0. 
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1 

/• 

r 

/• 





fo h ° (1> 

Jo 

Ai(2 1/3 C)d? + D* 

dr* - jf r-h 0 (l) 

/ Ai (2 

Jo 


-I 



[i 

1 


D 3 

dr ' - d i y 0 fe h o <i)r ' 

/ Ai(2 1/3 ?)dC + DI 

Jo J 


1/3, 


dr '+. . 


1/3 ~ 1/3 

Replacing now r' by £ r and dr ' by £ dr, the previous expression 
becomes r 


e 1/3 D, 


i / 1/3 

r 


f 

h o (1) / 

/ Ai(2 1/3 C)d^ + D' 

•A) 

Jo 


- r • [ 


dr 


Ai(2 1/3 £)d£ + D’ 


“ 1/3 % 


Ai(2 1/3 0dC + 


dr 


dr +. . . 


The limit process applied to the normality conditions gives us then 
1 


f rh () h 1 dr + D 1 h Q (l) f 
V 0 J 0 


/“ 

,/ 0 


Ai(2 1/3 ^)d? + 


dr = 0 . 


Defining the integral of the transcendentally small quantity appearing 
in the last bracket as being D^, we have 


J f r h Q h ± dr + D 1 D 4 h 0 (l) = 0 . 


This last expression can now be used to determine the value of k^. 

f r Vl + k ], f r h 0 dr + D 1 D 4 h 0 (1) = ° 

•/ 0 •/ 0 
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or 


k l = - n 


f r h 0 h 1 + Dl D 4 h 0 (l) 


With the determination of k^, the expression for h^ is completely 
specified and we will then conclude here the study of the modal 
analysis of the last structure. Several properties that we mentioned 
for the previous problem can also be observed in this case, but for 
reason of shortening will not be developed here. 

We conclude this chapter by recording briefly some results in 
the determination of the eigenfrequencies of the free rotating, circu 
lar membrane. The expansion found for the eigenfrequencies has the 
form 


2 3+v „2 ,2 

0) = — r— A 

a,m 8 a,m 


3+v n 2,,2 , .2 ■ . 

g ^ (Aq + A^ e +. . •) 


where 


e = 


2 E h 


3(l-v 2 )(3+v)a 4 y ft 2 


The value of A^ is given by the sequence defined by 
(n+m) (n+m-2) - k 2 m 2 = A 2 


where 


and 


. 2 -l+3v 

k = 


3+v ’ 

1 


f 


r V h-.h_ dr. 
m 0 0 
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The obtention of results from 

00 


h 0 < r > 



where the truncated series is defined recursively by 

a n (n+m) (n+m- 2) - (k^m^ + Aq) 

_ = — 

n-2 n + 2mn 
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CHAPTER 4 


LIAPUNOV STABILITY ANALYSIS 

4.1 Introduction 

Due to the particular orientation of this chapter, we should, at 
the outset, investigate the concept of stability, and state precisely 
what we mean by a stable motion for our problem. Ue will then in a 
first part review briefly the different ideas pertaining to the stabil- 
ity and cover to some extent a peculiarity appearing in the class of 
problem under investigation. In this first part, the concept of stabil- 
ity will be developed and the definition of nominal motion will become 
more specific. 'In the second part, the Liapunov stability analysis will 
be studied for flexible space craft when continuous models are used, and 
finally analytic criteria for the class of problems studied in chapters 
2 and 3 will be obtained. 

4.2 Concept of Stability 

We might introduce this section on concept of stability by 
defining formally the notion of Liapunov stability. 

If, for a given dynamical system, the equation of state 
X = F (X, t) 

admits the solution X = 0, and if for any £ > 0, there exists some 6 > 0 
such that if the norm |x(0) | £ 6 then |x(t) | < £ for all t > 0, then the 
motion characterized by X = 0 is said to be Liapunov stable; and if in 
addition 
lim 

|x(t)| - 0 
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then the solution X = 0 is asymptotically stable. If X = 0 is not 
Liapunov stable, we will call it unstable. 

One remark should be made about this definition: only local 

stability properties are examined. In order to explore the kind of 
instabilities occurring in the system, nonlinear simulations become 
necessary. The reason is that the characterization of stability is 
based on local properties of the motion, permitting no formal conclusions 
beyond the neighborhood of the nominal motion. 

If we, now, consider the consequences of a small perturbation 
applied to a rigid body, initially rotating, about a principal axis 
(major or minor) colinear with its inertial angular momentum vector, the 
subsequent motion is characterized by a rotation of its spin axis about 
the inertial direction defined by the new angular momentum vector — the 
sum of the initial angular momentum vector and its perturbation. The 
space craft is said to "wobble" and the angle between the spin axis and 
the new angular momentum vector is called the nutation angle or "wobble 
angle". Because the wobble angle resulting from a sufficiently small 
perturbation remains within any preassigned value, the motion is said to 
be Liapunov stable. 

On the other hand, heuristic energy sink methods — based on the 
assumed presence of a hypothetical non-moveable energy dissipator, and 
applied to quasi rigid bodies — indicate a decrease in the wobble angle 
if the spin axis is the axis of maximum moment of inertia; we would like 
to define our state variables in such a way that we can call this motion 
asymptotically stable. Conversely, energy sink analysis shows that when 
the spin axis is the axis of minimum moment of inertia, the wobble angle 
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grows and the motion becomes unstable. These last observations are 
generally referred as "the greatest moment of inertia rule". The fact 
that in the perturbed motion of a stable quasi rigid body, the spin axis 
converges to the new angular momentum vector brings out an interesting 
point: the system is not "completely" damped when the motion is 

described in terms of coordinates measuring the deviation from the iner- 
tial orientation existing prior to perturbation. We will focus our 

attention on this last observation, by referring to a development by 

13 

Barbera and Likins. 

In the application of asymptotic stability to our particular case, 
special attention is required in the choice of elements in the state 
variable X. If the state variables include, for instance, the devia- 
tion of the component of the angular velocity along the body spin axis, 
then the null solution will not be asymptotically stable for any vehicle 
configuration, since after perturbation, the vehicle will not return 
asymptotically to its nominal spin rate. For the same reason, we cannot 
include in X any deformations which vanish only in the unperturbed 
state, but we must define deformations relative to whatever deformed 
state corresponds to steady rotation about the body spin axis at the 
rate appropriate for the actual angular momentum, whether before or 
after, perturbation. Similarly, we cannot use in the state variables any 
set of inertial attitude angles which are zero only when the body spin 
axis is colinear with the angular momentum vector prior to perturbation; 
the perturbed angular momentum will have in general a different orienta- 
tion, and we wish to characterize the motion as asymptotically stable if 
the angular velocity vector of the central body approaches alignment 
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with the body spin axis and with the perturbed angular momentum vector. 
We must therefore choose attitude angles for the central body which 
vanish whenever the body spin axis is aligned with the actual angular 
momentum vector, whether before or after perturbation. 

Our concept of stability is now well defined and the constraints 
on our definition of the state variable made more precise. We are in a 
position where two theorems, usually referred to as being elements of 
Liapunov's direct or second method, can be mentioned. 

Theorem 1 : The null solution X(t) = 0 of the differential equation 

O 

X = F (X) is asymptotically stable if there exists a function L(X) in a 
region around the origin both positive definite and strictly decreasing 
for all solutions in that region except for X = 0; 

Theorem 2 : The null solution X(t) = 0 of the differential equation 

O 

X = F (X) is unstable if there exists a function L(X) in a region around 
the origin both negative semi-definite (or sign variable) and strictly 
decreasing for all solutions in that region except for X = 0. 

Although these theorems constitute usefuls tools for the 
determination of the necessary and sufficient conditions for asymptotic 
stability, they do not present procedures for the generation of a 
testing function. However, in a wide class of applications, the 
Hamiltonian is a convenient choice. More specifically, for "completely 
damped", systems, the Hamiltonian is a suitable testing function when the 
total energy of the system is free of explicit time dependence. For our 
purpose, the concept of "complete damping" implies a dissipation of 
energy for any possible motion, in the neighborhood of the origin of the 
coordinate space adopted, except for the nominal motion. However, as 
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mentioned in the preceding section, the damping of a freely spinning 
body with internal energy dissipation is not complete in terms of atti- 
tude angles measured with respect to the spin prior to perturbation, 
since the vehicle cannot return to its original state. As a result, for 
such systems, the Hamiltonian is not strictly decreasing in the neigh- 
borhood of the null solution and asymptotic stability cannot be pro- 

14 

claimed. In 1969, Pringle provided a method to circumvent this prob- 
lem. The procedure consists to constrain the attitude angles through 
the angular momentum integral so that they measure the deviation of the 
body spin axis from the instantaneous angular momentum vector, before or 
after perturbation from its nominal inertial orientation. The attitude 
angles are thus defined in such a way that complete damping is assured 
and asymptotic stability can be predicted. The requirement of complete 
damping also implies a judicious choice of deformation variables; they 
must vanish whenever the vehicle adopts a steady state configuration and 
rotates about the body spin axis (at whatever rate is appropriate with 
the actual angular momentum) . 

4.3 Stability Analysis 

We specified in the previous section our objective and the path to 
follow in our stability study of flexible structures. In the following 
derivation, the analysis is kept as general as possible. We will 
restrict the class of problems covered by the study when needed and 
finally, specialize in a class of systems, including the flexible struc- 
tures studied before, when analytic criteria will be developed. 

We assume, at the outset, that the freely spinning flexible 
spacecraft under investigation is composed of two parts: a rigid part 
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and a flexible part. Later on, the flexible part will be assumed to lie 
in a plane perpendicular to the spin axis and passing through the center 
of mass of the whole structure when at steady-state deformation. 

The kinetic energy of the system, denoted by T, is derived from 


the general expression: 
2T = | p*jj dm 


- / i- 

•'AjB 


where £ is a inertial generic position vector and the capital letters 
A and B indicate that the integration is carried over the flexible 



In our force-free case, the system mass center is inertially 
fixed, and £ can be written as the sum c_ + p where £ is the. position 
vector directed from the system center of mass CM to a point N fixed in 
B, and p is a generic position vector directed from N. The point N is 
chosen to coincide with the CM when the structure is steadily spinning 
at the rate consistent for the actual angular momentum and hence elasti- 

J 

cally deformed by forces induced by spin. Thus 2T becomes 
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2T “/<£ + £)•(£+£) dm 

*/A,B 

= . c_ • J ( £ + £) dm + / p_ • (c + p) dm 

A,B yA.B 

The first term vanishes by definition of mass center, i.e., 

* (c' + £) dm = Cl for every t, 

JA, B 


so 


-l i 


= _c / £ dm + / (£ • £) 

/a.b */a.b 


2T = ff £ • (A + £) dm = _c I £ dm + I (p • p) dm . 
'A,B A, B A,b 

By definition of the center of mass, we also have: 


L 


c + | £ dm = £ for every t, 

,B 


where represents the total mass of the system. 
The kinetic energy takes, then, the form: 


2T = - (£ • £) + / (p • p) dm + / (p • P) dm 

JA~ ~ J B 


/ ~ 
J B 


or 


2T = - (c • c) Ji + 00 • 1^ 

— — D 


CO +■ 


/*(£ • £ 

JA 


p) dm 


(4.1) 


where u is the inertial angular velocity vector of vector basis {_b} 

N 

fixed in B and is the inertia dyadic of body B about point N. 

Expressing the right hand side of equation (1) in matrix notation 
with respect to the vector basis {£} fixed in body B, we obtain for the 
first term 
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i/M (c • c) = ,A( [_c + (to x O • ^ + (co x j:) ] 

= [c^ • &_ + 2 c_ • (a) x c) + (a) x c) • (a) x c) ] 

T t ~ T ~T ~ 

= -M [c • c - 2d c u + to c cto] 

where («) represents time differentiation with respect to the vector 
basis Os} fixed in B and the tilde (~ ) operation is defined as 



The last two definitions can be extended to any vector used in the 
sequel . 

The second term becomes simply 

t N T t N 

(0 . I < (i) = u I u, 

— B — -b 

where 1^ is the inertia matrix in vector basis {b}, corresponding to the 
B 

dyadic 1^. 

To study the third term we express p as the summation of JT + _u, 
where _T is the position vector from N fixed in B, to the element of mass 
dm in the "steady state" position. The "deformation" ju measures the 
displacement relative to B of the element of mass dm from the location 
that it would occupy if the vehicle was rotating about the body spin 
axis at a steady state configuration at the rate appropriate for the 
angular momentum. With these definitions, we have successively: 
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X £ - £ 


dm = 


+ to x (r + u) . u + to x (F + u) ] dm 


+ _u • r + u 

. 

/. i 

J k 

/ (_u • u^) dm + 2 / u • [ oj x ( 1^ + _u) ] 

A 


dm 


+ / [to *(£ + _ u) • to x (T + u) ] dm 

A “ ~ 


/ u T u dm +2 / u^to(r+u)dm + / [ (r+u) to] ^( r+u) to dm 

/*u T u dm - 2 /d 1 

A A 


(T+u) dm + to x I I (F+u) T (r+u)dmlto 


Lf 


The Hamiltonian can be written as 


H - ~ to T 1^ to + - y"u T u dm + W T [f T T u dm + u T ? dmj ( 
+ | to T j^y*u T G dmj to - j^y* T (r+u) dmj to 


[c T c - 2c T c to + to T c T c to] + V + C Q 


(4.2) 


where the system potential energy has been defined as V plus a constant. 
We introduce the definition 


ij? = + / r " r dm 

u u 


ft 


i 


where the term j T T dm represents the inertia of the undeformed 

N 

appendage. The expression Iq defines precisely the inertia of the com- 
plete system about point N when undeformed. If the vector basis {bj} are 

N 

assumed parallel to the principal axis, 1^ is diagonal, i.e. 
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( A 0 0 

0 B 0 

0 0 c 

We are now in a position where we can compute the angular momentum 
integral. By definition, we have: 


h = # ( c^ + p) x ( £ + p) dm 

'A,B 


L 

L 


(c + p) x (c + p) dm + / (c + p) x [w x (c + p) ] dm 


= c x 


/< 

-/A,B 


(c + p) dm - c x I p dm + 


/ 

•/A,B 

/ £ dm + / 

.B A, 


p x p dm 


/ 

A.: 


+ c x |(i) x / (£ + p) dm + 

/a,b “ J 


/ £ * (“ 

A,b 


a,b 

x p) dm 


/ 

A.] 


- (to x c) X I £ dm . 

'A,B 

The first and fourth terms vanish by mass center definition. Defining 

N 

the inertia dyadic of the complete system about point N, I , we have 


h = I N • to + J( 


c x £ + #£ x 

Jk 


p dm 


The last two terms may be written in the vector basis {b} as 


•M C X c = [c + (U) X c) X c] = {b} ^Ai. (c c + c c to) 


/ £ x £ dm = / ( £ + £) x u. dm = {_b / ( 

Jk Jk Jk 


(f u + u u) dm 


which allows the matrix representation of h in {b}; 


t 


^h = 1^ to + / (f + u) u dm + (h. c + c c to) . 
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In order to constrain the Hamiltonian through the angular momentum 
integral, we now have to solve the angular velocity components in terms 
of the attitude angles. 

Let us define an inertial vector basis {n} and its transformation 
with respect to {b.} as 0, i.e. 

{n} = 0 {b} 

such that n^ is colinear with the angular momentum vector h after per- 
turbation from its nominal inertial orientation. As the system is 
assumed torque free ti must remain colinear with n^; thus, subsequent to 
perturbation the following must prevail: 

h • n = 0 , 

h • n 2 = 0 

_h • = h 

where h is the magnitude of the perturbed angular momentum vector. With 
h written in vector basis {_b} as prescribed before, the above equations 
become: 


0 




+ u) 


u dm + (c c + c c co) 



or 


[i N + , 4 ( c cl 

CD 

H 

O O 

- / (f + u) u dm - *y(( c c 

L J 

.Nh' 

Jk 


Let us, now, determine the transformation 0. We have 
{b} = 0 T {n} 

We will use in the three rotations needed to define completely 0, the 

T 

sequence of rotations 1-2-3. As a result the matrix 0 may be 
21 

written: 
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0 T = E- 8+ y08+-|eu0 

where E is the 3*3 identity matrix 

8 is a skew-symmetric matrix defined as 



0 is the diagonal matrix defined as 
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from its quadratic approximation, as long as all coordinates are 
included. 

The solution of (3) for the angular velocity components is then 
substituted into the Hamiltonian and limited to terms of order no higher 
than second. 

The inertia matrix is defined for the flexible appendage with the 
{b} vector basis as 




T T 

(P P E - p p ) dm 


where E is the 3*3 unit matrix. Replacing p by its value given in 
matrix form as 

p = T + u 


we have 


.N 


k 


(r T + u T ) <r + u) e - <x + u ) cr T + u T ) dm 


/ (r T te - rr T ) dm + /(u'rE - u r> dm 
Jk A 

J' (T T u E - T u T ) dm + y*(u T u E - u u T ) 


dm 


N 

The inertia matrix I of the complete matrix, computed with respect to 
the vector basis {is}, consists of terms independent of deformation as 
well as terms both linear and second order in the deformations. 


i N = + /(r T r e - rr T ) dm + /(u x rE - u r T ) dm + / (r T uE - r u T )dm 

B 4 Jk Jk 


k 

f u T U 

Jb. 


/<■ 


k 


E - u u ) dm 
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N 

The sum of the first two terms has been defined earlier as 1^, 
prescribed to be principal and having as diagonal elements A,B and C. 
We also define 


T T 

f E - u T ) dm + 


/<-■ 

A 

f (u T u 

A 


fr 1 u 

•A 


E - T u T ) dm = 2A 


T A 

E-uu) dm = U 


where we notice that A represents a symmetric matrix, so that 
N N 

I' - Iq + 2 A + U. 


The definition of the center of mass or 
,M c_ + / £ dm = 0 

Jk, B 

may be written as 


tAl + I £_ dm + / £ 

Jk A 


dm = 0 




// c + 


At dm + /u dm + / £ dm = 0 . 

A A A 


By definition, the point N has been chosen as the center of mass of the 
complete system when the structure is undeformed. This allows us to 
consider 


/£ dm + / £ 

A A 


r dm + J p dm = 0 . 

'b~ 

The vector c^ is then expressed with respect to the deformation u as 




dm. 


so that 
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= -J / 5 


dm 


r N 


The matrix [I J r l M c c] expands to: 




2A + u + i[jf * ] [/“ d 4 


The inversion of the previous matrix becomes easy if we consider the 

N 

fact that the matrix I Q is diagonal, and if we recognize the fact that 
only terms up to the second order are needed for the purpose of our 
derivation. We obtain the following expressions: 



Taking now, into account the fact that the matrix A contains only first 
order terms in the assumed small variables, and similarly, the matrix U 
contains only second order terms in the same variables, we have: 



+ iM cc 


-1 




AI 


N-l 

0 



UI 


N-l 

0 


-i /o ' 1 [/HE/ 5 '“K 1 




Introducing this last expansion into the expression previously found for 
the components of the angular velocity vector, and performing the indi- 
cated multiplication, we obtain as final expression, where only terms up 
to the second order have been kept 
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dm - Iq kM c c. 


When the expressions derived before for the angular velocity components 
are substituted into the Hamiltonian (4.2), a series of simplifications 
and combination of terms reduces the previous expression to: 


r N 1 „ t N _1 a 1 „ 1 r 


N _1 . 1 A 1 


H c- T< 00h ) *6 - 2I o AI o - 2I 0 0 + 4 x o AI o AI o 

-i l-i -i r /* IT r 1 -1 -1 -1 ~ 

N N IN l~ #~ N N N 

- x 0 u r o -jt T o r dm J u dm x o + 4 r o A x o 0 

LA JL A J 

— ■' ■ ~ ~ / 0 \ 

+ l“ [66 + 8U6] -6 ijj 0 ( 0 ) 

U \h/ 

+ Y \f u £ r dm Ig fv u dm + y J u T u dm 

[/A l^A j A 

- -bu [f k ° T<to ] l 1 dm j + v + c o 


(4.4) 


where H stands for the constrained Hamiltonian, 
c 

In this last expression, the potential energy must include the 
contribution of the steady state deflection of the flexible appendages 
induced by spin at a rate concordant with the angular momentum. In the 
following text, this spin rate will be identified as being This 

last observation expresses the fact that the expansion of the potential 
energy of deformation, in terms of the deformation variables, is not 
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accomplished with respect to the minimum energy of deformation state. As 
a result, the expansion consists not only of terms quadratic in the 
deformations, but also of terms linear in the deformation variables, 
plus a constant term. We limit the approximation of the potential 
energy of deformation to those terms because they represent a good 
indicator for the local behavior of the system. 

U d = L 0 + f k h \ ■ udm+ 2 

where represents the energy of deformation, and and K are made 
of functions of position. 

If we, now, consider our definition of u, the deformation with 
respect to the steady-state deformation induced by spin, we have to con- 
clude that the equilibrium position is defined when the deformation 
variables are zero. But for a conservative system, the equilibrium 
position is obtained when the dynamic potential, composed of the ener- 
gies of position and deformation and also of T^- the part of the kinetic 
energy not a function of the time derivatives of the state variables of 
the system — is minimum with respect to the state variables, when eval- 
uated at equilibrium. Only the deformation energy and the kinetic 
energy are functions of the deformation variables, and equilibrium has 
been defined when u and 0 are zero and when the spin rate is given by 
The kinetic energy has been determined earlier and the minimization of 
Tq - with respect to the deformation variables, evaluated at equilib- 
rium, is equivalent to the requirement that the linear terms in the 
deformation variables u appearing in accommodate themselves when 

the spin rate is Wgg or 

u dm + u T dm w - 
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The potential energy of deformation may then be written, after the 
observation that the terms in brackets are nothing else than 2A, in the 
following way 


U d * L o + “ss 4 “ss + 


\ /“ tk u 


.N 


dm . 
-1 /O' 


If we recognize the fact that togg is 1^ ^0 J and notice that 

-1 /o\ 


(0 0 h) !q 6 ^0 1= 0 , 


we obtain the following expression for H^: 


H c = ~ (0 0 h) 


N - "*" N - '*' N ^ N - ^ N ^ 1 n'^T fi 

“o 4I o 4I o - *0 ul o *0 [// d ” 

r — v, _"j ~ n - ^- 

e + i“ ee + 0 u 0 - 0 ij e 


1-1 -1 

N N N ~ 

dm I 0 + 4 I q il 0 ~ 


0 

0 

h/ 


\k 

+ i[^u T f dm] ig [jff u dm] + \ fi T u dm - ijjf dmj \Ji dm] 

+ ~ J ' u T Kudm . (4.5) 


This last expression contains only terms second order in the deformations 
and attitude angles. In the expansion of equation (4), several constant 
terms appear in addition to the arbitrary constant C q previously identi- 
fied. Since C q is arbitrary, it has been chosen in such a way that the 
summation of the constant contributions is zero, allowing the 
Hamiltonian to vanish at the origin — a necessary requirement for a 
Liapunov function. .. 

As H c in the presence of complete damping is strictly decreasing, 
then by Theorem I the nominal motion is asymptotically stable for 
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positive definite, and by Theorem II this motion is unstable for H 

c 

either negative semi-definite or sign variable. We assume that energy 
is dissipated within the flexible appendage whenever time-varying defor- 
mations occur, as in the case for any physical object. We have care- 
fully chosen our variables so that energy is dissipated for any admis- 
sible solution of the equations of the motion in the neighborhood of the 
nominal motion, except for the nominal motion itself, so we are assured 
of complete damping. 

Since for asymptotic stability the complete function must be 
positive definite, then it is clear that the following must be satisfied 
for stability 



2 


C-A 

AC 


. 1 h 2 fl 2 
+ 2 h ®1 


C-B 

BC 


> 0 . 


This result provides the familiar necessary stability criterion pre- 
dicted by energy sink methods for spinning bodies having an internal 
energy dissipator, i, e. 

C > A and C > B. 

Thus, by inspection of the Hamiltonian we can formally conclude that, in 
the presence of damping, the spin axis must be the axis of maximum 
moment of Inertia. However, our objective is to extract additional 
stability criteria, and this requires the determination of conditions 
for positive definite H^. 

Now, that the constrained Hamiltonian has been found in a 
completely general way, we are forced to consider more specific problems 
if we want to come up with analytic criteria. In order to introduce 
those particular case, some useful relationships have been derived in 
Appendix V. To obtain literal closed-form stability criteria we shall 
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restrict our flexible appendage model to lie in a plane containing the 
center of mass and normal to the spin axis. 



Figure 13. Restricted Model. 

Following the definition of our particular case, we have to 
consider 



Under this last assumption, the three scalar equations expressing the 
rotational motion of the rigid core are given by considering equation 
(1) in Appendix V for this particular case or: 


0 * Afi i - “ SS “2 (B - C) * 2 A 23 “SS + 2A 13 “SS 


/ r A + / r 2“: 


dm 


0 B(i 2 + W 1 ^ A-C ^ + 2A 13 + 2 A 23 + 2 U 3 dm ” Jj l U j 


dm 


0 " C “3 + 2A 23 





dm 


where we used the fact that the matrix 1^ has been chosen to be diagonal 
and where the angular velocity vector at the "steady-state" deformation 
has only one non-zero component, Wgg, along the spin axis. 
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Under the assumption of a flexible appendage lying in a plane 
orthogonal to the spin axis and passing through the center of mass, the 
kinetic energy may be written as: 


T = T 1 + T 2 , 


where 


T 1 ■ 1 *“1 + 1 Bu 2 + i f *. 3 d “ - “ss “l ^ r i u 3 d " - “ss “2 J a -2 “3 


i f*l d “ - “SS “l / F 1 U 

+ “l / r 2“3 dm - “2 J}lS d ” - -tit If" 


A 


dm 


T 2 " 2 C 4s + Ca> 3 “SS + 


■1.2.1 f. 2 . , 1 f. 2 

2 C “3 + 2 7 u l dm + 2 / u 2 *» 


+ 2 4 / (u l + 4 dm + <4 jf (r 2 U 2 + F 1 u l } dm 


+ 2U 3 U SS jf (r 2 U 2 + r i U l )dm " (W 3 + W SS ) A 2*1 

•'A •'A 


dm + 


+ <M 3 + “ss> j(' r i' 1 2 dm - “SS J[ (i l u 2 ’ “2 V d " - TM f" 1 

A dm l 1 + jr If A dm l f A d “l “ss - f A -1 f A 


1 2 
2,ytf W SS 


A H f/ u 

[A HMArl 


dm 


dm 


0). 


SS 


If we, now, notice the fact that 
2 A 13 ‘ " f* 1 u 3 dm 


2 A 23 ' ' f * 2 u 3 dm 


2 A 33 ' 2 J [ (T 1 U 1 + r 2 u 2 ) dm ’ 
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it becomes obvious by looking at the previous set of equations that the 
equations describing the rotational motion of the central body separates 
The equations expressing the change in the angular velocity components 
0 )^, 0)^ are coupled with the transverse vibration of the flexible append- 
age and form what is called the "wobbling" motion. Also the third com- 
ponent tOg is coupled only to the components u^ and u^ of the deformation 
and constitutes the "spinning" motion. 

If the same separation is being observed in the deformation 
equations, we could conclude that the general motion is divided into two 
different types of motion One, where we observe a coupling between 
transverse vibration and the first two components of the angular veloc- 
ity, and the other where we find a coupling between the in-plane vibra- 
tion and the angular velocity component along the spin axis. But we 

.1 

just observed that for our particular case (T^=0) the kinetic energy can 
be expressed as the summation of and where in appear only the 

variables of the wobbling motion and in appear only the variables of 

the in-plane motion. 

As a result the separation suggested before is effectively 
accomplished for our particular case if the potential energy of deforma- 
tion is such that U, can be written as the summation of U. and U_ where 

d 1 z 

the aforementioned separation prevails. 

We will now assume that the in-plane and out-of-plane stiffness 
elements are uncoupled which is the same as expressing that the matrix 
K of the energy of deformation is partitioned such that 


K 11 ^12 


K = 


K 


21 

0 


K 


22 
0 K 


0 

0 

33 


2-170 



Under this last assumption, the motion separates and the Hamiltonian 
itself can be written under the form: 

H c = H 1 + H 2 

where H^contains the variables appearing in the wobbling motion and 
those constituting the spinning motion. 


H -IK 
H 1 “ 2 2 
C 


+ 


BC 

1_ 

2A 


A [^1 “3 

»3 - - -K * 2 (!?) - -K » 2 (i?) 


12 



2 


fv 

u„ dm 


B 

[A 2 

O 



0 2h 2 

AC" J A ‘ 1 U 3 


/ r i u : 


dm 


[X r 2 S d "]-a[X r i 4 3 dm J + l j} 

- ~hi \f* 3 dm + 2 / K 33 u 3 dm • 


dm 


H 

2 c 2 


+ 


■' [f^i r i + u : 

|6 


2 

dm + 


r 2 ) dm 


_ h_ 


2C 2 


2 


iu, dm 


^L_ 

\ 1 
A 


2C 


;) dm 


+ ”2 
“2 - r 2 d l> 


dm 


. 1 f,.'2 . .2. , 1 

+ 2 / (u l + U 2 ) d ” - 2 M 

A 

[X 4 1 dm 

2 

+ 

X 4 2 

2 







+ 


i M 


K i;l + 2u^ u 2 K 12 + u 2 K 22 ) dm 


We have proven that H = 4- H 2 where the variables appearing in 

are uncoupled to those appearing in H 2 through the equations of the 
motion. We conclude that if both and H 2 are positive definite, we 

are assured of asymptotic stability, and if either or H 2 is sign 
variable or negative definite, the nominal motion is unstable. 
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Experience suggests that stability criteria emanating from H~ are 
of little practical importance; the variables in describe deformation 
in the plane of the appendage, which is normal to the spin axis, and 
generally those criteria express simply that the structure has to be 
sufficiently stiff to avoid destruction by "centrifugal forces." We 
will thus consider the case where the stiffness elements orthogonal to 
the spin axis are assumed infinitely large (u^ = u^ = 0) so that the 
structure is allowed to vibrate only in the u^ direction. The cited 
assumptions allow stability criteria extracted from to establish 
characteristics of the entire system. To formulate the restricted 
Hamiltonian in a more useful form, we have to introduce some new notions. 

The class of problems covered by all the previous assumptions 
includes as particular cases all the problems and particular structures 
we consider in chapters 2 and 3. We now consider that the deformations 
are expressed in terms of modes of the structures; this last consider- 
ation will enable us to be more specific about the stiffness . 

The deformation will be expressed in linear approximation as a 
linear combination of modes 


u 3 t 
V=1 


3 <f> 
v 


We will become more explicit later, concerning the exact definition and 
choice of modes. Let us just remark that we introduce in the last defi- 
nition N new independent generalized coordinates — the 3^'s. The expres- 
sion cj>^ is related to the shape of the "displacement" in the mode V and 
is thus defined as being a function of position only, the real 
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deformation for this mode being then equal to the product of <j> by the 
harmonic function of time 

Introducing the last definition into the energy of deformation, we 
obtain: 

U d * L 0 + “ss A 33 + \ fXf *v K 33 £ S li dm • ' 

A v y 

For our particular case, = 0, so the previous expression may be 

written: 

u d “ L 0 + 1 8 1 16 

where we define $ as being the N x 1 vector 



and K as being the symmetric N x N matrix defined by 
K - [k vp '■ X+v K 33 *u d " ■ 


Expanding the deformation u^ in terms of the previous N modes, the 
equations of deformation are readily obtained in their Lagrangian 
f orm : 



where T and U, are the kinetic and potential energies, 
d 

For this particular case = 0, u^ = 0, u ^ = 0) the kinetic energy is 
reduced to the following combination: 
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1 2 , 
iAWi + 2 


■f r u 

x 2 + i ^SS W 2 A 1 

+ 7 BW 2 + 2 J* 3 




dm 


+ W 


X 

ff ti~ dm 2^\.A 3 

A i z *° ' “ 2 j 

f 1 3 

A 


+ -» CW SS 


1 


+ i CUE, * 

+ CW 3 W SS 2 3 ^, pS 6 , *e have, 

. o£ the variables 
a in terms or 

Expressed 


= & 1 4 > 


vector 


the S * 1 

where ve a—- 

4 : Ta ... *** M 


4 > as 


we define 

1- 4 [h " 


Vie also 


have 


/ -3 ' 1 


4 > 4 . T & dm = & 1A 


[/w 1 4 • ■ eT “ 


f r i U 3 

A 

1 y 

fy 


T A- <j> dm = & ^1 

_ f T A 1 

dm - X 1 

£ B 1 A ? 


dm - j[V** 
•'A. 


rr> R 1 fT 0 4 > d ® 

?<b dm " P J A 2 


dm 


- f A 2 



x h 


,. ric matrices 

symmetric 
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/' 

A 


T A 
<p<p dm = M 


4> dim 


[f “1 ! 


m 


and the N x l vectors and by 


/ r i 


4> dm = A 1 


I 1 


r 2 <p dm = A 2 


The kinetic energy then becomes: 

T = ~2 ^1 + "2 BW 2 + *2 ^ ~ “ss “l ® A 1 " “ss “2 B A 2 

+ ^ f A 2 ' “ 2 8 1 A 2 - 2^,8 T m S + j + c™ 3 o. ss + j Cu* . 

The deformation equations can, now, be written in a completely explicit 
form: 


M 


6 - e + K e ■ “2 A i - *1 A 2 - “i “ss A i - “2 “ss A ; 


If we add to this set of equations the rotational equations of the 


central rigid body under the assumptions = 0, u^ = 0, u.^ = 0 we get 


0 ■ *“1 - “2 “ss (B " C) + “ss I • 2 “3 

■'A 


f T 2 U 3 * + jC r 2 Sj 


dm 


0 ■ b “ 2 + “i “ss (A - C) - “ss 


Si 1 u 3 dn ' fh S 


dm 


0 = Cto„ 


or, introducing some of the previous definitions: 

0 - Au x - w 2 W ss (B-C) + Wg S 3 T A 2 + g T A 2 ✓ 


0 - Bl 2 + Uj Uss (A-C) - 4 e T Aj - e T Aj 
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In order to study the stability of the previous system, we now 

need to define how the modes have to be computed. In order to do so, 

we have to keep in mind that we are interested in obtaining a system 

which is suitable for a mathematical analysis. Similarly, for the 

common procedure used in the "hybrid-coordinates" approach of flexible 

22 

space craft and developed to some extent by P.W. Likins, , we have to 
use some practical simplifications in order to obtain the solution. 

Even though, in the equations of the deformation, the forcing 
terms, linear in the U)'s, are in fact coupled to the flexible appendage 
deformation variables, 8^, through the central body motion, we define 
our eigenmodes as being obtained only by the consideration of the homo- 
geneous part of the equations in 8's. Unfortunately, the modes are 
still coupled, through the motion of the center of mass, motion result- 
ing from the appendage .deformation. The center of mass motion is imbed- 
ded in the matrix m of those equations. In order to decouple the equa- 
tions, we must add the assumption that the eigenmodes are defined when 
the system's center of mass is assumed fixed, but the structure deformed 
by the deformations induced by the steady-state spin. 

This last assumption could appear too restrictive but the coupling 
between the rotational equations and the deformation equations is 
obtained through the vectors and A^. For a symmetric structure, with 
homogeneous physical characteristics, only antisymmetric modes introduce 
nonzero elements in the vectors A^ and So, we remark that the cou- 

pling between the two sets of equations is obtained by means of anti- 
symmetric modes and those, precisely, do not affect the motion of the 
center of mass. 
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We could also, as several authors did, consider that the 

rotational motion of the central body of a symmetric vehicle is only 

affected by the anti-symmetric modes, and from there consider in the 

deformation variables, only those corresponding to anti-symmetric modes. 

1 ** 

With this approach, the term m 3 would not be present in our equa- 
tion. 


In what follows, we focus attention on symmetric vehicles, in 
order to take advantage of the possibility of separating eigenmodes into 
symmetric and anti-symmetric classes. 

With the modes completely specified, we know that if we let the 
structure vibrate in one of its modes, the frequency of the vibration 

■ • s. 

2 

would be given by 0 )^, the eigenfrequency corresponding to the aforemen- 
tioned mode, and we will anticipate a periodic motion (no damping has 
been introduced when deriving the eigenmodes and eigenfrequencies . ) For 
an orthonormal system of modes, the matrix M becomes the unity matrix E 
and we thus obtained as equations for the deformation variables 


E 3* + K 3* = 0 

where the variables 3* have been substituted for 3 V in order to differ- 
entiate them from the equations obtained with the forcing terms. 

This last equation may also be written with the substitution 



where the K matrix becomes a diagonal matrix the elements of which are 
composed of the squares of the modal frequencies. Those frequencies are 
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the loaded natural frequencies of the appendage accounting for the fact 
that the stiffness elements are modified by spin. That is, the tension 
in the structure induced by spin alters the appendage stiffness charac- 
teristics . 

With all the previous assumptions and new definitions, we are in a 

position where stability criteria could be extracted from and where 

the stability characteristics of the entire system could be established. 
Introducing the following definitions 

[/ r i “3 -to'] 2 ■ \f aJ 2 £ b 1 y 
[/ r 2 U 3 dm ] 2 ■ [ bI A ] 2 ' ^ V 

where the matrices II ^ and II^ are the N x N symmetric matrices satisfying 


1,1 4 A x A X T 


A T 

n 2 = ^2 A 2 


the general expression for the Hamiltonian becomes in a matrix nota- 


tion 


_ 1 2 2 C-A 1 2 2 C-B 
H l~2 h e 2 AC~" + 2 h 8 1 BC 


h 2 ° h 2 0 

-Air 6 A i + tst *\ 


* i * T 


1 *T 
+ j B 


2 2 

n 2 + n + n. 

AC BC 


e - r n, - ~ n, - ^m 


3 


A 2 B 1 2 Jl 

We can now write the stability criteria as being given by the 
positive definiteness of or 


H > 0 
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\ 


where it is to be understood that > 0 means positive for all values of 
0^ 0^, 3 and g in the neighborhood of the origin 0^ = 02 = 8 V =3 V = O, 
except equal to zero at the origin itself. 

If we, now, group the eigenmodes into two categories, the first 
one including symmetric modes and the second one regrouping the antisym- 
metric modes, it is straightforward to see that in the development of H^, 
the symmetric modes contribute to H^, only through 


1 o T 

2 


~ 2 

„ ' 1 A T 

1 


+ 2 

E s iiM m s 


where the subscript "S" has been introduced to distinguish this part 
from the rest of the development of H^. This last observation is in 
fact a consequence of the remark made earlier: only anti-symmetric 

modes have nonzero contribution to the vectors and A^. It is evident 
that the first term is always strictly greater than zero. In order to 
prove that the second term 


1 A T 

2 


E s - 


lOi m s 


is always strictly greater than zero, let us consider the following 
expression 



where u^ s stands for the symmetric part of the deformation u^. Applying 
the Scharwz inequality. 


P(v,w)| 1 2 < P(v,v) .P(w,w) 


to the term 


[As *"] 2 - 


we obtain 
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./“ 3s dm J s $ s ' ° A X“ 3s 


dm 


where m^ stands for the mass of the flexible appendage only. 

J , 

As a result, we have 

2 


/ U 3S * - 2^ _/ u 3S ^ (- 1 * -ht) / U 3S 


dm 


but 


m. 


UK- m A ... V 2m B 
2^r ^ 


i - 


2*A( 


> 0 


where«y^has been written as the summation of m., the mass of the flexi 

A 

ble appendage, plus m_, - the mass of the rigid core. 

From the last observation, we conclude 

1 2 


/ U 3S dm - 2 hi [{ u 3S d “] 


> 0 


Expressing the deformation u^ s in terms of the symmetric mode, we have 


N 


3S 


i e vS ^VS’ 


V° 


The last inequality becomes: 


3c 


E„ - 


s 2<M m s 


3g > 0 


We conclude from the last relationship that the matrix 

E s ' 


is positive definite and from these that 


1 • T 
~2 


E s 2^M m s 


3g > 0 
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So we just proved that the variables of deformation corresponding to 
symmetric modes separate from the remaining terms and that their con- 
tribution to the Hamiltonian is always strictly greater than zero. 
The stability criterion becomes finally 


*}(&)- 


h 2 e h 2 e 


4e T 


1 *T 


b 2 + n + n 

AC BC 


E - i ” 2 - s n i 


8 > 0 


where only the contributions of anti-symmetric modes have to be taken 
into account. We observe that in the last inequality the matrix m is 
no more present. 

In the last development of the Hamiltonian, the last term is 
uncoupled from the remaining terms. We conclude directly from here a 
requirement for the asymptotic stability of our system, i.e., the matrix 


s-l'V-iK 


has to be positive definite. The author believes that due to physical 
properties of realizable systems, the previous matrix is always positive 
definite. He was however unable to prove this last statement for a 
general case, but for a wide class of flexible spacecrafts, it is veri- 
fied. 


For the class of flexible spacecraft where the rigid core presents 
an inertia symmetric with respect to the spin axis, or where the princi- 
pal axis of the rigid core coincide with those of the. flexible append- 
ages, the principal axis of the entire system are the same as the prin- 
cipal axis of the flexible appendage alone. If, together with this last 


2 - 1.81 



property, the flexible appendage presents a symmetry not only with 
respect to the center of mass but also with respect to the principal 
axis of inertia, then the previous statement is indeed true. 



Let us consider the following example where all the previous 
requirements are met. The figure represents the body when spinning at 
steady state. From previous consideration, we know that only anti- 
symmetric modes are of interest. We, now, consider in the flexible part 
two elements of mass dm of coordinates (T^, and (-1^, r ) when at 
steady-state deformation. Because of the assumed symmetry of the flex- 
ible appendage, such a pair always exists. 

We now decompose the anti-symmetric part of the deformation — our 
main interest — into two parts: 


U 3A ^l’ = U 3A ^1’ T 2 ) + U 3A ^1’ T 2* 


r r 

1 2 

where the deformations u^ and u^^ are defined in the following way 
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U 3A <V r 2> + U 3A '- r i> r 2> * 2 “3A <r i- r 2> 

and u 3A r 2 ) - u 3A 1 (r x , r 2 > - uj 2 <r r r 2 > 

We, also, have then 

U 3A (_r i' T 2 > " U 3A *" r i’ r 2 ) + U 3A (_r i’ f 2^ 

with uj 1 (-Fj, F 2 ) - U3/ 1 (r lt r 2 ) 

and “sa" 2 <‘ r i> r 2> - - “sa' 2 < +r r ; r 2> 

We then have for the two elements of mass dm, of coordinate (T^, T^) and 

<- r r r 2> 


F 1 U 3A (r i* r 2 ) r i U 3A ( ~ r i* V dm 


( r i U 3A " (r i- r 2> + r i "J 2 <r i> r 2> - r i “Sa" 1 <- r r V 

- r i U 3A <- r i> r 2>) d “ ■ 2F 1 ^ < r i> V d ” 
and similarly 

r 2 U 3A < r i r 2 ) d ” + P 2 U 3A <- r r r 2> dl " ' 2F 2 U 3A ' < r i’ r 2 > d " 
From there, we may conclude 


/ r i “3A * ’ f/l “JA^ 
[ T 2 U 3A ' / r 2 U 3A 1 


dm 


Let us now consider the expression: 
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/ r 'l d "l[/ r 2 d "l f/ u2 3A*]-[/ r i 


dm 


Jh “3A 


dm 


“ 

" 

r 

J£ 

,2 j 
2 dm 



— 



II 1 “ 3A 


dm 


We have 


/“ 3 A d " * ') + (“ 3 A 2 ) 


r i M 

+ 2 U 3 A U 3 A ) dm 


(4.6) 


‘ K "* 1 ) 2 * + / f ^ 2 ) 2 * 


r, r. 


f -J i _ 

where we observed /u_ 4 u_. dm = 0 

3 A 3 A 

Introducing the equality previously found, the expression (4.6) is also 
equal to 

r\v 2 rl r„V2| 




dm 


I: 


Y 2 2 dm 


/(•»■/ * * 


dm 


dm 


h J 

U 3 A dm 

2 

/ 

f T 2 dm 

•'A 

m 




/ r i U 3 A 


dm 


Applying the Scharwz inequality, we have simultaneously 



) 

% 

■V, 

to 
i 


h dm 


> 

X r2dm . 

X(^) 2 r 

> 


l y 2 u 3a 1 


dm 


3 A 


dm 


Introducing those two inequalities into the development of the expres- 
sion (4.6), we conclude that for the particular case into investigation, 
we have for every anti-symmetric deformation u.^, the following inequal- 


ity 
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d "l [X U 3A *1 - [ft dm ] \X r 2 U 3A 


dm 




J* 

f T? dm 


/r, u„. dm 

7 A 2 


7 a 1 3A 





> 0 . 


With this last inequality, we are now able to conclude that for every 
anti-symmetric deformation u,^, the following inequality is true 


AB 


/* 2 

/*r 


2 

1 

- 

dm - B 

J 

J k 

2 U 3A dm 

- A j 

[A 

1 U 3A d " 


> 0 


(A. 7) 


where A, and B are the inertia of the total system. By definition, we 
have 


A = A 1 + (v 2 dm 

*A 1 


* b1 + I T 1 


dm 


where A"* 1 and B^ are the inertia of the rigid core. 

Expanding the left hand side of the expression (A. 7), we have 

1 + fr\ dm A 1 + fv\ dm dm - B 1 + fl\ dm fv » JA 


dm 


ll+ / F 2“” Zb “3A 


dm 


The use of the Scharwz inequality and the consideration of the previous 
inequality enables us to conclude the inequality (A. 7), after the obser- 
vation that for any physical realisable system we have 


■' A 


dm > 0 
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So we just proved that for the class of problems where the principal 
axis of the flexible appendages and of the entire system coincide, we 
have 

ro.Tr 12 T /• 12 

> 0 


AB 

«• 

f u 2 dm - B 

4 3A 

'4 

'l U 3A 

2 

- A 


"l U 3A *■ 




_ 



_ 


Expressing u^ A in terms of the anti-symmetric modes, we have 
N, 


“ 3A e vA * VA 


or 


ab e T Eg - B3 T n 2 3 - A3 T n 1 3 > o 


Dividing through by AB, we obtain finally: 


3 


E - X "2 - 5 n i 


3 > 0 


We thus conclude that the matrix 


E - S " 2 - 5 "i 


is positive definite. 


The requirements of positive definiteness of are now reduced to the 


consideration of 


H i * i h2 0 2 (s*) ♦ \ h2 9 i (fr) - ^ ® T A i 


h 2 0 

+ A 2 


+ \ 3 t ft 2 + n + ^ n 3 > o 

AC BC 


In matrix notation, we have 
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The N deformation variables can now be truncated to a single mode. 


identified by index 1; thus the total number of coordinates is reduced 
to three. Accordingly, the N xi vectors and A^ reduce to scalars. 
Similarly the N x n matrices and JI^ reduce to scalars. 


Implementing the above truncation allows the stability condition to be 
written as 
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The sign character of the above quadratic function is determined by 
testing the sign character of its corresponding symmetrix matrix. 
Sylvester's Theorem assures that the necessary and sufficient con- 
ditions for the previous matrix to be positive definite is that all 
principal diagonal minors be simultaneously positive. If this test 
fails, H| is not positive definite and is either negative semidefinite 
(or sign variable) , implying instability, or positive semidefinite. 

If we exclude this latter limiting case (as for an axisymmetric 
vehicle, with C = B or C = A) , necessary and sufficient conditions for 
asymptotic stability of the restricted planar appendage model are given 
by : 


(C-B) 

BC 


2 

hi 

2 

hi (C-B) 
2 BC 


.(C-AJL 

AC 


> 0 


> 0 


ill (C-A) 
2 AC 


r 2 

w i 


tL h 2 + jL a 2 
2 A 1 2 Jl 2 


2AC 


2 A 


_2 

BC 


2BC 


h 4 A 2 


(C-A) 
4 BC AC 


£ A. 

4 2 2 
A C 


> 0 


The combination of the first two conditions, as predicted by energy 
sink methods, requires that the spin axis be the axis of maximum 
moment of inertia, i.e. 


C > A and C > B . 


The requirement of the third condition emerges and can be written 


2 ^ A 1 

W 1 2 

C^(C-A) 


h 2 A 2 


C (C-B) 
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By replacing h by its zeroth order approximation C Q , where Q is the 
nominal spin frequency, the above condition simplifies to the following 



Thus a stability criterion arises which explicitly bounds the first 
modal frequency of the structure. 

To conclude this derivation and to demonstrate the interrelation- 
ship between Chapter 4 and the preceding chapters, it might be useful 
to work out a simple example. An elementary but meaningful example is 
given by the beam pair configuration, as in Figure 15. 



Figure 15. Simple Example. 


For this particular structure, we have I* = Y^ = 0. The 
stability criteria for asymptotic stability of spin are given by C > A 
and C > B and also 
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The consideration of = 0 reduces the last criterion to 



Let us concentrate on the first mode, as our single mode. From our 
previous derivation, we conclude that only the anti-symmetric modes 
are relevant to our study. From the derivation accomplished in 
Chapter 2, we have the following results, when considering the first 
mode. 

(r) - l + 2 - 12 ^ • 

also 

4> 1 = h Q (x) + Je h 1 (x) +... 


For the first mode, h^(x) represents a linear function k^x, so we 
obtain for the following expressions 


A 


1 



dm 



dm 


where the anti-symmetry of the deformation has been used. For a uni- 
form beam, we have dm = y dx, so the previous expression becomes 



2y 



hp(x)dx + 2y 


/e 



x h^(x) dx +. . . 


For the first mode, the orthogonality relationship between h^ and h^, 

cancels the second term and we are left with 

L 

A 1 = 2P /V 2 dx = 3 y k 0 L3 ' 
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The modes have been defined in such a way that 


ft 


dm = 1 


or 


so 


u 


2 I k“ x“ y dx = 1 

'o 


2 = 3 __L_ 

0 ~ 2 3 

yL 


The final expression for is then 


.2 2 T 3 2 T 2 

A 1 ’ 3 " L ’ 3 V 


where m^ = p L stands for the mass of one beam. The final expression 


for the last criterion is thus 


1 + 2.12 > 4 ^ 


3 C-A ’ 

If we observe that the inertia of a uniformly distributed beam about 

1 2 

the core is given by — m^L , the last criterion is also 


1 + 2 12 Je > ( Inertia the two beams) 

C-A 


This last result has to be compared with the stability criteria 
developed by L. Meirovitch and R.A. Calico, Reference (23), where they 
obtained stability criteria for a spacecraft characterized by a rigid 
core having attached to it flexible booms. Their results, when applied 
to radial rods only, become 

( 

which shows a perfect analogy with the result of Chapter 4. 


^1 | - (Inertia of the rods) 

ny c-a 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 

In the work developed, two different problems have been 
approached: A modal analysis, covered in Chapters 2 and 3, and an 

attitude stability study, presented in Chapter 4. Each part represents 
a different aspect to investigate for the preliminary design of flex- 
ible appendages to be attached to a spinning spacecraft; they comple- 
ment each other under the general question of stability. 

In the first and more extended part, general inferences con- 
cerning the magnitude and character of the influence of spin on the 
natural frequencies and mode shapes of some rotating structures are 
done. One of the objectives of this part was also to address ourselves 
to the level of sophistication that a derivation using a continuous 
model for an elastic appendage attached to a rigid base which is 
constrained to rotate with a constant angular speed about a body 
axis, fixed in inertial space, could bring us. . We decided to try to 
get away from more standard procedures — based on energy consideration 
and transforming the eigenvalue problem to a minimization one — by 
approaching them by a method unfamiliar to those particular appli- 
cations. The method of matched asymptotic expansions has been proved 
to present a very powerful tool for the particular problems we looked 
at. Even though, we didn't attempt to consider all the meaningful 
structures where a singular perturbation method could be appropriate, 
we believe that we cover several problems where the usefulness of the 
method has been found. In contrast with the classical Rayleigh-Ritz 
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or Galerkin approaches , the accuracy of the perturbation method is 
determined by the number of terms taken in the various expansions and 
may be improved step by step by the investigator without the necessity 
of repeating the analysis with an augmented number of terms in the 
expansion of the solution. 

The literal stability criteria developed in the last part of this 
present work represent anticipated results in the search for closed 
form conditions for attitude stability of spinning flexible spacecraft. 
They are, in fact, a natural extension of those derived in a similar 
study where another way of modeling the flexible part of the space- 
craft was adopted. They represent (except for a few recognizable 
singular cases) necessary and sufficient conditions for stability for 
any spacecraft characterized by the planar appendage model, such as 
a spacecraft containing solar panels and/or radial booms. More pre- 
cisely, they are necessary and sufficient for that portion of the 
system representative of the wobbling motion. In terms of the 
composite motion, these conditions can only be classified as necessary. 
Also this analysis does not reveal the system behavior at any time, 
but only the quality of the motion in the neighborhood of the dynamic 
equilibrium configuration consisting of the body spinning at a high 
angular velocity about the spin axis. 

Although the results of this study could be used in preliminary 
design of flexible appendiages to be attached to spinning spacecraft, 
we have to be very cautious in extending the previous work to more 
general types of configuration. We mentioned earlier that a modal 
analysis of an idealized elastic structure on a rotating base requires 
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the derivation of the linearized equations of small vibrations of the 
mathematical model from its steady state of deformation, induced by 
spin, and the transformation of these equations into a system of 
uncoupled equations of motion representing each normal modes. For a 
continuum model, the number of normal modes could be infinite. For a 
completely general structure, this step represents a major obstacle. 
Even in the very elementary structure consisting of a spinning central 
body with a cantilevered beam oriented in the direction of the spin 
axis, the influence of the basis rotation is manifested in the form 
of centripetal acceleration and Coriolis acceleration. The latter 
term couples the equations, and provides an obstacle to modal analysis 
except for the particular case where the linear density is constant 
throughout the beam and both transverse inertias of the cross section 
of the beam are the same constant. And, even in this last case, the 
modes are given by a complex function. 

It is to circumvent the last obstacles, that several authors 
consider in their modal analysis an attractive alternative, employing 
for the flexible appendage the coordinates which are normal mode 
coordinates when the basis is inertially fixed. Although this last 
shortcut may be acceptable when the motion is a small perturbation from 
a rigid body displacement at a relatively slow spin, it is more diffi- 
cult to justify for a general case. They also look, at the outset, 
at the spinning stiffness matrix K as composed essentially of the 

A 

summation of the nonspinning stiffness K and a matrix representing the 
contribution of spin — sometimes referred to as the geometric stiff- 
ness. But here again, they face another source of problems. In 



general, the eigenmodes of the nonrotating structure are not orthogonal 
for the rotating structure; as a result, the spinning stiffness matrix 
is no longer diagonal for this last choice of eigenmodes. As a result, 
the nonspinnihg stiffness K is well defined and made of the unloaded 
natural eigenfrequencies along the diagonal, but there is ho real 

A 

physical basis to determine the correction to add to the last K matrix 
in order to obtain the complete stiffness matrix and this last effect, 
configuration dependent, may or may not seriously modify the stiffness 
property of the structure. 

In view of those several difficulties, we might conclude that the 
development of elastic continuum models could yield very useful results 
when applied to a small class of special cases, but this approach 
lacks the general utility and tractability of distributed-mass finite 
element models, since in the latter case, the governing equations are 
always linear, constant coefficient ordinary differential equations. 
Finally, the generalization of the previous results can only be con- 
sidered for that small class of problems of elastic models where the 
steady state deformation can be solved first, and where the vibration 
equations are obtained by a linearization procedure; and even in this 
case implementation of the general theory may be very difficult. 
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APPENDIX I 


ORTHOGONALITY OF THE MODES 


Our objective is to prove that the values previously found for 
2 2 

Xq and also satisfy: 


Y 

J.-Y 


<j) 0 dx = 0 for a ^ B 

p 


The asymptotic expansions for the inodes are given by: 


<j> a (x,e) = h 0a (x) + 71 

h la^ + ’ “ 

<j>g(x»e) = h Q g(x) + 7e h 1 g(x) +... 

The problem becomes equivalent to: 

h„ h_ Q dx = 0 

, Oa OB 
-k 

y 

j (h Oa h lB + h 0B h la )dx ° 

-k' 1 

where h n and h, _ satisfy 
Oa la 

[(l-x 2 )hy ' + X^ a h 0o - 0 

[a-* 2 )h ia j' + h la . -4 h 0a , 

and similar expressions for h^g and h^g. 

The first equality is easy to show. Multiplying the differential 
equation defining hg a by hgg and integrating over the whole span, we 
have : 



AI-1 



-1 


f 
/ -1 


t . (1-x2)h 0a 1 ' h 03 dx + X 0a 


-k 


-1 


f 
/ -1 


h. h dx 

Oa 03 


= 0. 


-k 


After an integration by parts, the use of the transcendental equation 
2 2 

defining Xq^ and X^, gives us: 


-1 


-1 


/ - 

-k -1 


x2>h Oa h oe dx ■ A 0a 


/ * 

i- 1 


0a h 03 dx 


f 1 

By symmetry the previous equation is also equal to X nD / h n h „ dx, 

Up OCt Up 

from where the first equality results. 

The equality 
, -1 


/ / h 0a h 13 + h O0 h la )dx 
-k _1 


is more difficult to prove. The differential equation defining h^ a 

is first multiplied by h no and integrated over the whole span. The 

Up 

orthogonality of h^ and hg^ is then used and we have: 


-1 


f 

-1 


t(l-x^)h’ ] h nQ dx + X^ 

1 ' la 03 0a 


-1 


f 

/ -1 


h. h^o dx = 0. 

la 03 


-k -k 

Integrating the first term by parts and using the boundary conditions 


on h no , we have 
UP 


-1 


-J a- 

V 1 


x2)h ia h 0S dx + A 0a 


-1 


/ 

J -1 


h., h_.„ dx = 0. 

la 03 


-k 
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Integrating by parts once more, and after the use of the differential 


equation defining h no , we have: 

Up 


-(! X ) h ^ gh la 


-1 


-1 


-k 


-1 


+ (X 0cf W 


• i ,- 1 


, h- D dx 

la 03 


= 0 


A similar expression can be obtained starting from the differential 


equation defining h 1G or 

lp 


< 1 - >h ia h ie 


-1 


-k 


-1 


+ < X 0a - *0B> 


-1 


L 


h „h dx = 0 
13 0a 


-k 


Adding up those two relations, we have 


a -^ b 'o*\e 


-1 


-k 


-1 


- 


-1 


-1 


+ (X 0a " X 03 ) 


-k 

-1 




(h 13 h 0a + h la h 03 )dx ° ' 


-k 


The two first terms can be written: 

(1 - k ' 2) [h ta <k ' 1,h 1B (k ' 1> - h ia ( - k " 1,h lS ( - k ‘ 1, - h i6 (k " 1>h la (k ' 1) 

- h i B ( - k ' 1)h lo« ( - k ' 1)1 ' 


and by direct substitution, this parenthesis comes out to be 
identically zero. We then conclude 
,-1 


f 

A ' 1 


(h l3 h 0a + h la h 03 } dx = °* 
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APPENDIX II 


SOLUTION TO THE REDUCED EQUATION OF THE TAUT MEMBRANE 

In the consideration of the second linearly independent solution 
of the previous equation, two cases have to be considered. 

CASE A : Both roots of the indicial equation are the same and 

equal to zero or 

1 - y = 0 or y = 1 * 


This case corresponds physically to m=0 or to radially symmetric vibra- 
tions of our clamped disk; To obtain the second linearly independent 
solution, we shall assume that 

CO 


H 

c 


c , , n+c 

x + /. b x 
, n 
n=l 


where c is retained as a parameter. The last expression does not 
satisfy the differential equation defining H^, but if we write the 
differential equation defining H^ or 


d^H 


x(l-x) 


0 1 ^0 

7? + u- 2 *> + r 7 H o * 0 

dx k^ 


under the form 


- 0 

where SB represents the linear differential operator 

0 2 

2 A 

\ - x(1 ~ x) 7T + (1 - 2x) 57 + 4 7 ’ 

0 dx k^ 


we have 
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A 

5? [H ] = c 2 x C c 2 x 2 +b 1 (c+l)cx C 4 ^ x C -cx C + b n (c+1): 

0 k l 


n=2 L 


b (n+c) (n+c-l)-b .. (n+c-1) (n+c-2)+b (n+c) 

n n-J. n 


-2 b n (n+c-1) + b 1 
n-i 4k 2 n-1 


n+c-1 


x 


We are now able to make all but the first term of the right hand mem- 
ber of the above equation vanish without selecting c. To do this, we 
must take 

A 2 

-c 2 + b (c+l)c + - c + b (c+1) = 0 

1 4k^ 


and 


2 o 

b (n+c) - b . (n+c-1) (n+c) + b . — - 0 . 
n n-1 n-1 4fc 2 


So far, we have found an expression 


H = x + 
c 


IX 


n+c 


n=l 


that reduces the previous expression to 

K [H i = c 2 x C 1 , 

H 0 C 


no matter what value c takes on. But this last expression not only 
shows that the case c=0 is solution of the equation defining H^, but 
also (dH^/dc) evaluated at c=0. The solution, found before, corres- 
ponds to the case where c=0 so our interest lies mainly in finding 


dH 

c 

dc 


for c=0 , 
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or 


dH 


c 2‘ V"' , n+c, 

— = x log x + / i b^x log x + 

n=l 


3b 


OD 


n+c 


n=l 


We finally obtain the second linearly independent solution 

co co 

3b 


H 20'( 1 + £' b n| “j lo < ,+ 2 IT 


n 


x 


c=0 


or 


E dC *n 
■57- 


n=l 


n 


c=0 


The complete solution is then given by: 

V' ^ 

H 0 = a 0 F ( a » e »Y.x) + b 0 F(a,3,y;x)log x + b Q 2-# 


n=l 


n 


x 


c=0 


Due to the singularity at the origin, we have to take b^ = 0, and the 
solution reduces to: 


H 0 = a o F ( a »^>Y;x) .. 


Case B: In this case, the roots of the indicial equation are 

different and from there, m is different from zero. In order to 
determine the other solution, two different cases are to be considered 
depending on the values of a and 3. 

a) a or 3 is an integer between 1 and m. Let us try a particu- 
lar solution starting with the index -m, or 


OO 


H 


02 


-E 

n 


n-m 

a x 
n 


By substituting the last expression into the differential equation 
defining H^, we obtain as a coefficient for a^, the indicial equation 
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which is satisfied because the particular solution starts with the 
index -m. The other part gives us the recurrence relationship between 
the a^s or more explicitly 

= (n-ct) (n-g) 
n (n-m)n n-1 

By looking at this particular solution, and knowing that m is different 

from zero, but taking on positive integer value, the only way such a 

solution could exist is if the numerator of a is also equal to zero, 

m 

but 

_ (m-a) (m-B) x . (m-q-1) (m-3-1) 
m (n-m)m (m-l-m) (m-1) ’’’ a 0 

The last particular solution is then possible only if a or 3 satisfy 
one of the following equality: 

a is an integer between 1 and m 
3 is an integer between 1 and m. 

If a or 3 is an integer between 1 and m, one of the a^s, in the pre- 
vious recurrence relationship becomes zero and the particular solution 
is represented by a truncated serie. This particular solution com- 
bined with the result found earlier represents the general solution of 
the equation. Unfortunately, this last solution has to be dis- 
regarded for. its singularity at the origin. 

b) If a and 3 are not integers lying between 1 and m, then the 
second linearly independent solution for can be taken under the 
form 

h 02 = H oi lo e x + H 0c • 
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Introducing this last expression into the differential equation defin- 
ing Hq, we are left with 

d 2 H n dH dH m 

x(l-x) — ~ + ty- (0+0+1) x] ■ - a& H_ + 2(l-x) — — 

/ dx uc dx 

H 

- H Q1 + [y-(ct+3+l)x] = 0 . 

We now consider for H. a solution of the form 

Oc 

OO 

H = 53 d n x n+C where d * 0 . 
n=0 

Introducing for Hq^, the solution we found before, we have for the 
equation defining d^, two choices: 

1) c(c-l)d Q + y c d Q - a Q + y a Q = 0 and c=0 
where the terms in d^ compensate for the term in a^ or 

2) c(c-l)dg + y c dg = 0 and c < 1 

where the terms in d^ balance themselves. The first case is impossible 
as can be seen directly and we are thus left with 

c - 1 + y = 0 
or 

c = -m ' ■ 


Knowing that c = -m, we could now compute the recurrence relationship 
between the d^s. But/ once more, the solution is to be disregarded 
for its singularity at the origin, and from there, the solution found 
before in Chapter 3, represents the general solution of our equation. 
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APPENDIX III 


BEHAVIOR NEAR THE ORIGIN 

In order to shorten the text, the behavior of the solution near 
r=0 has been introduced as an Appendix. It is indeed deduced that 
the solution previously found is valid throughout the neighborhood of 
the origin and doesn't present any singularity. As explained earlier, 

the reason to consider the neighborhood of the origin is justified for 

2 - ■ 

the singularity of the operator at the origin. 

Boundary Layer Expansion Valid Near r=0 

Close to r=0, let us define a suitable boundary layer variable 


or 


* 

r = 


0(e) 


where 0(e) -*■ 0 when £ 0 


The corresponding asymptotic expansion valid near r=0 is taken as 
4> a ,m( r »e) = y Q ( £ ) g Q (r*) + y l (e) g x (r*)+... 


In the new variables, the operator becomes: 


2 ? ? 9 2 

v 2 _ _d_ 1 _d _ n" d 1_ _d_ _ m = JL * 

m , 2 + r dr 2 _ 2*2 2*,* 2 *2 2 m 

dr r 0 dr 0 r dr 0 r 0 

*2 


where V stands for 
m 


d 1 d 

*2 + * * 
dr r dr r 


2 

m 

*2 


The complete equation for the deformation becomes, with the use of the 
previous definition: 
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-e 


■4 C C (>'o (e)s 0 Cr * ) + > J l (e)g l <,: * )+ -") 

O r . 

*(l-k 2 a 2 r* 2 ) (ti 0 (e)g 0 (r*) + y 1 (e)g 1 (r*)+. . 

~^2~ [l-k 2 a 2 r* 2 ] ^y Q (e)g c '(r*) + ^(e) g 1 (r*)+. . . j 
Q + ..) lu 0 (e>gg(r ) + u 1 (e)g 1 (r )+...! = 0 . 


_1 d 

2 * * 
a r dr 


+ (A 


The thickness of the boundary layer is obtained by considering: 


e ~ 1 j * r 

—7 — — or a = /£ and r = 

a 4 a 2 


The dominant boundary layer equation is then obtained by considering 
the terms of order y Q (e)/e in the previous development or 

2 * 

*2 *2 * Id , * d 8 0^ r ' m * 

-v- v- ) + 4 -tt 8 0 (r > * — h — g o (r > - 0 ■ 

r dr dr 


'm m e 0' 


This expression can be written shortly: 

*2 *2 * *2 , *. 

- V V Bn< r > + V 8 n (r ) = 0 • 

m m U mu 


*2 


Writing this last expression in terms of the operator V m , we have 

*2 *2 , * 

V/ C - l)g n (r ) = 0 
m m U 

or 

*2 *2 * 

07 - 1) V g (r ) = 0 . 

m m u 

*2 *2 

We observe a permutativity between the. two operators (V^ - 1 ) and . 
As a result, we know from the theory of linear differential equation, 
that the obtention of the four linearly independent solutions of the 
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previous equation can be done by solving the two systems 

(V * 2 * * - 1) g (r*) =0 
m u 

*2 , 

V > ‘ 0 • 

The general solution of the complete equation is then obtained by 
taking a linear combination of the four solutions generated. The 
determination of the constants appearing in the general solution are 
then determined through the matching process near the origin. Let us 
look first at the solutions of 



(V* 2 - l)g Q (r ) = 0 or 


2 2 
a _1 d_ ( * _m_ 

*2 s 0 + * dr S 0'‘ r ) *2 

dr r r 


g Q (r *) 



This equation is a Bessel equation which solution is given in terms of 

the modified Bessel functions I and K . But both of these solutions 

n n 

have to be canceled, one for its exponential growth and the other for 
introducing an unbounded displacement at the origin. The other inde- 
pendent solutions are determined in solving 
*2 , 

7 m 2 g 0 (r)-0. 

This equation can easily be integrated by the use of the substitution 
* *k 

g^(r ) = r leading to the characteristic equation, with two roots or 

2 2 

k = m or k = ± m . 


The solution k =-m introduces a singularity at the origin and will 

therefore be neglected and we are left with 
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This expression represents precisely what the solution found in the 

central area becomes in the neighborhood of the origin and the 

determination of the constant Cq is obtained directly. 

Due to the fact that there is no boundary condition at the 

2 * 

origin and that both solutions found for (V - l)g.(r ) = 0 have to be 

m U 

canceled for being improper for matching, we deduce that the solution 
found earlier is valid in the neighborhood of the origin, and that, 
for this problem, there is no need for a boundary layer at the center 
of the membrane. 



* 

Again, the solution in Jin r has to be rejected for its singularity at 
the origin, and we are left with 
g 0 (r*) - C r 


which represents the displacement of the membrane at the origin. The 
determination of is again obtained by matching with the solution 
found previously. Once more, there is no need for a boundary layer in 
the neighborhood of the origin. We conclude then, by saying that the 
solution found earlier is valid in all the central area of our problem. 
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APPENDIX IV 


DETERMINATION OF THE CONSTANT C, 


In the matching process valid near r=l, it seems natural to con- 
sider as an order of magnitude for v^(e) the value 


V-jXe) = e 


1/3 


Going back to the expansion done for the differential equation defining 


4>, 


= 1/3 


The 


a m> an indicated choice for <^(e) is given by K^(e) = e 
differential equation defining h^(r) is then given by the following 
expression 

2 

2 ^ ^1 1 2 ^1 1 2 2 2 2 2 2 
(1-r ) 1 + ~ (l-3r^) [m "(mV + Apr^ =-A‘h Q . 

dr r 

We can see from a straightforward observation that the homogeneous part 

of the previous equation is the same as the differential equation 

defining h^, and that hg appears as a forcing terms. As a result, the 
2 

eigenfrequency A^ of the previous equation, can be obtained once more 
by considering only finite displacement at the rim of the disk. Let us 
express: 


h l = £ b n 


n+m 


n=0 


The index of the last equation is obviously the same as the one found 
for fig. We then have if we replace h^, by its development 


h o - £ V 


n+m 


n=0 


where all a^'s are known. 
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oo 


or 


^ v [ (n+m) (n-hn-l)b 
n=2 n 


(n+m-2) (n+m-3)b „+ (n+m)b 

n-2 n 


„ . 2, ,.2 2 .2,, , n+m-2 

- 3 (n+m-2) b „-m b + (k m + A„)b „]r 
n-2 n 0 n-2 

CO 

, 2 S" 7 ' n+m-2 

= " X 1 a n-2 r 


(n+m-2) (n+m)- (k^m^ + A^) 

b = 5 — b „ 

n 2 0 n-2 

n + 2mn 


n^+2mn 


n-2 


The recurrence relationship contains two terms, one has the same ratio 

as the one found for the a^'s and the second term is coming from the 

2 

forcing terms. But if A^ is such that the serie defined for h^ is 

truncated starting at a ,, we have that a =0. But in the 

n n +z 

expansion of h^, the term 


b n' 


! (n ' ^+2n * m) 


n’-2 


represents a nonzero contribution. 

As a result, we conclude that the rest of the expansion for h^ 

represents a diverging expression for r=l, as it was proven for h^. 

We, then have to conclude that k^(e) is not of an order of magnitude of 
1/3 

e , and by the orthogonality relationship between the terms of the 

1/3 

asymptotic expansion valid for the mode, that v^(e) is not e . This 
result justifies the choice 
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APPENDIX V 


EQUATIONS FOR THE ROTATIONAL MOTION 


In this appendix, we will rederive some basic equations whose 
consideration are necessary for the mentioned particular case. First 
of all, let us derive the equations of motion for the attitude angles 
of the central body, for a completely general problem. 

In the derivation done in Chapter 4, we obtained the following 
expression for the angular momentum vector with respect to the center 
of mass. 


h = I N • to 


c * c + / 

Jk 


Q_ x £ dm . 


The equation defining the angular velocity vector for a freely spinning 

flexible spacecraft is then obtained by expressing the nullity of the 

■" : "I ■ 

applied torque acting about the center of mass or 


• N . N N 

h = 0 = I • to + to x i . to + I • to + 


to x^ j " £ x £ dm + J " jd x dm 


where we limited our derivation to the linearized equation of the 
motion. Written in matrix form, the previous equation becomes 

„ _N • ~ _N °N t ~ j / ~~ 

0 = 1 to + u I w + I (0 + to / p p dm + / PP dm . 

"'a •'a 

The two last terms can be expanded in the following way: 


/ p p dm = / T u dm + / u u < 

J A J A 

f~ fi- f~ 

I pp dm= / ru dm+ fuu 
J h -'a •'k 


dm . 
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Expanding similarly the matrix of inertia of the complete system into 


N N 

I = i" + 2A + U 


where 2A represents the first order terms of the matrix of inertia and 
similarly, -U represents the second order terms. 

The linearized equation of the rotational motion are given by 


0 = 1^ a) + 2 i ss A o) ss + 2 A 0) ss + oi ss 


/ ? “ 

Jk 


dm 


+ fru 


dm 


(V.l) 


where has been introduced to express the angular velocity vector 
when the system is spinning at equilibrium at a rate compatible with 
the angular momentum. 
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